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Background  and  Objectives  of  the  Investigation 

Recent  geological  and  geodetic  studies  have  suggested  that  the  region  surrounding 
Vandenberg  AFB  is  undergoing  active  crustal  deformation,  with  important  implications  for 
both  the  geodetic  stability  and  the  seismogenic  potential  of  the  Western  Test  Range  (WTR) 
[Feigl  et  ai,  1990].  Part  of  the  evidence  for  significant  deformation  was  obtained  from 
GPS  measurements  over  a  broad  area  of  central  and  southern  California,  which  we  carried 
out  in  cooperation  with  other  university  and  government  scientists  from  1987  through 
1991.  Although  useful  in  defining  the  regional  tectonic  setting,  these  measurements  are  of 
insufficient  spatial  and  temporal  density  to  answer  many  important  questions  about  the 
seismogenic  potential  of  Vandenberg. 

In  1989  we  received  funding  from  AFOSR  under  the  Commander's  Reserve  Fund 
(Grant  AFOSR-89-04(X)),  with  matching  funds  from  MIT,  to  purchase  five  GPS  receivers 
and  to  begin  a  series  of  measurements  designed  to  determine  the  magnitude  and  spatial 
distribution  of  deformation  in  the  Santa  Maria  Fold  and  Thrust  Belt  (SMFTB),  a  region 
encompassing  the  major  faults  and  folds  within  50  km  of  Vandenberg.  During  the  course 
of  this  grant  we  have  carried  out  four  GPS  surveys  of  the  SMFTB.  Together  with  the 
earlier  measurements,  these  surveys  have  allowed  us  to  determine  the  relative  velocities  of 
stations  within  the  region  with  accuracies  of  a  few  millimeters  per  year.  Two  additional 
receivers  were  purchased  in  May,  1992,  and  have  since  been  installed  in  continuously 
operating  GPS  stations  at  Vandenberg  and  the  China  Lake  Naval  Weapons  Center.  They 
are  operated  as  part  of  the  Permanent  GPS  Geodetic  Array  (PGGA)  in  California, 
providing  the  ability  to  monitor  not  only  interseismic  deformation,  but  any  transient  motion 
which  might  occur  prior  to,  during,  or  after  an  earthquake.  Observations  from  the 
Vandenberg  station  allowed  us  to  observe,  for  the  first  time,  far-field  displacements  from 
the  June,  1992,  Landers  earthquake  [Bocketal.,  1993]. 

The  primary  scientific  goal  of  our  investigation  is  to  characterize  the  partition  of  strain 
among  the  faults  and  folds  in  the  vicinity  of  Vandenberg.  But  GPS  observations  in  the 
vicinity  of  Vandenberg  cannot  be  understood  in  isolation  from  the  rest  of  southern 
California.  Only  by  determining  (vector)  motions  of  the  stations  in  the  Vandenberg 
network  relative  to  stations  outside  the  network  can  we  assess  the  relative  importance  of 
dilatation  and  rotation  in  the  pattern  of  deformation.  Moreover,  it  is  essential  to  include  in 
the  analysis  observations  from  a  continental-  or  global-scale  network  in  order  to  reduce  the 
errors  due  to  unmodeled  motions  of  the  GPS  satellites.  Thus,  we  have  analyzed  all  of  the 
Vandenberg  data  simultaneously  with  both  GPS  and  VLBI  data  from  central  and  southern 
California  and  a  global  network  of  stations.  Our  approach  and  the  results  of  our  analysis 
for  data  acquired  through  1991  are  described  in  a  comprehensive  paper  published  in  the 
Journal  of  Geophysical  Research  [Feigl  et  al.,  1993]  and  included  as  a  supplement  to  this 
report.  We  have  now  completed  the  analysis  of  the  1992  and  1994  Vandenberg  surveys 
and  included  these  data  in  our  estimates  of  station  velocities. 
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The  tectonic  setting  and  geological  descriptions  for  the  Vandenberg  region  are  described 
in  Feigl  et  al  [1990],  Feigl  et  al.  [1993],  and  Cooray  [1994]  and  will  not  be  repeated  in  this 
report. 


GPS  Measurements 

Since  January  1987,  we  have  carried  out  nine  GPS  surveys  involving  the  stations  of 
the  Vandenberg  network  (Figure  1).  Four  of  these  surveys  (February  and  September 
1990,  March  1992,  and  February  1994)  included  most  or  all  of  the  current  network,  the 
others  three  or  four  stations  at  the  comers  and  sometimes  a  single  station  (GRAS)  in  the 
center  (See  Table  1).  Our  primary  concern  in  densifying  our  network  was  bridging  the 
major  tectonic  structures  of  the  SMFTB.  In  1990  we  added  two  stations  (LIND  and 
ALAM)  on  or  across  the  Orcutt  and  Point  San  Luis  anticlines.  In  1992  we  added  two 
stations  (FARF  and  RDRK)  that  span  the  Lompoc-Purisima  anticline  and  one  (FIGP)  at  the 
eastern  extreme  of  the  network  to  provide  a  stronger  tie  to  a  potentially  stable  region  outside 
of  the  SMFIB. 

A  subset  of  our  measurements  were  made  within  10  km  of  the  VLBVPGGA  site  on 
South  Vandenberg  in  order  to  establish  that  the  monuments  there  are  not  subject  to  local 
motion  due  to  soil  instability.  We  surveyed  these  stations  in  1990,  1992,  and  1994,  and 
NASA  surveyed  nearly  the  same  set  in  1991. 

Prior  to  our  1990  SMFTB  survey,  GPS  measurements  in  California  had  been 
performed  using  one  or  two  observers  per  instrument  with  attended  operation  over  6-8 
hours.  The  availability  of  relatively  secure  sites  allowed  to  reduce  demonstrate  a  much  less 
costly  mode  of  operation,  with  two  observers  managing  nearly  continuous  observations  by 
6  instruments.  Our  scheme  became  the  prototype  for  similar  surveys  in  other  parts  of 
California  and  the  world  during  the  past  three  years. 
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Table  1.  Vandenberg  Network  Surveys  1987-1994 


Stations 

1987 

1987 

1988 

1989 

1990 

1990 

1991 

1992 

1994 

Jan 

Oct 

Mar 

Mar 

Jan 

Sep 

Feb 

Mar 

Feb 

VNDN 

South  Vandenberg 

• 

• 

• 

• 

• 

• 

• 

• 

• 

LOSP 

Mount  Lospe 

• 

• 

• 

• 

• 

• 

• 

• 

• 

MADC 

Madre  Peak 

• 

• 

• 

• 

• 

• 

• 

• 

• 

GAVI 

Gavlota  Peak 

• 

• 

• 

• 

• 

• 

• 

GRAS 

Grassy 

• 

• 

• 

• 

• 

• 

ALAM 

Alamo 

• 

• 

• 

• 

LIND 

Linda 

• 

• 

• 

FARF 

Farfield 

• 

• 

RDRK 

Red  Rock 

• 

• 

FIGP 

Figueroa  Pine  Tree 

• 

• 

ALVA 

Alvado 

• 

• 

• 

• 

RUSl 

Rustad 

• 

• 

• 

• 

Figure  1.  Tectonic  map  of  Vandenberg  and  vicinity  (Santa  Maria  Fold  and  Thrust  Belt)  showing  the 
location  of  faults  (thick  lines)  and  folds  (thin  lines)  and  of  the  GPS  stations  (dots  with  4-letter  codes) 
included  in  our  study.  Solid  dots  indicate  stations  which  have  been  surveyed  both  in  1990  and  1992;  open 
dots  indicate  stations  surveyed  for  the  first  time  in  1992. 


Results 
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Our  GPS  estimates  of  motions  of  stations  in  the  SMFTB  network  and  surrounding  area 
are  shown  in  Figure  2.  The  dominant  motion  for  the  northernmost  sites  in  Figure  2  is  the 
simple  shear  due  to  the  San  Andreas  Fault  (SAF)  system,  which  contributes  a  few 
millimeters  per  year  even  to  sites  as  far  west  as  Vandenberg. 

In  order  to  understand  better  more  local  sources  of  deformation,  we  show  in  Figure  3 
residual  velocities  after  removing  a  model  of  velocities  associated  with  the  SAF  system 
{Feigl  et  al.,  1993;  Cooray,  1994].  With  the  component  of  motion  due  to  shear  on  the  San 
Andreas  removed,  the  pattern  of  deformation  shows  two  main  features:  There  is  regional 
compression  normal  to  the  San  Andreas  and  residual  shear  parallel  to  the  San  Andreas. 
The  largest  gradients  in  shearing  are  in  the  vicinity  of  Vandenberg  (VNDN).  Sites  Santa 
Rosa  Island  (SNRI)  and  Center  (CENT)  on  the  northern  Channel  Islands  show  no 
resolvable  motion  relative  to  VNDN,  while  sites  Gaviota  Peak  (GAVI)  and  La  Cumbre 
Peak  (LACU)  on  the  north  side  of  the  Santa  Barbara  Channel  show  motion  to  the  southeast 
relative  to  VNDN.  Southeasterly  residual  velocities  are  also  observed  at  sites  Piedras 
Blancas  (BLAN),  Black  Hill  (BLHL),  and  Lospe  (LOSP),  all  in  the  vicinity  of  the  Hosgri 
fault  zone.  The  most  straightforward  interpretation  of  these  southeasterly  velocities  relative 
to  VNDN  is  that  strike-slip  motion  on  the  Hosgri  fault  cuts  through  the  Santa  Maria  basin 
east  of  VNDN,  joining  the  zone  of  compression  extending  from  the  central  Santa  Barbara 
Channel  [Larsen  et  al,  1993]  through  the  eastern  Ventura  basin  [Donnellan  et  al,  1993]. 
There  is  no  clearly  recognized  geological  structure  to  accommodate  this  motion.  But  a 
similar  situation  existed  for  the  Eastern  Mojave  Shear  Zone  before  the  Landers  earthquake 
[e.g..  Savage  et  al,  1990] 

In  the  vicinity  of  VNDN,  a  number  of  stations  exhibit  residual  motions  which  are 
significant  at  95%  confidence.  The  largest  residual  velocities  are  at  Grassy  (GRAS)  and 
Madre  (MADC).  We  believe  that  the  Grassy  motion  is  anomalous,  possibly  related  to  oil 
withdrawal  or  an  unusual  earthquake  (see  below).  Madre’s  southwesterly  motion  relative 
to  station  Linda  (LIND),  3  ±  0.5  mm/yr,  indicates  significant  compression  in  the  San  Lucia 
and  Sierra  Madre  ranges.  The  motions  of  Lospe  (LOSP)  and  Alamo  (ALAM)  suggest  a 
small  amount  of  shortening  (1  ±  0.5  mm/yr)  on  the  anticlines  to  the  south. 

We  have  calculated  the  (one-sigma)  uncertainties  in  our  estimated  velocities  by  scaling 
the  formal  errors  from  our  analysis  by  a  factor  of  2.0  to  reflect  the  goodness-of-fit  (chi- 
square)  in  our  combined  analysis  (see  Feigl  et  al.  [1993]  for  a  complete  discussion).  In 
order  to  assess  the  reliability  of  the  velocities  for  the  SMFTB,  however,  it  is  necessary  to 
examine  the  time  evolution  of  relative  position  for  stations  in  the  network.  In  Figure  4,  we 
plot  time  series  for  four  key  interstation  vectors.  Figures  4a  and  4b  show  relative  positions 
between  South  Vandenberg,  Mt  Lospe,  and  Madre,  the  three  stations  with  the  longest 
series  of  observations.  For  three  of  the  four  horizontal  components,  the  normalized  root- 
mean-square  (nrms)  scatter  (i.e.,  the  square-root  of  chi-square  per  degree  of  freedom)  is 
higher  than  the  value  (2.0)  used  to  scale  the  uncertainties  from  our  analysis.  For  the  north 
components,  however,  the  high  nrms  is  due  largely  to  outliers  for  the  surveys  near  the 
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iiiddle  of  the  span  (1990),  which  have  little  influence  on  the  estimated  velocities.  For  the 
east  component  of  the  VNDN-LOSP  vector,  the  best-fit  straight  line 
the  endpoints,  and  the  nrms  is  less  than  2.0.  The  east  component  of  is  he 

most  wVrrisome  since  a  shift  of  either  the  earliest  (1987)  or  latest  (1994)  values  by  the 
amount  each  deviates  from  the  best-fit  line  would  change  the  velocity  estimate  by  severa 
times  the  scaled  uncertainty.  If  this  shift  were  10  mm  (a  reasonable  upper  teund  from  the 
plot),  the  velocity  change  would  be  about  1.5  mm/yr.  Nevertheless,  for  these  three 
stations,  measured  eight  or  more  times  over  seven  years,  we  conclude  that  the  relative 
motions  and  uncertainties  shown  in  Figures  2  and  3  are  reliable. 


The  estimated  motions  of  stations  in  the  center  of  the  network  are  more  problematical. 
There  are  six  measurements  between  Grassy  (GRAS)  and  Vandenberg  spanning  almost 
seven  years,  but  the  variation  of  estimated  position  is  large  and  somewhat  systematic 
(Figure  4c).  The  horizontal  components  of  the  interstation  vector  (baseline)  could  it 
equally  well  with  the  slopes  shown  (indicating  northwesterly  motion  of  ~3  mm/yr)  or  wi 
no  motion  but  a  discontinuity  between  September,  1990,  and  February,  1992.  The  vertic 
motion  is  also  systematic,  suggesting  uplift  of  20  mm/yr  between  October,  1987,  and 
September,  1990.  From  a  geodetic  point  of  view.  Grassy  is  unusual  because  we  measured 
its  position  (with  respect  to  Vandenberg)  only  in  "small"  surveys  involving  simultaneous 
observations  from  3  or  4  SMFTB  stations,  whereas  Madre  and  Lospe  were  usually  part  of 
regional  surveys  involving  12  or  more  stations.  If  we  were  to  remove  points  from  the  most 
suspect  of  these  small  surveys-October,  1987,  and  September,  1990-the  honzontal 
velocity  estimate  indicated  by  the  upward  sloping  lines  Figure  4c  would  seem  more 
reasonable.  Grassy  is  also  anomalous  from  a  geophysical  point  of  view.  First,  it  is 
situated  above  an  anticline  where  there  is  currently  pumping  of  oil,  as  well  as  injection  o 
water  (personal  communication,  R.  Hauser,  Calif.  DOGGR,  1994).  Second,  an  unusual 
"slow"  earthquake  occurred  on  the  Orcutt  anticline,  10  km  north  of  Grassy,  m  January, 
1991  [Kanomori  and  Hauksson,  1992].  While  the  expected  displacement  of  Grassy  from 
this  earthquake  is  too  small  to  explain  our  observations,  it  might  be  that  other,  subseismic, 
displacements  were  associated  with  this  event 

The  final  time  series  shown  in  Figure  4  is  for  the  vector  between  Vandenberg  and  Linda 
(LIND),  for  which  the  rate  of  change  has  been  estimated  from  only  three  measurements. 
Our  estimate  of  nearly  zero  shortening  between  these  two  stations  could  be  changed  by 
nearly  2  mm/yr  (in  either  direction)  by  removing  the  first  or  last  of  the  three  values.  From 
the  VNDN-LOSP  and  VNDN-MADC  plots  (Figures  4  a  and  b),  the  more  suspect  value  is 
the  first  one.  from  the  September,  1990,  survey.  If  the  north  components  of  position  for 
Lospe  Madre,  and  Linda  estimated  from  this  survey  are  all  low  by  about  7  mm,  then  Linda 
would  be  moving  south  by  2  mm/yr.  implying  that  shortening  across  the  SMFTB  is  more 
evenly  partitioned  north  and  south  of  Linda  rather  than  primarily  between  Lmda  and  Madre. 
We  cannot  rule  out  this  possibility  from  our  current  analysis. 


238*30'  239*00'  239*30'  240*00'  240*30'  241*00'  241*30'  242  00' 
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Figure  3.  Residual  velocities  with  respect  to  Vandenberg  (VNDN)  of  sites  in  the  SMFTB  and 
surrounding  area,  after  removal  of  a  model  of  shear  strain  due  to  the  San  Andreas  fault  system  [Cooray, 
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Figures  4a  and  4b  Observed  north,  east,  op,  and  length  components  of  the  baselines  between  South 
Vandenberg  and  Mt  Lospe  (in  the  sense  LOSP  -  VNDN)  Oeft)  and  between  South  Vandenberg.and  Ma^ 
Peak  (MADC  -  VNDN)  (right).  Each  value  represents  an  average  over  a  2-  to  6-day  set  of  measurements, 
each  8  to  24  hours  long.  The  uncertainties  shown  for  each  value  are  formal,  before  scaling,  and  thus 
should  be  scaled  by  a  factor  of  2.0  to  obtain  our  best  estimate  of  the  one-sigma  uncertainty. 
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Figures  4e  and  4d  Observed  north,  cast,  up,  and  length  components  of  the  baselines  between  South 
Vandenberg  and  Grassy  (in  the  sense  GRAS  -  VNDN)  (left)  and  between  South  Vandenberg  and  Linda 
(LIND  -  VNDN)  (right).  Each  value  represents  an  average  over  a  2-  to  6-day  set  of  measurements,  each  8 
to  24  hours  long.  The  uncertainties  shown  for  each  value  are  formal,  before  scaling,  and  thus  should  be 
scaled  by  a  factor  of  2.0  to  obtain  our  best  estimate  of  the  one-sigma  uncertainty. 
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Summary  and  Future  Work 

GPS  measurements  in  the  SMFTB  and  surrounding  area  reveal  a  zone  of  both  strike 
slip  and  compressive  deformation  not  readily  identifiable  with  known  faults.  Stimulated  by 
the  recent  Landers  (June,  1992)  and  Northridge  (January,  1994)  earthquakes,  the  Southern 
California  Earthquake  Center  (SCEC)  has  begun  a  concerted  study  of  the  potential  for  a 
damaging  seismic  event  in  areas  of  significant  deformation.  Although  the  SCEC  study  has 
focused  primarily  on  the  highly  populated  Los  Angeles  basin,  the  SMFTB  is  an  important 
region  for  continuing  study  because  of  its  similar  tectonic  setting.  In  particular,  "blind" 
thrust  faults  likely  serve  as  the  principal  architects  of  the  topography  of  both  regions. 
Annual  or  biennial  GPS  measurements  under  SCEC  auspices  will  allow  us  to  refine  our 
estimates  of  deformation  in  the  SMFTB  network. 

The  measurement  techniques,  data,  and  analysis  tools  developed  under  this  grant  are 
among  the  most  significant  contributions  the  use  of  GPS  in  tectonic  studies.  We  have 
documented  these  in  the  three  journal  papers,  three  MIT  Ph.D.  thesis,  and  one  MIT  Senior 
thesis  listed  below.  A  more  detailed  analysis  of  the  results  presented  in  this  report  is  being 
prepared  for  submission  to  the  Journal  of  Geophysical  Research.  We  will  forward  a  copy 
to  AFOSR  when  it  is  submitted. 
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Space  Geodetic  Measurement  of  Crustal  Deformation 
in  Central  and  Southern  California,  1984— 1992 

Kuril  Feigl'--  Duncan  C.  Agnew,-*  Yehuda  bock.-’  Danan  Dong,'-’  Andrea  Donnellan,’-'' 
Bradford  H.  Hager,'  Thomas  A.  Herring,'  David  D.  Jackson, Thomas  H.  Jordan,' 
Robert  W.  King,'  Shawn  Larsen,^.’  Kristine  M.  Larson,-’-"  Mark  H.  Murray.'-'' 
ZhengkangShen,’*  and  Fr-ank  H.  Webb’-’ 

Wc  cslimalc  ihe  velocity  field  in  central  and  southern  California  using  Global  Positioning  System  (GPS) 
obsciA-ations  from  Id86  to  l‘)')2  and  very  long  baseline  interferometry  (VLBI)  observations  from  l<f84  to 
1991.  Our  core  network  includes  12  GPS  sites  spaced  approximately  50  km  apart,  mostly  in  the  western 
Transverse  Ranges  and  the  coastal  Borderlands.  The  precision  and  accuracy  of  the  relative  horizontal 
velocities  estimated  for  these  core  stations  are  adequately  described  by  a  95%  confidence  ellipse  w  ith  a 
semiminor  axis  of  approximately  2  mm/yr  oriented  roughly  north-south,  and  a  semimajor  axis  of 
approximately  .3  mm/yr  oriented  e.-ist-we.si.  For  other  stations,  occupied  fewer  than  5  times,  or  occupied  during 
experiments  w  ith  poor  tracking  geometries,  the  uncertainty  is  larger.  These  uncertainties  are  calibrated  by 
analyzing  the  scatter  in  three  types  of  comparisons:  (1)  multiple  measurements  of  relative  position 
("repeatability”).  (2)  independent  velocity  estimates  from  separate  analyses  of  the  GPS  and  VLBI  data,  and  (.3) 
rales  of  change  in  baseline  length  estimated  from  the  joint  GPS-sVLBI  solution  and  from  a  comparison  ol  GPS 
w  ith  irilaleralion.  The  dominant  tectonic  signature  in  the  velocity  field  is  shear  deformation  asstKialcd  w  ith  the 
San  Andreas  and  Garlock  faults,  which  wc  model  as  resulting  from  slip  below  a  given  locking  depth,  _ 
Removing  Ihe  effects  of  this  simple  model  frtim  the  observed  velocity  field  reveals  residual  delormation  that  is 
not  attributable  to  the  San  Andreas  fault.  Baselines  spanning  the  ea.siern  Santa  Barbara  Channel,  the  Ventura 
basin,  Ihe  Los  .Angeles  basin,  and  the  Santa  Maria  Fold  and  Thrust  Bell  arc  shortening  at  rales  of  up  to  5  ±  1. 

5  ±  I,  5  ±  I,  and  2  ^  1  mm/yr,  respectively.  North  of  the  Big  Bend,  some  compre.ssion  normal  to  the  trace  of 
the  San  Andreas  fault  can  he  re.solvcd  on  both  sides  of  the  fault.  The  rales  of  rotation  about  vertical  axes  in  the 
residual  geodetic  velocity  field  differ  by  up  to  a  factor  of  2  from  those  inferred  from  paleomagnetic 
declinations.  Our  estimates  indicate  that  the  "San  Andreas  discrepancy”  can  be  resolved  lo  within  the  .3  mm/yr 
uncertainties  by  accounting  for  deformation  in  California  between  Vandenberg  (near  Point  Conception)  and 
Ihe  westernmost  Basin  and  Range.  Strain  accumulation  of  1-2  mm/yr  on  structures  offshore  of  Vandenberg  is 
also  allowed  by  the  uncertainties.  South  of  the  Transverse  Ranges,  the  deformation  budget  must  include  5 


Ceulo^icul  and  SeLsinological  Indicators 

Dclbrmalion  in  the  southern  Coast  Ranges  (SCR  in  Figure  1)  is 
characterized  primarily  by  strike-slip  motion  on  the  San  Andreas, 
Hosgri.  Rinconada.  and  other  parallel  faults  [e.g..  Dibblee.  1977], 
Away  from  the  San  Andreas  fault  (SAF),  the  deformation 
includes  a  comprcssional  component  oriented  perpendicular  to  the 
trace  of  the  SAF  as  evidenced  by  subparallel  fold  axes  [Foge, 
1966.  1981]  and  thrust  faulting  earthquake  focal  mechanisms 
[Dchlingcr  and  Holt,  1988],  The  rale  of  shortening  has  been  esti¬ 
mated  at  7-13  mm/yr  from  a  balanced  cross  .section  extending 
from  the  SAF  lo  an  offshore  point  west  of  the  Hosgri  fault 
[Namson  and  Davis,  1990].  On  the  other  (northeast)  side  of  the 
SAF.  there  is  also  evidence  of  compressional  strain,  notably  the 
anticlinal  structures  .associated  with  t>il  production  [Callaway, 
1971]  and  the  1983  Coalinga  earthquake  [.V/cm  and  King,  1984]. 

Farther  south,  the  part  of  the  Santa  Maria  Basin  lo  the  northeast 
vif  Point  Arguello  forms  a  tectonic  transition  zone  between 
probable  strike-slip  motion  on  the  San  Gregorio-Hosgri  fault 
system  [Hall.  1978,  1981],  and  compre.ssion  in  the  western 
Transverse  Ranges  (WTR)  and  Santa  Barbara  Channel  (SBC)  lo 
the  south  [Crouch  el  al.,  19,84;  Namson  and  Davis,  1990]. 

The  Santa  Barbara  Channel  is  undergoing  north-.south  shorten¬ 
ing.  as  indicated  by  earthquake  focal  mechanisms  [YcrJ-cs  and 
i.ee.  1979]  and  geological  investigations  of  ft'Iding  and  faulting 
[Yc-o/.v,  19,81.  19,83],  The  average  rate  of  shortening  is  2-9  mm/yr. 
c.'-timalcd  from  1.8  km  over  the  last  0. 2-1.0  m.y.  [Ych/.v.  1983]. 

There  is  substantial  deformation  of  Oualcrnary  structures  ac¬ 
commodating  convergence  across  the  Transverse  Ranges.  For 
example.  Namson  and  Davis  [19,88]  propo.se  an  average  conver- 


mm.'yr  between  the  offshore  islands  and  the  mainland. 


Introduction 

Determining  the  velocity  field  in  the  vicinity  of  the  Pacific- 
North  America  plate  boundary  in  central  and  southern  California 
(Figure  1)  is  a  long-standing  problem  in  tectonics.  While  most  of 
the  motion  between  these  plates  occurs  on  the  San  Andreas  fault, 
the  deformation  extends  for  a  substantial  distance  on  either  side 
of  this  structure.  Such  off-fault  deformation  is  evident  in  geologic 
.structures,  seismicity,  paleomagnetic  declinations,  and  geodetic 
networks.  Measuring  that  deformation  with  space  geodesy  is  the 
primary  objective  of  this  study,  which  seeks  to  quantify  the  veloc¬ 
ity  field  in  this  intercontinental  plate  boundary  zone. 
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Fig.  I.  Generalized  map  of  California  showing  traces  of  major  Quaternary  faults  (solid  lines)  and  Neogene  fold  axes  from 
Sicin  and  Yeats  [1989]  (dashed  lines).  The  GPS  stations  (triangles)  are  labeled  with  the  four-character  codes  listed  in  Table  1. 
The  tectonic  domains  (in  italics)  include  the  Santa  Maria  Fold  and  Thrust  Belt  (SMFTB),  the  southern  Coast  Ranges  (SCR),  the 
western  Transverse  Ranges  (WTR),  the  Ventura  Basin  (VB),  the  Santa  Barbara  Channel  (SBC),  the  southern  Borderlands 
(SBL),  and  the  Eastern  California  Shear  Zone  (ECSZ).  Major  faults  include  the  San  Andreas  (SAF),  the  Elsinore  (ELS),  San 
Jacinto  (SJC),  and  the  Garlock  (GAR).  Stations  not  labeled  here  are  labeled  in  Figures  5  and  6. 


gence  rate  of  18-27  mm/yr  over  the  past  2-3  m.y.  for  a  cross 
.section  between  Ventura  and  the  White  Wolf  fault  in  the  Great 
Valley.  Yeats  [1983]  proposes  accelerated  deformation  near 
Ventura  during  the  past  0.2  m.y.  with  22  mm/yr  occurring  acro.ss 
this  region  of  the  Ventura  basin  alone. 

In  the  region  between  Ventura  and  Los  Angeles,  T.  L  Davis  et 
al.  [1989]  find  convergence  rates  south  of  the  SAF  of  8-16 
mm/yr.  Along  the  frontal  fault  of  the  San  Gabriel  Mountains 
(between  VB  and  JPLl  in  Figure  1),  Weldon  and  Humphreys 
[19.86]  conclude  that  the  rate  of  convergence  is  negligible,  while 
Bird  and  Rosenstock  [1984]  propose  convergence  rales  of  the 
order  of  10  mm/yr. 

The  deformation  in  the  Borderlands  south  of  the  Santa  Barbara 
Channel  is  much  less  well  known.  The  overall  structural  grain  and 
most  earthquake  slip  vectors  trend  norihwe.st  to  southeast  [Bent 
and  Helmberger,  1991],  but  the  July  1986  Oceanside  earthquake 
[Hauksson  and  Jones,  19.88;  Pacheco  and  NdHlek.  1988]  and  the 
San  Clemente  lault  [/.egg  el  ai.  19.89]  do  not  fit  this  pattern. 

Although  the  San  Andreas  fault  is  generally  considered  to  be 
the  major  seismic  hazard  in  this  region  (indeed,  it  ruptured  in  a 
major  (M  =  7.8)  earthquake  in  1.857).  most  of  the  seismic  moment 
relca.se  in  this  century  has  not  been  associated  with  rupture  of  the 
SAF  [Ellsworth,  1990],  This  is  apparent  in  the  map  of  all  events 


of  magnitude  (M^)  greater  than  4.0  in  the  past  60  years  (e.g.. 
Figure  4  of  Hutton  et  ah  [1991]).  Between  1812  and  1987,  four 
large  earthquakes  contribute  70%  of  the  moment  release  to 
produce  an  average  slip  rate  of  17-20  mm/yr  [Ekstriim  and 
England.  1989].  Of  these  four  events,  two,  the  1872  Owens 
Valley  (A/,v  =  7.5-7. 7)  and  the  1952  Kern  County  =  7.3-7.5) 
carthquake.s.  were  not  associated  with  the  San  Andreas  fault 
[Ellsworth,  1990].  The  recent  Landers  (My,-  =  7.3)  earthquake 
[Sieh  et  al.,  1993]  provides  another  example  of  substantial 
.seismic  moment  release  not  a.ssocialcd  with  the  SAF. 

Large  (approximately  30°)  palcomagnetic  declinations  ob- 
.served  in  po.st-Miocene  rocks  in  the  western  Transver.se  Ranges 
sugge.st  clockwise  rotations  about  vertical  axes  [Euyendyk,  1991], 
When  averaged  over  the  last  5  to  20  m.y..  these  data  imply 
rotation  rates  of  the  order  of  67m. y.  (0.1  urad/yr).  which  can  be 
explained  by  models  that  accommodate  simple  shear  on  rotating 
blocks  [Euyendyk,  1991]  or  by  “booksheir'  or  ■'collapsing  ladder" 
faulting  [y«c/:,w/;  and  Molnur,  199(1],  Although  these  models  in¬ 
voke  finite  rotations  over  several  million  years,  the  instantaneous 
rale  of  rotation  appears  to  he  of  the  same  order  of  magnitude 
[Jackson  and  Molnar.  1990]  and  can  be  resi'lved  by  our  geodetic 
network. 
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The  San  Andreas  Discrepancy 

While  ihere  is  abundant  evidence  for  regional  deformation  off 
the  SAF,  the  most  quantitative  argument  is  the  deviation  from  the 
vector  equality  required  by  rigid  plate  tectonics  [Minster  and 
Jordan,  1984].  The  starting  point  is  the  vector  velocity  of  the 
Pacillc  plate  with  respect  to  the  North  America. plate,  which  rep¬ 
resents  the  integrated  rate  of  deformation  across  the  “wide,  soft 
boundary”  [Atwater,  1970]  separating  the  two  tectonic  plates. 
From  this  vector,  one  subtracts  a  vector  for  the  direction-and"  rate 
of  slip  on  the  SAF  in  central  California.  The  resulting  nonzero 
vector  has  been  named  the  “San  Andreas  discrepancy”  [Minster 
and  Jordan,  1984].  A  large  part  of  this  residual  motion  may  be 
explained  by  extension  in  the  Basin  and  Range  Province  [A///is/er 
and  Jordan.  1987],  estimated  from  very  long  baseline  interfer¬ 
ometry  (VLBl)  to  be  9  i  1  mm/yr  at  N29  *  4°W  [Argus  and 
Cordon,  1991].  Subtracting  this  value  from  the  “discrepancy” 
yields  a  remainder  (the  “modified  discrepancy”)  estimated  to  be 
6*2  mm/yr  at  N20  ±  17°W  [.Argus  and  Cordon,  1991].  This 
motion  must  be  accommodated  to  the  west  of  the  Basin  and 
Range,  most  likely  in  onshore  California  [Sauber,  1988;  Argus 
and  Cordon,  1990;  Suvoge  <?/ a/..  1990;  Ward,  1990]. 

The  di.screpancy  is  reduced  by  an  additional  2  mm/yr  when 
recent  revisions  in  the  magnetic  time  scale  are  taken  into  account. 
The  current  plate  motion  model,  NUVEL-1  [DeMeis  el  al..  1990], 
averages  velocities  over  the  time  since  the  normally  magnetized 
chron  in  the  2A  sequence  of  palcomagnctic  anomalies,  taken  to  be 
at  3.03  Ma  in  the  time  scale  of  Harland  et  al.  [1982].  Recent 
revisions  in  the  time  scale  suggest  that  this  chron  is  several 
percent  older  [Hilgren.  1991;  McDougall  et  al.,  1992;  Cande  and 
Kent,  1992].  A  reasonable  approach  is  to  u.se  an  average  value  of 
3.16  Ma  from  these  three  studies,  which  means  that  the  NUVEL- 
1  rates  must  be  scaled  by  a  factor  of  0.959  (C.  DeMets,  personal 
communication,  1993).  Throughout  the  rest  of  this  paper,  we  shall 
use  these  slower  rates  and  refer  to  them  as  “rescaled  NUVEL-1.” 

Previous  Ceodetic  Studies 

The  triangulation  netw'orks  established  along  the  California 
coast  in  the  late  1800s  form  the  basis  of  much  of  the  subsequent 
geodetic  work  [Bowie,  1924;  1928].  Resurx'cys  led  to  the  eventual 
detection  of  a  “slow  drift”  of  roughly  50  mm/yr  of  the  crust  west 
of  the  SAF  [Whitten,  1956].  The  historical  survey  data  have  been 
used  to  infer  dextral  angular  (engineering)  shear  .strain  rates  of  at 
least  0.2  prad/yr  for  most  of  the  western  Transverse  Ranges  [5;w_>’ 
el  al.,  1983,  1986,  1987].  Many  ofthc.se  historical  stations  have 
been  reoccupied  with  the  Global  Postioning  System  (GPS)  as  part 
of  our  study. 

Trilatcration  networks  monitored  by  the  U.S.  Geological 
Survey  (USGS)  straddle  the  major  segments  of  the  SAF  sy.stem. 
As  summarized  by  l.isowski  et  al.  [1991],  the  velocity  field  de¬ 
termined  by  trilatcration  in  central  and  .southern  California  is 
dominated  by  right-lateral  shear  a.ssociatcd  with  the  SAF  sy.stem. 
Indeed,  they  can  explain  all  the  deformation  observed  in  their 
networks  by  shear  alone,  requiring  remarkably  little  dilatation. 
Observed  changes  in  line  lengths  in  the  USGS  netwrnks  have 
been  modeled  as  due  to  strike-slip  motion  on  dislocations  buried 
at  depths  of  tens  of  kilometers  beneath  the  major  faults.  Inferred 
displacement  rales  for  the  SAF  north  of  the  Big  Bend  are  32- 
36  mm/yr  [Eherharl-l’hillips  et  al..  1960;  l.isowski  et  al..  1991], 
consistent  with  the  geological  inference  of  34  mm/yr  [Sieh  and 
Johns.  1984],  Given  the  success  of  these  models,  the  shear  com¬ 


ponent  of  the  SAF  motion  in  central  California  may  be  regarded 
as  a  relatively  well-understood  signal. 

Using  VLBl.  Clark  el  al.  [1987]  estimated  that  the  velocity  of 
Vandenberg  with  respect  to  Mojave  was  some  7  mm/yr  slower 
than  the  48  mm/yr  predicted  by  the  NUVEL-1  plate  motion 
model  [DeMets  et  al.,  1990].  This  result  implied  deformation  off¬ 
shore  of  Vandenberg  or  east  of  Mojave.  By  excluding  Mojave 
from  the  stable  North  America  plate.  Ryan  et  al.  [1993]  estimated 
the  motion  of  Vandenberg  relative  to  the  North  American  plate  to 
be  within  1  mm/yr  of  the  value  predicted  by  the  rescaled 
NUVEL-1.  VLBl  data  also  allowed  Ward  [1990]  to  attribute  ap¬ 
proximately  6  mm/yr  to  deformation  in  the  Coast  Ranges,  west  of 
the  SAF.  Since  this  motion  is  northward  with  respect  to  North 
America,  it  would  appear  to  include  a  fault-normal  component  as 
well  as  right-lateral  shear. 

The  southern  Coast  Ranges  are  actively  deforming,  as  has  been 
observed  by  triangulation  [^oi'age  and  Burford,  1973],  trilatera- 
tion  [King  et  al.,  1987],  and  a  comparison  of  triangulation  and 
GPS  [Shen,  1991;  .SVie/i  and  Jackson,  1993],  Most  of  this  motion 
appears  to  be  right-lateral  shear  describable  by  a  simple  model  of 
a  dislocation  in  the  plane  of  the  fault  [c.g..  l.isowski  ct  al..  1991], 
Other  models,  allowing  motion  normal  to  the  fault  [Harris  and 
Segall,  1987;  Segall  and  Matthews.  1988],  additional  tectonic 
blcxrks  [Cheng  et  al.,  1987]  or  time  dependence  [Z.i  and  Rice, 
1987],  suggest  that  the  deformation  away  from  the  San  Andreas  is 
more  complex  than  simple  shear  from  constant  slip  on  a  buried 
«:rew  dislocation. 

Another  deviation  from  simple  shear  occurs  in  the  western 
Transverse  Ranges,  where  several  geodetic  studies  have  measured 
active  compression.  In  the  area  south  of  Santa  Maria,  a  significant 
component  of  compression  has  been  estimated  from  a  comparison 
of  historical  triangulation  surveys  [Bowie.  1924,  1928;  Savage 
and  Prescott,  1978]  with  GPS  obserxations  [Feigl  et  al.,  1990; ' 
Shen  and  Jackson,  1993].  The  Santa  Barbara  Channel  is  also 
actively  shortening,  as  shown  by  comparisons  of  GPS  data  with 
historical  triangulation  measurements  [Webb,  1991]  and  1970s 
trilatcration  [l.ar.sen,  1991;  Larsen  et  al.,  1993].  In  the  Ventura 
basin,  comparison  of  GPS  observations  in  1987  with  triangulation 
from  the  1950s  has  measured  rapid  shortening  [Donnellan,  1992; 
Donnellan  el  al.,  1993a],  but  at  rates  less  than  half  the  geological 
estimates  for  the  last  0.2  m.y.  [Yeats,  1981,  1983]. 

Our  Ceodetic  Studies 

GPS  and  VLBl  observations,  especially  when  combined,  offer 
several  advantages  over  classical  terrestrial  techniques.  First,  they 
retain  high  precision  (less  than  10“”)  over  lines  longer  than  30 
km.  Second,  they  measure  the  Cartesian  vector  between  two 
stations  rather  than  the  distance  or  direction  only.  Third,  they 
measure  with  respect  to  a  single  reference  frame.  Taken  together, 
the.se  improvements  in  geodetic  technique  allow'  us  to  estimate  a 
precise,  self-consi.stent  relative  velocity  field  over  most  of  central 
and  .southern  California. 

The  southwest  United  Stares  has  been  the  ideal  location  to 
study  tectonic  motions  with  CPS  measurements  [Di.ron,  1991; 
Hager  et  al..  1961  ]  becau.se  the  Department  of  Defense  optimized 
the  initial  (Block  I)  satellite  constellation  to  provide  the  best  accu¬ 
racy  for  testing  in  this  region.  Measurements  made  in  California 
as  early  as  June  1986  have  shown  both  short-term  repeatability 
and  agreement  with  VLBl  for  horizontal  coordinates  at  the  sub¬ 
centimeter  level  for  intersite  distances  up  to  400  km  [Dong  and 
Bock.  1989;  Blewiil.  (989;  ./.  /,.  Davis  el  al..  1989;  Di.xon  et  al.. 
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1990;  Larson  and  A^ncw,  1991;  Murray,  1991],  This  accuracy 
implies  lhat  we  can  dclcrminc  velocities  at  the  level  of  a  few 
millimeters  per  year  in  5  years  with  annual  measurements.  The 
sites  for  which  w'c  determine  velocities  using  repeated  GPS  and 
VLSI  are  shown  in  Figure  1.  The  following  section  describes  the 
techniques  for  collecting  the  various  types  of  data;  VLBI  from  the 
global  network,  GPS  in  the  California  region,  and  GPS  al  globally 
distributed  stations. 

Data 

Between  1984  and  1992,  over  1700  VLBI  experiments  were 
performed  using  a  global  array  of  radio  telescopes  [Clark  et  ai, 
1985]  under  the  sponsorship  of  the  National  Aeronautics  and 
Space  Administration  (NASA)  [Coates  et  ai,  1985]  and  the 
National  Oceanic  and  Atmospheric  Administration  (NOAA) 
[Carter  et  al.,  1985].  Our  analysis  includes  1618  of  these  exper¬ 
iments.  dc.scribed  by  Ryan  et  al.  [1993j.  About  170  of  the.se 
experiments  have  included  at  least  one  fixed  and  one  mobile 
antenna  in  central  or  .southern  California. 

Most  of  the  GPS  observations  u.scd  in  our  analysis  were 
obtained  in  over  20  experiments  between  1986  and  1992.  Twelve, 
stations  make  up  the  “core"  network  extending  along  the 
California  margin  from  San  Simeon  to  San  Clemente  Island 
(Figure  1).  These  core  sites,  listed  in  Table  1,  were  all  occupied  at 
least  live  limes  in  the  5  years.  During  the  six  experiments  involv¬ 
ing  the  core  network,  GPS  receivers  also  occupied  three  or  more 
VLBI  sites  within  California  to  provide  a  regional  anchor  for  the 
network.  In  addition  to  the  core  experiments,  we  also  conducted 
more  than  a  dozen  smaller  experiments  to  increase  the  spatial 
density  of  the  network  in  regions  of  tectonic  interest.  These  small 
experiments  often  included  sites  also  measured  by  VLBI,  trilat- 
eration,  or  historical  iriangulation.  A  complete  list  of  .stations  with 
identifying  codes  and  approximate  positions  is  given  in  Table  1. 
The  field  experiments  are  listed  in  Table  2.  The  configurations  of 
the  tracking  sites  available  for  each  experiment  are  listed  in  Table 
3. 

Between  1986  and  1989,  all  of  the  field  observations  were 
made  using  Texas  Instruments  (TI)  4100  receivers;  between  1990 
and  1992,  most  were  made  usingTrimble  4000  SST  receivers.  As 
part  of  our  March  1990  campaign,  we  occupied  .seven  stations 
with  both  TI  4100  and  Trimble  receivers  in  two  successive  4-day 
experiments.  The  receivers  used  at  the  permanent  tracking 


stations  were  more  diverse  and  changed  over  lime  (Table  3).  They 
included  TI  4100,  MiniMac  2S16AT,  Trimble  4000  SST,  and 
Rogue  SNR-8  receivers. 

Altogether,  the  GPS  data  set  includes  useable  observations  for 
over  100  days  between  June  1986  and  May  1992.  Subsets  of  these 
data  arc'  dc.scribed  in  detail  in  several  preliminary  analyses. 
Experiment  0  is  described  by  Blewitt  [1989]  and  Dixon  et  ai 
[1990],  and  experiment  3  by  Dong  and  Bock  [1989].  The  first  2.7 
years  of  data  in  experiments  0.  2,  7,  8,  11.  and  14  are  analyzed  by 
/.orso/i  [1990a. b].  This. analysis  is  extended  to  include  a  1991  oc¬ 
cupation  of  .several  .siie.s  in  the  Channel  Islands  (experiment  SBl) 
by  Larson  and  Webb  [1992]  and  Larson  [this  issue].  Details  on 
the  occupations  of  specific  sites  are  given  by  Murray  [1991]  for 
the  core  sites,  by  Donnellan  [1992]  and  Donnellan  el  al.  [1993a] 
for  the  Ventura  basin,  by  /-V/g/  [1991]  for  Vandenberg,  by  Shen 
[1991]  for  the  Coast  Ranges,  and  by  Larson  [1990a]  and  Larsen 
[1991]  for  the  Channel  Islands. 

To  improve  the  accuracy  of  the  coordinates  of  the  global  GPS 
tracking  sites  used  to  analyze  the  California  campaigns,  we  al.so 
include  data  from  two  additional  CPS  data  sets.  A  global  GPS 
campaign  was  conducted  for  23  days  in  January-February  1991 
under  the ‘auspices  of  the  IniernaTional  Earth  Rotation  Service 
(lERS)  with  the  coordination  of  the  Jet  Propulsion  Laboratory 
(JPL).  We  analyze  the  data  obtained  from  the  21  Rogue  receivers, 
which  provide  a  sufficiently  strong  global  network  to  improve 
significantly  the  coordinates  of  six  of  the  sites  (Algonquin, 
Tromso,  Wcttzell,  Kokce,  Scripps,  Pinyon,  and  JPL)  used  in  our 
California  experiments  [Herring  et  ai,  1991 ;  Blewitt  et  ai,  1992]. 
In  addition,  we  have  included  1 19  days  of  data  obtained  between 
October  1991  and  May  1992  from  35  global  sites  as  part  of  the 
operations  of  the  Permanent  GPS  Geodetic  Array  (PCGA)  in 
California  [Bock  and  Shimoda.  1990;  Bock,  1991;  Lindqwister  et 
al.,  1991;  Blewitt  et  al.,  1993;  Bock  et  ai,  1993].  Four  of  the 
global  sites  (Algonquin.  Tromso,  Wetlzell,  and  Kokee)  and  three 
PGGA  sites  (Scripps,  Pinyon,  and  JPL)  were  also  ob.served  during 
several  occupations  of  the  core  network. 

Data  Analysis 

We  analyze  the  GPS  and  VLBI  obserx'aiions  in  two  steps.  In 
the  first  .step,  we  perform  separate  least  squares  analy.ses  of  the 
GPS  phase  and  VLBI  group  delay  data  in  each  individual  day 
(“session’).  In  these  single-.se.ssion  solutions,  we  e.siimate  the  sta- 


TABLE  I.  List  of  Stations 


ID 

PID” 

Stamping  or  Localion* 

Lat°N 

Dw'W 

Height^  m 

ALAM 

DZ1334 

ALAMO  1^25 

54.79S5 

120.2568 

460 

ALVA 

DZI732 

ALVADO  I<333 

34.5927 

120.6170 

290 

BLAN"' 

FV1009 

NAVY  DEPT  12  NAVAL  DISTRICT 

35.6646 

121.2845 

•25 

BLHL'' 

BLACK  HILL  1S81  (R.M) 

35.35S7 

I20.S3I7 

170 

BLUF'' 

TZI974 

BLUFF  1933.  Sun  CIcmcnic  Island 

32.9268 

n  8.5  1 85 

300 

BOLD 

TZ1946 

BOULDER  1933 

32.8958 

1 18.4682 

560 

BOUC 

DX50SI 

BOUCHER  2  1975 

20.3347 

1  16.9193 

1660 

BPA3 

BP  ARIES  3.  Owens  Valley  Radio  Ohs 

37.2320 

It  8.2836 

1180 

BRSH 

DY2I50 

BRUSH  Ifi76,  Catalina  Island 

33.4070 

I  18.4049 

450 

CATO 

Castro  Peak 

34.0858 

1  18.7858 

830 

CATW 

Catalina  Island  West 

33.4598 

1 18.5687 

510 

CENT-' 

DY3159 

CENTER  1934,  Santa  Cruz  l.sland 

33.9948 

1  19.7529 

390 

CHAF 

EW786I 

CHAFFEE  2  1923  1941 

34.3006 

It  9.33  to 

310 

COTR 

EW6129 

Colar.  Pt.  Mugu 

34.1202 

11 9. 1540 

■30 

CSTL 

FVI421 

Castle  Mount 

35.9389 

120.3403 

1330 

DEVL 

ENH'KOVO 

DEVILS  PEAK  2  1951,  Santa  Cruz  Island 

34.0291 

11  9.7844 

7(X) 

ECHO 

EW7224 

ECHO  ROCK  C.  A. 

34.2249 

1  18.0550 

1720 

ELMO 

EW7230 

MERCED  R.  S.  1923  NO.  4  1972 

34.0302 

1  18.0951 

ISO 

FIBR'' 

FU1972 

A  364  1953,  Buiionwillow 

35.3985 

1  19,3940 

50 

FTOR 

NASA  CSFC  7266.  Ft.  Ord 

36.6698 

121.7733 

20 
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TABLE  1.  (continued) 


-  (D 

PID 

Slampinn  or  Localu>n 

Lal*N 

L(>n*'W 

Hciuhl.  m 

FORO 

NASA  GSFC  7421  1  WO,  Ft  Ord 

36.5894 

121.7721 

250 

FOR2 

NASA  GSFC  7421  RM2,  Ft.  Ord 

36.5897 

121.7716 

250 

CAVI 

DZ1256 

Cavioia 

34.5018 

120.1988 

710 

GRAS 

DZ1327 

GRASSY  uses  1959 

34.7306 

120.4141 

330 

HAPY 

HAPPY  1959 

34.3580 

1  18,8501 

670 

HAP2 

HAPPY2  1992 

34.3280 

1  18.8771 

350 

HOPP 

HOPPER  1941 

34.4777. 

-  118.8655 

1340 

JPLM 

Mesa.  JPL,  Pa.sadena  (PCGA) 

34.2048 

118.1732 

420 

JPLl 

EW1949 

JPLl  ARIES  1  1975,  Pasadena 

34.2047 

118.1710 

■440 

LACU-' 

EW8022 

La  Cumbre«  Santa  Barbara 

34.4944 

119.7139 

1160 

LIND 

BUREAU  OF  RECLAM  UNDA  1955 

34.9599 

120.2997 

ISO 

LOSP'' 

DZI559 

Ml  Lospe,  Vandcnbcfg  AFB 

34.8937 

120.6062 

460 

LOVE 

LOMA  VERDE  RESET  1961 

34.4963 

1  18.6687 

730 

MADCy 

MADRE  ECC1980 

35.0756 

120.0671 

960 

MILL 

DZ126I 

MILLER  1974 

34.5101 

120.2297 

90 

MOJl 

FTI572 

NCMN  1983  RM  1,  Mojave  Station 

35.3316 

1 16.8908 

930 

MOJA'' 

TI  4100  phase  center,  Mojave  Station 

35.3316 

1 16.8882 

900 

MOJF  - 

• 

FRPA“1  phase  center,  Mojave  Station 

35.3316 

116.S882 

900 

MOJM 

MiniMac  phase  center,  Mojave  Station 

35.3316 

116.8882 

900 

MONU 

DC1438 

MONUMENT  PEAK  NCMN  19S3 

32.8918 

116.4228 

1840 

MPNS 

MT  PINOS  USCiGS1941 

34.8128 

1 19.1454 

2660 

MUNS 

MUNSON  (TPC)  (USCE)  1971 

34.6358 

-  119.3006 

2110 

NIGU 

DX5266 

NIGUEL  A  1884  1981  - 

-33.5145 

1  17.7303.. 

240 

OCOT 

DB1234 

OCOTILLO  NCMN  1982 

32.7901 

1  15.7952 

0 

OVRO 

MOBLAS  7114  1979.  Owens  Valley  R.  O. 

37.2326 

1  18.2938 

1180 

PARG 

DZ1175 

PTARGUELLO  1933 

34.5549 

120.6160 

-10 

PEAR 

Pcarhlo-ssum  NCMN  WS3 

34.5121 

1  17.9224 

890 

PINY 

DX3617 

PINYON  FLAT  NCMN  1981 

36.6092 

1  16.4588 

1270 

PINl 

Pinyon  I  PGGA 

33.6122 

116.4582 

1260 

P1N2 

Pinyon  2  PGGA 

33.6121 

116.4576 

1260 

PL9A 

EW7395 

PICO  L9  A  1957 

34.3295 

1 18.6007 

1100 

POZE 

FV0810 

K  66  1927 

35.3474 

120.2955 

770 

POZO 

FV081 1 

LS6I  1957 

35.3460 

120.2987 

730 

PTDU 

EW4215 

POINT  DUME  RESET  1947 

34.0016 

118.8067 

60 

PVER'' 

PALOS  VERDES  ARIES  1976  1980 

33.7438 

118.4036 

70 

ROKY 

FV1829‘^ 

Rocky  Butte  2,  RM  1 

35.6653 

121.0596 

1050 

RUSI 

DZ1778«' 

RUSTAD  1933  RMl 

34.5708 

120.6270 

180 

SAFE 

PICOL9C 

34.3304 

118.6013 

1100 

SBA2 

EW7997 

SANTA  BARBARA  2  1956.  S.  Barbara  Is. 

34.4041 

119.7160 

140 

SBIS 

DY3066 

SANTA  BARBARA  ISD  2  1940 

33.4721 

119.0413 

160 

SC  LA 

SANTA  CLARA  1898 

34.32S7 

119.0392 

660 

SCRE 

EW8055 

Santa  Cruz  East 

34.0547 

119.5647 

60 

SCRW 

EW8085' 

Santa  Cruz  West  2.  RM  ! 

34.0732 

1 19.9180 

ISO 

SIOI 

DC212I 

Scripps  I  PGGA 

32.8678 

1  17.2523 

10 

5102 

Scripps  2  PCGA 

32.8675 

1  17.2524 

10 

SIVP 

Sierra  Vista  Park 

34.0660 

118.0120 

90 

SJOS 

FV1440 

SAN  JOSE  1884  1956 

35.3152 

120/2696 

1150 

SJUA 

DX4280 

San  Juan  (1886) 

33.9138 

1  17.7381 

540 

SLUl 

FV1464 

San  Luis 

35.2778 

120.5618 

870 

SMIG 

DZ1512'' 

NEW  SAN  MIGUEL  RM  2  1934 

34.0396 

120.3866 

210 

SNPA 

EW7538 

SANTA  PAULA  NCMN  1981 

34.3879 

1 18.9988 

ISO 

SNP2 

SANTA  PAULA  1941 

34.4404 

1 1 9.0096 

1480 

SNRI 

DZ1207 

SOLED  AD  1872  1^34.  Santa  Rosa  Island 

33.9509 

120.1057 

440 

SNTZ 

LA  COUNTY  COVINA  O  RM  NO  1 

34.0125 

117.8837 

360 

SOU 

EW7886 

SOLIMAR  1974 

34.2983 

119.3427 

-10 

SOU 

DC  1849 

Mt  Solcdad,  La  Jolla 

32.8399 

1 17.2525 

220 

SYNZ 

SANTA  YNEZ  1 1917  1990  CSC  DET  1 

34.5305 

1 19.9860 

1220 

TEPW 

Tepusquel  Witness 

34.9100 

120.1867 

P50 

TWIN 

DY2177 

TWIN  ‘564,  San  Nicholas  Island 

33.2318 

1  19.4790 

200 

VAND 

VLBI  STA  7223,  Vandenberg  AFB 

34.5561 

120.6164 

-10 

VNDN"' 

VLBI  STA  7223  RM  I  1983  DET  1  GSS 

34.5563 

120.6162 

-10 

VSLR 

TLRS  STA  7880  Vandenberg  AFB 

34.5560 

120.6164 

-10 

VICE 

DYIOII 

E  788  1946 

33.7419 

1 18.4107 

40 

WHIT 

Whitaker  Peak 

34.5574 

1  18.7428 

1220 

WHT3 

Whitaker  Peak 

34.5675 

118.7427 

1220 

WORK 

DY023(F 

WORKMAN  HILL  RESET  1978 

33.9917 

1 18.0029 

420 

YAM2 

uses  ELEV  2749  FT 

34.8525 

1  19.4844 

810 

YUMA 

Yuma  (Arizona)  NCMN  1^83 

32.9391 

1  14.2031 

238 

^  PID  is  the  “Permanent  IDcniincation*’  number  assigned  by  the  Nalii>nal  Geodetic  Survey. 
^  Stampings  arc  listed  in  uppercase;  lower  or  mixed  case  gives  ItKalion  or  description. 

^  Coordinates  are  geodetic  with  respect  to  the  NAD  83  ellipsoid. 

^  Core  site. 

A  reference  mark  (RM)  has  been  u.scd  and  the  PID  refers  to  the  main  monument. 
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TABLE  2.  List  of  Experiments 


_  Exp.  . 

Dales 

Days 

Revrs 

Area  or  Objective 

California  Sites  Observed 

0 

Jun.  1986 

4  ■ 

VLBI  sites  and 

BLUF 

BOLD 

BOUC 

CATW 

HATC 

LAJO 

MACA 

Channel  Islands 

MOJI 

MONU 

NIGU 

OTAY 

PINY 

PVER 

SANT 

sou 

TWIN 

VNDN 

YUMA 

1 

Dec.  1 986 

5 

9 

Greater  Lx)s 

BLHL 

FTOR 

LOSP 

MOJA" 

OVRO 

POZO 

POZE 

•  ■  -  •'  ' 

'  --  • 

Angeles  Basin 

PVER 

SNPA 

WDN 

2 

Jan.  1987 

5 

12  . 

core  network 

BLAN 

BLHL 

CENT 

COTR 

FIBR 

FTOR 

LACU 

LOSP 

MADC 

MOJA 

OVRO 

VNDN 

3 

Jan.  1987 

5 

12 

Hi.storical  and 

BRSH 

CHAF 

DE\T 

FIBR 

FTOR 

GAVI 

LACU 

core  silcs 

MILL 

MOJA 

OVRO 

PVER 

SNRI 

SOU 

VNDN 

•  .  - 

-  ■  ' 

VSLR 

- 

- 

4,  5,  6 

Jan.  1987 

VLBI  .sites 

BLAC 

COTR' 

lActr 

MONU" 

NIGU" 

OCOT” 

PEAR" 

PINY" 

PVER* 

SAND" 

SOU" 

YUMA" 

7 

May  1987 

4 

6 

central  network 

BLHL 

(BTF 

MOJA 

OVRO 

PVER 

VNDN 

8 

Sep.  1987 

4 

7 

central  network 

BLHL 

CENT 

FTOR 

MOJA 

OVRO 

PVER 

SCRE  ■ 

SCRW 

VNDN 

9 

Sep.  1987 

7 

historical  sites 

CHAF“ 

CSTL' 

DEVL" 

ECHO" 

ELMO" 

GRAS 

MOJA 

NIGB“ 

PVER 

SBA2" 

SJUA" 

SLUl" 

SOU" 

TEPT 

TEPW 

VNDN 

• 

Los  Angeles  and 

CATO 

DELT 

ECHO" 

ELMA" 

ELMO" 

HAPY 

HOPP 

Ventura  Basins 

LOVE 

MOJA 

PTDU" 

PVER 

SAFE 

SCLA 

SIVP" 

SNP2 

SNTZ' 

VICE" 

WORK" 

IP 

Mar.  1988 

.  4 

Channel  Islands 

BLUF 

BRSH 

CENT 

MOJA 

NIGU 

PVER 

SBIS 

SMIG 

SOU 

VNDN 

n 

Mar.  1988 

4 

12 

cone  network 

BLAN 

BLHL  ■ 

CENT- 

FIBR  •- 

FTORr- 

LACU 

LOSP . 

MADC 

MOJA 

OVRO 

PVER 

TWIN 

VNDN 

12 

Mar.  1988 

4 

VLBI  sites 

BLAe 

BOUC“ 

JPLI" 

MACA" 

MOJA" 

MONU" 

PEAR" 

PINY' 

SNPA” 

YUMA" 

13 

Mar.  1989 

4 

Channel  Islands 

BLUF 

BRSH 

MOJF 

MOJM 

NIGU 

PVER 

SOU 

TWIN 

VNDN 

14 

Mar.  1989 

4 

11  ■ 

core  network 

BLAN 

BLHL 

CENT 

FIBR 

JPLt 

UCU 

LOSP 

M.ADC 

MOJF 

MOJM 

OVRO 

PVER 

VNDN 

15 

Apr.  1989 

3 

4 

Santa  Maria  Basin 

CAVI 

GRAS 

MILL 

MOJM 

PARG 

VNDN 

16 

Apr.  1989 

3 

Ventura  Ba.sin 

CATO 

HOPP 

MOJF 

MOJM 

MUNS 

PVER 

SAFE 

SNPA 

YAM2 

VFI 

Feb.  1990 

10 

2 

Santa  Maria  Basin 

ALAM 

ALVA 

GAVI 

GRAS 

LOSP 

MOJM 

RUSI 

SYNZ 

VNDN 

17 

Mar.  1990 

4 

7 

core  network 

BLHL 

CENT 

JPLI 

MAJX 

MOJA 

OATtO 

PVER 

(TI4100) 

VNDN 

18 

Mar.  1990 

4 

17 

core  network 

BLAN 

BLHL 

BLUF 

BPA3 

BRSH 

CENT 

CENl" 

(Trimble) 

CEN2 

CEN3 

FIBR 

FORO 

FOR2 

JPLI 

LACU 

LOSP 

MADC 

MOJM 

OVRO 

PIN2 

PVER 

SI02 

VNDN 

VB1 

Jun.  1990 

9 

6 

Ventur*i  Basin 

CATO 

COTR 

HAPY 

HOPP 

LACU 

LOA'E 

MPNS 

MOJM 

MUNS 

PL9A 

PVER 

SAF3 

SAFE 

SCLA 

SNP2 

SNPA 

SOU 

WHIT 

YAM2 

VF2 

Sep.  1990 

14 

3 

Santa  Maria  Basin 

ALAM 

ALVA 

GAVI 

GRAS 

UND 

LOSP 

MADC 

MOJM 

RUSI 

VNDN 

20 

Feb.  1991 

4 

16 

core  network 

BLAN 

BLHL 

BLUF 

BRSH 

CENT 

FIBR 

LACU 

LOSP 

MADC 

MOJM 

NIGU 

OVRO 

POZO 

PVER 

SOU 

TWIN 

AWN 

SBl 

Jun.  1991 

4 

8 

Channel  Islands* 

CENT 

DEVL 

GAVI 

LACU 

MOJM 

OATIO 

PINl 

PVER 

ROCH 

SI01 

SNRI 

SOU 

VAND 

A/NDN 

VB2 

May  1992 

7 

15 

Ventura  Basin 

CATO 

GOLD 

HAP2 

HAPY 

HOPP 

JPLM 

LACU 

LOVE 

MOJM 

MPNS 

MUNS 

PINl 

PVER 

SAFE 

Slot 

SNPA 

SNP2 

SOLI 

WHT3 

YAM2 

a  Site  observed  but  not  included  in  ihe  solutions. 
b  SBl  included  22  Callrans  .silcs  not  lisicd. 


lion  coordinates,  atmospheric  parameters.  Earth  orientation 
parameters  (for  VLBI),  orbital  elements  (for  GPS),  and  pha.se 
ambiguities  (for  GPS).  In  the  second  step,  we  estimate  station 
velocities  in  “multisession”  solutions  which  combine  the  esti¬ 
mates  and  covariance  matrices  from  all  se.ssions.  Station  coordi¬ 
nates  and  ttrbital  parameters  are  also  estimated  in  the  multi.session 
solutions  in  a  consistent  and  non-redundant  manner. 

As  detailed  in  the  appendix,  this  two-step  process  has  two 
advantages.  First,  it  allows  us  lo  handle  the  data  easily.  The 
single-day  solution  condenses  ihc  information  in  the  large  (up  to 
I  Mbylc/station/day)  data  .set  into  a  few  compact  files  which  may 
then  be  used  to  perform  multi.session  solutions  easily  and  quickly. 


Second,  it  affords  a  rigorous  solution  to  the  problem  of  an 
inhomogeneous  tracking  network,  where  the  set  of  stations 
changes  from  day  to  day  and  year  to  year.  Since  this  (“fiducial") 
network  determines  the  frame  to  which  the  estimaled  vectors  are 
referred,  naively  comparing  a  vector  estimated  on  two  days  w  ith 
different  networks  can  lead  to  an  inaccurate  estimate  of  its  rale  of 
change.  The  magnitude  of  the  error  can  reach  1  part  in  10’  for  the 
early  observations  of  our  network  [l.ar.son  d  oL.  lOd]].  As 
di.scu.ssed  in  the  appendix,  our  approach  minimizes  the  effect  of 
Ihe  shifting  fiducial  geometry  by  imposing  the  constraints  on  Ihe 
coordinates  in  a  consistent  manner.  This  is  particularly  important 
for  GPS  tracking  stations  that  have  been  used  only  a  few  limes. 


TABLE  3.  Global  Tracking  Sites 


Expcrinicni  IS  iiicluilcd  daia  also  from  a  Tl  receiver  on  monumenl  867001  al  lire  Pacific  Geosciences  Cenlrc,  Pal  Bay,  British  Columbia, 
The  SB  I  analysis  used  global  orbils  estimated  with  data  from  three  southern  hemisphere  sites  as  well  as  those  shown.  i 
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The  velocities  eslimalcd  I'rom  the  joint  VLBI  and  GPS  solution 
arc  given  in  Table  4  and  arc  discussed  in  the  appendix.  Because 
most  of  our  sites  are  west  of  the  SAP,  we  transfer  the  North 
American  frame  in  which  velocities  are  estimated  to  the  Pacific 
frame  using  the  rcvi.scd  NUVEL-1  model.  These  velocities  are 
shown  in  Figure  2  with  respect  to  the  Pacific  plate.  In  Table  4  and 
the  following  discussions,  we  quote  the  uncertainly  as  one 
standard  deviation.  As  discussed  in  the  appendix,  lhc.se  uncer¬ 
tainties  are  obtained  by  scaling  the  formal  values  by  a  factor  of  2 
R)  reflect  the  scatter  in  the  position  estimates  and  the  presence  of 
sy.stematic  errors.  For  one-dimensional  quantities  quoted  in  the 
text  and  Table  4,  these  .scaled  standard  deviations  .should  be  mul¬ 
tiplied  by  an  additional  factor  of  1 .96  for  testing  hypotheses  at  the 
95%  confidence  level.  In  the  maps  of  the  velocity  fields,  the 
ellipses  denote  the  area  of  95%  confidence  in  two  dimensions, 
after  scaling. 

TTie  precision  and  accuracy  of  the  relative  horizontal  velocities 
estimated  for  the  core  stations  are  adequately  de.scribcd  by  a  95% 


TABLE  4.  Velocities  with  Respect  to  the  Pacific  Plate 


Station 

Ea.st 

Nonh 

Cor. 

Obs. 

Res. 

Unc. 

Obs. 

Res. 

Unc. 

BLAN« 

4.0 

1.4 

1.2 

-2.0 

-3.4 

0.9 

0.091 

BLHL"  - 

2.6 

-0.0 

0.9 

-1.8 

-2.5 

0.6 

0.064 

BLKB7269 

24.3 

2.6 

1.5 

-32.4 

-12.7 

2-2 

-0.092 

BLUF“ 

0.4 

0.7 

1.3 

0.4 

2.9 

0.9 

-0.066 

BRSH 

-0.4 

-0.2 

1.3 

-2.2 

0.8 

0.9 

0.079 

CATO 

0.3 

—0.6 

0.9 

-6.4 

-2.3 

0.8 

-0.091 

CENTF 

1.9 

2.5 

0.8 

-2.4 

-0.2 

0.6 

-0.014 

COTR* 

0.5 

0.2 

1.4 

-2.0 

1.8 

I.l 

0.153 

DEAD7267 

25.1 

1.8 

.3.5 

-33.8 

-1.3.7 

5.1 

0.166 

DEVL 

-0.2 

0.3 

1.2 

-0.7 

1.4 

0.9 

0.198 

FIBR'J 

19.6 

1.0 

1.0 

-23.6 

-6,5 

0.7 

0.120 

FTOR 

-1.4 

-3.3 

1.4 

-0.1 

0.0 

1.6 

-0.047 

GAVI 

-1.0 

-1.4 

1.3 

-3.4 

-2.6 

0.9 

0.114 

GRAS 

-1.7 

-2.9 

1.4 

-0.8 

-0.5 

I.l 

-0.072 

HAPY 

4.3  . 

2.2 

1.3 

-6.0 

-1.0 

1.0 

-0.087 

HOPP 

6.4 

.3.1, 

0.9 

-12.0 

-6.4 

0.7 

-0.149 

JPLl 

2.7 

-0.7 

1.1 

-11.5 

-6.3 

0.9 

-0.047 

LACU'' 

-0.6 

-0.8 

0.7 

-6.3 

-3.7 

0.6 

0.000 

LOSF« 

1.4 

-fl.l 

0.8 

-.3.1 

-.3.2 

0.6 

0.095 

LOVE 

4.1 

-0.6 

I.l 

-1.3.7 

-7.6 

1.0 

-0.064 

MADCA 

4.0 

0.3 

1.1 

-6.5 

-4.0 

0.8 

0.035 

.MOJA“ 

23.9 

-2.0 

0.3 

-26.6 

-7.8 

0.3 

-0.067 

MONP7274 

-0.8 

-2.6 

0.7 

-8.6 

-6.8 

0.9 

0.019 

MUNS 

I.l 

-1.1 

0.9 

-II. 1 

-5.1 

0.8 

-0.1.30 

NIGU 

-1.3 

-1.6 

1.2 

^.3 

-0.8 

0.8 

0.029 

OVRO" 

20.1 

-3.0 

0.5 

-28.0 

-7.1 

0.5 

0.029 

PEAR7254 

14.4 

-0.9 

1.6 

-14.6 

-2.6 

2.4 

0.0.36 

PIN  I 

1.2.3 

2.5 

I.O 

-1.3.2 

-.3.0 

0.9 

-0.063 

POZO* 

4.5 

-0.1 

1.2 

-5.5 

-2.9 

0.7 

0.083 

PRES7252 

8.7 

-9.0 

0.8 

-1.5.9 

1.2 

i.l 

-0.008 

PVER" 

0.7 

0.4 

0.6 

-5.2 

-1.7 

0.5 

-0.0.33 

SAFE 

3.0 

0.1 

0.8 

-9.9 

-4.8 

0.7 

41.107 

SCLA^ 

2.9 

1.6 

2.8 

-7.1 

-2.4 

2.0 

-0.242 

SNP2 

4.9 

2.7 

1.3 

-12.7 

-7.4 

1.0 

41.180 

SNPA 

1.9 

0.1 

0.8 

-8.7 

-3.7 

0.6 

0.024 

SNRi 

-0.3 

0.2 

1.4 

-1.8 

41.5 

1.0 

0.276 

SOLI 

0.0 

-0.2 

0.9 

-8.0 

-4.1 

0.7 

0.0,86 

SOU 

-2.0 

-2.5 

2.4 

-5.9 

-4.1 

1.7 

41.127 

TWIN 

-1.2 

-0.6 

1.2 

-0.8 

1.3 

0.8 

41.065 

VNDN" 

0.3 

-0.4 

0.4 

-1.0 

-1.0 

0.4 

41.044 

\'ekx.'ilics  arc  in  millimeters  per  year. 

Ohs.,  observed;  Res.,  residual:  Unc..  uncertainty  after  scaling  by  2. 
Cor.,  correlation  coelTieient  between  east  and  north  componenLs. 

"  core  site  observed  at  least  5  limes.  _ 

site  observed  only  2  limes,  or  .1  limes  with  Vx’ / >  -■ 


confidence  ellip.se  with  a  .semiminor  axis  of  approximately  2 
mm/yr  oriented  roughly  north-.south,  and  a  semimajor  axis  of  3 
mm/yr  oriented  cast-west.  As  discussed  in  the  Appendix,  this 
level  of  uncertainty  is  supported  by  short-  and  long-term  .scatter 
and  good  agreement  in  velocity  with  independent  VLBI  and  GPS 
.solutions.  Further  calibration  is  provided  by  the  rates  of  shorten¬ 
ing  estimated  by  comparing  our  estimated  line  lengths  with  those 
estimated  from  Irilaleralion.  The  velocity  estimates  are  not 
sensitive  to  errors  in  lies  between  GPS  and  VLBI  monuments, 
because  these  lies  were  not  used.  On  the  other  hand,  they  may  be 
sensitive  to  the  orbital  reference  frame  established  by  the 
available  tracking  stations.  This  sensitivity  is  most  pronounced  in 
the  ea.st  velocity  component  for  stations  observed  only  two  or 
three  times.  For  these  .stations,  noted  in  Table  4,  the  actual  uncer¬ 
tainties  in  ea.stward  vekxrity  may  be  as  large  as  5  mm/yr. 

Removal  of  a  Model  of  the  Effects 
OF  THE  San  Andre.\s  Fault 

In  Figure  2,  the  dominant  feature  in  the  velocity  field  is  the 
simple  .shear  associated  with  the  San  Andreas  fault  system. 
Indeed,  at  this  scale  it  is  difficult  to  discern  any  other  features. 
Since  the  SAF  motion  is  not  the  primary  object  of  our  study,  we 
choose  to  remove  a  reference  model  for  velocities  associated  with 
the  fault  and  to  examine  the  residual  vefociUes. . ' 

A  simple  (and  conventional)  model  attributes  this  shear  to  slip 
at  depth  on  the  SAF  and  associated  faults  below  an  upper. locked- 
portion  of  the  fault.  Our  GPS  network  is  text  sparse,  with  loo  little 
coverage  near  the  faults,  to  reliably  estimate  parameters  in  such  a 
model.  The  best  constraints  on  the  rates  of  slip  and  locking  depths 
come  from  the  more  densely  spaced  irilatcration  surveys  of  the 
U.S.  Geokrgical  Survey  (USGS).  For  the  SAF  northwest  of  its 
intersection  with  the  San  Jacinto  fault,  the  rale  of  slip  estimated 
from  inversion  of  the  geodetic  data  (32-36  mm/yr  [t'ber/iarl- 
F/iiUips  et  ai,  1990:  Lisowski  et  ai.  1991])  is  consistent  with  the 
value  of  34  mm/yr  inferred  from  geology  [Sieh  and  Johns,  1984], 
which  we  adopt  here.  The  locking  depth  in  the  model  is  25  km 
between  Parkllcld  and  San  Gorgonio  Pass  but  shallower  else¬ 
where  (Table  5  and  Figure  3).  For  an  infinitely  long  fault  slipping 
at  veltKiiy  vj  beneath  a  locking  depth  d,  the  predicted  velocity  v 
is  parallel  to  the  fault  and  increases  with  distance  x  from  the  fault 
in  an  arctangent  sigmoid  curve:  v  =  (iN;/a:)  ian-'(.r/(/)  [e.g.. 
Savage  and  Burford,  1973].  Owens  Valley  Radio  Observatory 
(OVRO)  is  our  farthest  regional  site  (-  250  km)  to  the  northeast 
of  the  SAF;  TWIN  on  San  Nicholas  Island  is  our  farthest  site 
(-200  km)  to  the  southwe.si.  For  a  locking  depth  of  25  km.  almost 
32  of  the  total  34  mm/yr  velocity  due  to  deep  slip  on  the  SAF  is 
expected  to  accumulate  between  OVRO  and  TWIN. 

While  this  simple  model  is  useful  for  a  rough  estimate  of  the 
amount  of  deformation  associated  with  the  SAF,  the  actual  fault 
geometry  is  more  complicated.  For  example,  the  locking  depth 
probably  varies  along  strike  [Lisowski  et  a!..  1991],  In  addition, 
the  SAF  lakes  a  left  step  near  its  intcr.scction  with  the  Garlock 
fault  in  the  region  of  the  Big  Bend  and  forms  three  splays  in 
southern  California.  To  include  the  effects  of  this  relatively  well- 
mapped  complexity,  we  use  a  model  that  includes  slip  on  the  San 
Andreas,  San  Jacinto.  Elsinore,  and  Garkx.-k  faults  (Table  5).  We 
choose  not  to  include  the  Eastern  California  Shear  Zone  {Dokka 
and  Travis,  199()a.b:  Savage  el  a!..  190(1]  in  our  model  because 
the  liK'us  of  the  shear  is  not  well  determined. 

We  calculate  the  relative  site  velocities  using  Okada's  [19S5] 
expressions  for  velocities  due  to  slip  on  a  buried  planar  disk'ca- 
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95%  confidence,  scaled  a 


122”W  121 'W  120‘W  119“W  118‘W  117'W  116‘W  115’W  114'W 

Fig.  2.  Observed  velocity  of  stations  relative  to  the  Pacific  plate  estimated  from  the  combined  CPS  and  VLBl  data  .set.  The 
velocities  arc  estimated  in  the  reference  frame  fixed  on  North  America  as  described  in  the  appendix  and  then  transferred  to  the 
Pacific  plate  using  the  rescaled  NUVEL-I  prediction  for  the  relative  motion  of  the  two  plates.  The  ellipses  denote  the  region  of 
d5%  confidence,  after  .scaling  the  formal  uncertainties  by  a  facwr  of  2.  as  described  in  the  text.  For  clarity,  the  ellipses  are  not 
shown  for  the  sites  in  the  Ventura  Basin.  As  listed  in  Table  4.  they  arc  similar  in  size  to  the  ellipse  at  MUNS. 


lion  in  an  clastic  halt-space  with  uniform  clastic  moduli.  For 
simplicity,  wc  assume  a  Poisson  material  and  ignore  the  effects  of 
sphericity  and  variation  of  moduli  with  depth.  Figure  3  shows  the 
velocities  predicted  by  the  dislocation  model  plotted  by  assigning 
zero  velocity  to  the  point  shown  near  the  SW  corner  of  the  map. 
approximately  250  km  oll'shorc  and  375  km  from  the  SAF.  For  a 
point  this  lar  Irom  the  SAF.  the  velocity  predicted  by  the  arc 


tangent  function,  relative  to  points  much  farther  outboard  on  the 
Pacillc  plate,  is  still  almost  1  mm/yr. 

In  addition,  there  are  spatial  variations  in  the  model  velocity 
field  induced  by  the  Big  Bend  and  Oarlock  faults.  The  left  slep  of 
the  SAF  in  the  regiem  ol  the  Big  Bend  results  in  compression 
along  the  direction  ol  plate  motion,  with  extension  perpendicular 
to  this  direction.  The  net  cllect.  in  the  reference  frame  of  the 
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(dislocation  model) 


Fig.  X  VckKilies  predicted  by  our  a  priori  model,  which  includes  Ihe  relative  motion  of  the  Pacific  and  Nonh  America  plates 
as  well  as  deep  slip  on  the  fault  .segments  listed  in  Table  5.  The  faults  are  modeled  as  .slipping  at  the  indicated  rate  below  the 
indicated  depth.  The  reference  frame  is  the  .same  as  in  the  previous  figure,  with  the  offshore  point  at  (.'?.1‘'N.  1 21  'SO’W)  assigned 
zero  velocity.  The  modeled  faults  include  the  San  Andreas  (SAI-SA7),  the  Oarlock  (GAR),  the  San  Jacinto  (SJC),  and  the 
Elsinore  (ELS). 


Pacific  plate,  is  clockwise  rotation  and  fault-normal  velocity  in 
the  region  south  of  the  Big  Bend.  TTte  fault-normal  velocities  of 
the  three  most  northerly  sites  on  the  west  side  of  the  SAF  are  the 
result  of  the  combined  effect  of  the  Big  Bend  and  of  the  shallow 
locking  depth  used  in  their  vicinity. 

In  the  following  discussion,  we  consider  the  residual  velocity 
field:  the  difference  between  the  observ'ed  field  (Figure  2)  and 
that  predicted  by  the  dislcxation  model  (Figure  3).  Figure  4  shows 
the  residual  field,  relative  to  the  offshore  point.  These  residua! 
velocities  arc  thus  referred  to  a  frame  essentially  fixed  on  the 
Pacific  plate.  The  magnitude  of  the  overall  shearing  and  rotation 
is  now  substantially  diminished,  revealing  a  number  of  inlere,sting 
features. 

The  only  group  of  sites  with  residual  velocity  not  re.solvably 
dillci’cnt  from  zero  is  in  the  region  of  the  southern  Channel 
Island.s,  where  San  Nicolas  (TWIN),  San  Clemente  (BLUF).  and 
Santa  Catalina  (BRSH)  islands  all  have  residual  vcUicities  less 
than  2  ±  2  mm/yr.  consistent  with  their  being  on  the  Pacific  plate. 
Their  obserx'cd  velocities  to  the  south  or  SW  with  respect  to  the 
Pacific  plate  (Figure  2)  can  be  explained  almost  entirely  by  our 
model  for  the  effects  of  deep  slip  on  the  SAF  and  as.sociated 
faults  (Figure  4). 


There  is  substantial  convergence  normal  to  the  SAF  in  the 
region  spanning  the  .southern  Channel  Islands,  the  Los  Angeles 
basin,  and  the  Mojave  desert.  For  example,  the  residual 
shortening  of  a  line  oriented  S35°W  bctw'cen  MOJA  and  BLUF 
(on  San  Clemente  Island)  is  10.3  ±  1.0  mm/yr,  wdth  about  half  of 
this  iHtcurring  across  the  Los  Angeles  basin  between  Pasadena 
(JPLl)  and  Palos  Verdes  (PVER),  as  discussed  below.  The 
residual  velocities  of  stations  MUNS.  PEAR,  and  their  neighbors 
point  toward  the  SSW  (relative  to  the  Pacific  plate),  nearly 
perpendicular  to  the  trace  of  the  SAF. 

Similar  fault-normal  velocities  are  evident  in  the  residual 
velocity  field  northeast  of  the  Big  Bend,  where  the  fault-normal 
(S49'\V)  component  of  the  residual  velocity  is  8.2  *  0.6  mm/yr  at 
OVRO  and  4.9  i  0.8  mm/yr  at  FIBR.  Because  of  the  lack  of 
stations  in  this  area,  the  best  measure  of  the  rate  of  compression  is 
3.3  :t  1.0  mni/yr  for  the  fault-normal  component  of  the  residual 
velocity  between  OVRO  and  FIBR. 

Figure  5  shows  the  details  of  the  residual  velitcity  field  plotted 
relative  to  Vandenberg  (VNDN).  To  the  north  of  VNDN.  the 
veliKities  are  dominated  by  shear  parallel  to  the  trace  of  the  SAF 
(S41'E).  The  residual  fault-parallel  motion  of  LOSP.  MADC.  and 
FIBR  increases  systematically  to  the  NE.  with  rates  of  2.2  ±  0.4. 
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Fig.  4.  Residual  velocities,  showing  Ihe  observed  vclocily  field  (Figure  2)  minus  our  model  (Figure  .^).  Residual  vclocilics  are 
plotted  relative  to  the  offshore  point  (dot  at  1 21  °50’W)  on  the  Pacific  plate.  Confidence  ellipses  are  65%.  as  in  Figure  2. 


3.7  *  0.8,  and  6.7  t  0.6  mm/yr,  rc.speclivciy,  relative  to  VNDN. 
To  the  northwest,  sites  BLAN,  BLHL,  and  POZO  have  fault- 
parallel  residual  velocities  of  3.4  ±  1.0,  2.0  e:  0.6,  and  2.6  ±  0.8 
mm/yr  with  respect  to  VNDN.  These  are  the  residual  values  after 
removal  of  the  SAF  shear,  which  is  less  than  1  mm/yr  on  the 
coast  and  4  to  5  mm/yr  for  POZO  and  MADC,  respectively, 
according  to  the  reference  fault  model.  The  fault-normal 
component  of  the  residual  velocity  in  this  area  of  the  Southern 
Coast  Ranges  is  small  for  all  sites,  with  the  largest  value 
occurring  at  FIBR.  which  show’s  2.2  ±  0.8  mm/yr  of  residual 
fault-normal  compression  with  respect  to  VNDN.  There  is  a 
change  in  the  deformation  pattern  near  the  Santa  Maria  Fold  and 
Thrust  Belt.  Convergence  across  the  Santa  Barbara  Channel  and 
the  western  Transverse  Ranges  is  clearly  evident  in  the  LACU- 
CENT  and  MUNS-COTR  vector.s,  respectively.  The  motions  of 
LACU.  MUNS.  and  SOLI  with  respect  to  VNDN  suggest  clock¬ 
wise  rotation  in  the  western  Transverse  Ranges.  These  motions 
w  ill  he  discus.scd  in  their  regional  contexts  below. 

Figure  6  shows  the  residual  velocity  field  plotted  w  ith  respect 
to  Palos  Verdes,  including  interesting  residual  velocities  of  sites 
near  the  Ventura  and  Los  Angeles  basins.  The  residual  conver¬ 
gence  rale  across  the  Ventura  basin  reaches  5.1  ±  1.2  mm/yr  for 
SCLA-SNP2.  w  ith  rapid  clockwise  rotation  apparent.  In  the  Los 
Angeles  basin,  the  residual  convergence  of  site  JPLI  toward 


PVER  is  5.0  ±  1.0  mm/yr  at  SIZ'W  s  14°.  This  compressive 
regime  seems  to  extend  to  Catalina  Island,  w'here  BRSH  has  a 
nonhward  (N13°W  *  30°)  residual  velocity  of  2.6  ±  1.0  mm/yr 
relative  to  PVER.  A  large  group  of  sites  (BLUF,  BRSH,  TWIN, 
CENT,  and  COTR)  have  northerly  velocities  of  3-5  mm/yr 
relative  to  PVER.  This  motion  produces  large  velocity  gradients 
in  the  neighboring  regions.  The  velocity  of  this  group  of  stations 
w’ith  respect  to  the  mainland  (PVER.  NIGU,  and  SOU)  indicates 
shear  in  the  Gulf  of  Santa  Catalina.  Similarly,  rotation  of  the 
bliK'ks  to  the  .south  of  the  Ventura  basin  results  in  left-lateral 
strike  slip  on  the  Santa  Monica  lineament  between  COTR  and 
PVER.  while  the  Santa  Barbara  Channel  ab.sorbs  compression 
between  SOLI  and  CENT. 

Velocity  Gradient  Tensor 

We  have  calculated  the  velocity  gradient  tensor  L  [Malvern, 
1669]  in  each  of  the  triangular  subnetworks  spanning  the 
California  network.  The  triangles  arc  optimally  close  to  equilat¬ 
eral  and  thus  constitute  a  "Delaunay  triangulation”  [Dai'/,s.  1986: 
Gold.  1975;  IPo/.to/i,  1982].  The  symmetric  part  of  L  is  the  strain 
rate  tensor,  whose  eigenvalues  and  eigenvectors  are  written 
according  to  the  convention  of  Prescott  el  al.  [1979]  and  Peitfl  et 
al.  [1990],  The  antisymmetric  part  of  L  gives  a  local  measure  of 
the  average  vorticity.  or  rate  of  rotation.  The  obscrx'ed  values  of 
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Fig.  5.  Residual  velocities  plotted  relative  to  Vandenberg  (VNDN).  Note  that  the  velocity  scale  is  double  that  of  Figure  4. 
Conddence  ellipses  are  <15%,  as  in  Figure  2.  The  tectonic  features  include  the  San  Andreas  fault  fSAF),  the  Santa  ^irbara 
Channel  (SBC),  the  Hosgri  fault  (HSC),  the  Rinconada  fault  (RIN),  the  Santa  Maria  Fold  and  Thrust  Belt  (SMFTB),  and  the 
southern  Coast  Ranges  (SCR). 


the  velocity  gradient  tensor  are  given  in  Table  6.  TTicse  arc  domi¬ 
nated,  in  general,  by  the  shearing  and  rotation  a.ssociatcd  \s'ilh  the 
simple  shear  regime  of  the  SAF. 

The  strain  rates  calculated  from  the  residual  horizontal  velocity 
held  are  shown  in  terms  of  their  principal  axes  in  Figures  7  and  8 
and  li.stcd  with  their  uncertainties  in  Table  6.  Very  high  rates  of 
residual  strain  are  apparent  in  the  Ventura  basin,  while  lower,  but 
nonetheless  interesting,  rates  remain  in  the  southern  Coa.st  Ranges 
and  the  Santa  Barbara  Channel. 

Figures  9  and  10  .show  the  antisymmetric  part  of  the  velocity 
gradient  ten.sor.  the  average  rate  of  rotation  about  a  vertical  axis, 
calculated  for  the  residual  horizontal  velocity  field.  Before 
comparing  these  rotation  rates  to  tho.se  inferred  from  paleo- 
magnetic  declinations  [e.g..  Luyendyk,  IWl],  we  mention  .several 
details.  First,  the  rotation  rales  estimated  from  the  geodetic  data 
are  relercnced  to  the  Pacilic  plate.  Although  such  a  frame  is  not 
identical  to  the  palcomagnetic  frame  defined  by  the  apparent 
magnetic  North  pole,  the  difference  between  them  in  this  case  is 
small  enough  to  neglect.  Second,  we  emphasize  that  Figures  <? 
and  10  show  residual  rotation  rates,  after  removing  transient  shear 
in  our  simple  dishrcaiirrn  fault  model.  Such  a  model  is  nece.s.sary 
il  <ve  are  to  compare  rales  over  two  different  sampling  times.  The 
geodetic  measurements  span  a  period  of  5  years,  considerably 
shorter  than  a  seismic  cycle.  In  contrast,  the  palcomagnetic 


declinations  measure  Unite  rotation  accumulated  during  the  last 
6-16  m.y..  and  thus  represent  the  average  over  many  seismic 
cycles. 

With  these  caveats  in  mind,  we  note  that  the  geodetic  residual 
rotation  rates  do  not  agree  everywhere  with  the  paleomagnetic 
values  summarized  by  Jackson  and  Molnar  [19W].  In  much  of 
the  Transverse  Range.s.  the  average  rate  of  clockwise  rotation 
.since  the  Mkxrcne  is  5-6  7m.y.  [Luyendyk,  l<^<)1].  In  contrast,  the 
geodetic  residual  rotation  rate  varies  from  3  ±  27m.y.  counter¬ 
clockwise  in  the  eastern  end  of  the  Santa  Barbara  Channel 
(triangle  CENT-COTR-SOLI  in  Figure  ‘^)  to  7  i  27m.y.  clock- 
wi.se  in  the  Ventura  basin  (CATO-HAPY-SAFE  in  Figure  10). 
.Around  Santa  Cruz  Island  (CENT),  the  palcomagnetic  declina¬ 
tions  arc  larger  than  elsewhere  in  the  Transverse  Ranges,  indica¬ 
tive  of  rapid  clockwi.se  rotation.  Yet  the  geodetic  residual  rales 
are  less  than  0.5  *  T/m.y.  for  the  five  triangles  around  CENT 
(Figures  9  and  10).  In  the  coastal  part  of  the  Santa  Maria  Basin, 
neither  the  palcomagnetic  declination  nor  the  geodetic  residual 
field  yields  a  rrdalion  rale  resolvably  different  from  zero 
(Figure  9)  [Hornafius  el  a/..  19,86].  Thus  the  residual  geodetic  and 
paleomagnetic  ob.servations  of  rotation  rate  differ  by  an  timount 
which  varies  with  location. 

Why  do  the  two  types  ot  ohser\'alions  disauree?  The  most 
likely  possibility  is  that  rotation  rates  are  not  constant  in 
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geological  time.  This  might  occur  as  faults  rotate  out  of  the 
position  in  which  their  slip  is  favored  [Sco:ii  el  a!..  1991],  In  this 
case,  ihe  area  south  of  the  Ventura  basin  appears  to  be  the 
currently  rotating  element  [Yeats,  1983].  It  may  also  be  that  our 
elastic  dislocation  fault  model  overprcdicts  the  amount  of  rotation 
due  io  elastic  strain  and  leaves  too  little  to  inelastic  strain.  Yet 
another  possibility  is  that  our  geodetic  stations  provide  a  sample 
ol  the  velocity  field  which  is  too  coarse  in  space.  Rather  than 
speculate  on  these  po.ssibiliiies.  we  prefer  at  this  point  simply  to 
present  our  geodetic  estimates,  saving  detailed  modeling  for  a 
future  study. 

Tectonic  Interpretation 

Sonheast  of  the  San  Andreas  Fault 

There  is  substantial  residual  velocity  parallel  to  Ihe  SAF  for 
sites  to  Ihe  northeast  of  this  structure  (Figure  4).  The  sites  farthest 
to  the  northeast.  OVRO  and  MOJA.  have  residual  velocities  with 
fault-parallel  components  of  3.2  ±  0.4  mm/yr  (at  S4rE)  and 
4.5  ±  0.4  mm/yr  (at  ShS^E).  respectively,  relative  to  Ihe  Pacific 
plate.  Near  Ihe  lault.  site  FIBR  has  5.4  ±  0.8  mm/yr  of  residual 
\elocily  parallel  to  the  SAF  (at  S41''E).  Thc.se  rather  large 
residuals  could  be  reduced  by  modifying  Ihe  simple  dishvalion 
model  we  have  u.scd  to  describe  the  effects  of  the  major  faults  in 
ihe  region.  The  possible  changes  include  (1)  decreasing  the 


relative  velocity  of  the  Pacific  plate  relative  to  the  North 
American  plate,  which  provided  the  reference  frame  in  which  the 
velocities  were  determined:  (2)  increasing  the  rate  of  slip  v;  (3) 
reducing  the  locking  depth  d\  or  (4)  including  deep  slip  on  other 
faults.  An  additional  possibility,  di.scussed  in  a  later  section,  is 
that  the  simple  model  of  uniform  slip  at  the  geological  rates  on 
these  laulls  is  inadequate.  In  this  case,  time-dependent  effects  of 
strain  accumulation  during  the  seismic  cycle  must  be  considered. 

The  fir.st  two  factors  are  unlikely  to  be  the  most  important 
ones.  The  plate  rate  has  already  been  reduced  by  2  mm/yr  by  the 
rescaling  of  Ihe  NUVEL-1  motion  model.  Using  the  highest 
proposed  rates  of  slip  for  the  SAF,  36  mm/yr  [e.g..  IJsowski  el  at.. 
1991],  would  account  for  only  a  quarter  of  the  residual  fault- 
parallel  velocities.  In  addition,  it  seems  logically  inconsistent  to 
u.sc  a  slip  rate  other  than  the  geological  average  in  this  simple 
model. 

There  is  .some  evidence  from  our  measurements  that  the 
lockin^depth  on  the  SAF  should  be  shallower.  Station  .MADC. 
the  site  on  the  .southwest  side  of  the  SAF  nearest  to  FIBR,  has  a 
lault-parallel  residual  velocity  of  3.2  ±  0.8  mm/vr.  indicative  of 
some  2  mm/yr  of  iinmodclcd  shear  between  ihe.se  two  sites 
.straddling  the  SAF.  Near  the  fault,  the  residual  shear  could  be 
reduced  by  modeling  the  fault  with  a  shallower  kxking  depth,  but 
this  would  enlarge  Ihe  residuals  of  the  other  stations  southwest  of 
Ihe  fault. 
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TABLE  6.  Velocity  Gradient  Tensor 


Delaunay 
Triangle  . 

Extensive 
Eigenvalue 
(0.1  ppm/yr ) 

Compressive 
Eigenvalue 
(0.1  ppm/yr ) 

Rotation 

Rate 

(deg/m. y. 

) 

Compressive 

Azimuth 

(deg) 

Obs 

Res . 

Unc. 

Obs. 

Res . 

Unc . 

Obs 

.  Res . 

Unc. 

Obs. 

Res . 

Unc. 

BLAN-BLHL-LOSP 

0.5 

0.3 

0.5 

-0.3 

-0.5 

0.4 

-0.0 

-1.3 

1.6 

112 

102 

29 

BLAN-BLHL-POZO 

0.5 

0.1 

0.3 

-1.0 

-0.4 

0.4 

5.4 

1.9 

1.8 

-3 

-21 

13 

BLAN-FTOR-POZO 

0.4 

0.4 

0.1 

-0.5 

-0.4 

0.3 

0.1 

-  -1.7 

0.8 

-  49. 

.  11 

10  ■ 

BLAN-LOSP-VNDN 

0.2 

0.1 

0.3 

-0.7 

-0.8 

0.2 

4.1 

3.0 

0.8 

24 

21 

15 

BLHL-LOSP-POZO 

0.5 

0.1 

0.2 

-0.1 

0.0 

0.1 

3.8 

0.8 

0.7  . 

17 

80 

13 

BLKB-DEAD-MOJA 

2.2 

1.7 

3.4 

-1.0 

-0.7 

1.9 

-7.8 

-6.8 

11.0 

111 

109 

30 

BLKB-DEAD-PINl 

1.6 

0.4 

0.3 

-2.4 

-1.5 

0.7 

11.7 

4.0 

1.8 

9 

24 

4 

BLKB-MOKP-PINl 

1.6 

0.7 

0.3 

-0.8 

-0.2 

0.1 

13.6 

6.5 

1.1 

11 

59 

4 

,BLUF-BRSH-NIGU 

-0.0 

-0.1 

0.4 

-0.5 

-0.4 

0.2 

0.9 

0.7 

i.3 

24 

31 

21 

BLUF-BRSH-TWIN 

0.2 

0.1 

0.1 

-0.5 

-0.4 

0.2 

0.2 

-0.1 

0.9 

9 

8 

15 

BLUF-NIGU-SOLJ 

0.1 

0.3 

0.3 

-0.4 

-0.4 

0.2 

2.0 

2.0 

0.9 

48 

59 

14 

BRSH-NIGU-PVER 

-0.2 

-0-2 

0.3 

-0.8 

-0.7 

0.2 

1.9 

1.4 

1.1 

-6 

-3 

17 

BRSH-PVER-TWIN 

0.1 

0.1 

0.2 

-0.8 

-0.7 

0.3 

1.3 

0.7 

1.1 

-10 

-10 

12 

CATO-COTR-SCLA 

0.1 

-0.0 

0.5 

-1.6 

-1.3 

0.9 

7.7 

6.4 

3.7 

7 

14 

22 

CATO-COTR-TWIN 

0.7 

0.7 

0.4 

-0.4 

-0.5 

0.2 

4.7 

4.3 

1.3 

52 

60 

14 

CATO-HAPY-SAFE 

-0.0 

0.3 

0.5 

-0.6 

-1.0 

0.4 

9.1 

7.4 

1.9 

56 

68 

41 

CATO-HAPY-SCLA 

1.4 

1.2 

1.3 

-0.6 

-0.5 

0.8 

-  '3.0 

1.3 

-  4.1: 

-41' 

128 

-25- 

CATO-JPLl-PVER 

0.4 

0.2 

0.2 

-1.0 

-0.8 

0.2 

3.2 

1.5 

0.8 

9 

23 

5 

CATO-JPLl-SAFE 

0.3 

-0.0 

0.2 

-0.9 

-0.6 

0.4 

4.9 

3.0 

1.4 

-4 

14 

13 

CATO-PVER-TWIN 

0.2 

0.2 

0.1 

-0.5 

-0.3 

0.1 

1.2 

0.7 

0.6 

4 

12 

8 

CENT-COTR-SOLI 

-0.1 

-0.2 

0.3 

-2.4 

-2.2 

0.4 

-2.6 

-3.4 

1.8 

-4 

-7 

7 

CENT-COTR-TWIN 

-0.0 

0.1 

0.2 

-0.4 

-0.6 

0.2 

1.0 

0.1 

0.7 

127 

120 

■19 

CENT-GAVI-LACU 

0.4 

■  0.3 

0.2 

-1.0 

-0.8 

0.2 

1.2 

-0.4 

0.8 

■  28 

24 

7 

CENT-GAVI-SNRI 

0.7 

0.8 

0.4 

-0.3 

-0.3 

0.2 

0.9 

-0.4 

1.1 

3 

-1 

_  14 

CENT-LACU-SOLI 

0.4 

0.2 

0.2 

-1.1 

-1.0 

0.2 

1.8 

0.2 

0.7 

33 

33 

6  ' 

CENT-SNRI-TWIN 

0.7 

0.8 

0.5 

-0.3 

-0.3 

0.1 

2.5 

1.8 

1.0 

-14 

-21 

9 

COTR-SCLA-SOLI 

1.1 

0.8 

1.0 

-2.6 

-2.3 

0.7 

0.6 

-0.9 

3.8  - 

-10 

-11 

8 

DEAD-MOJA-PEAR 

1.0 

0.6 

0.3 

-0.2 

-0.3 

0.2 

4.9 

1.9 

1.2 

25 

45 

10 

DEAD-NIGU-PEAR 

1.0 

0.4 

0.3 

-1.2 

-0.4 

0.2 

9.3 

3.1 

1.1 

-2 

26 

5 

DEAD-NIGU-PINl 

1.2 

0.4 

0.1 

-2.8 

-1.5 

0.7 

6.0 

-  0.2 

1.6 

-7 

3 

4 

FIBR-LOVE-MUNS 

1.4 

0.2 

0.3 

-1.8 

-0.2 

0.  1 

9.6 

2.7 

0.8 

-14 

13 

2 

FIBR-LOVE-PEAR  • 

2.0 

0.8 

0.4 

-1.5 

-0.2 

0.2 

8.0 

-1.1 

1.3 

-22 

118 

4 

FIBR-MADC-MUNS 

1.2 

-0.0 

0.2 

-1.7 

-0.2 

0.1 

12.4 

2.0 

0.5 

-5 

2 

2 

FIBR-MADC-POZO 

1.8 

0.4 

0.2 

-1.2 

-0.1 

0.3 

10.5 

1.7 

1.2 

7 

39 

4 

FIBR-MOJA-OVRO 

0.3 

0.0 

0.0 

-0.2 

-0.1 

0.0 

0.8 

0.3 

0.1 

12 

60 

4 

FIBR-MOJA-PEAR 

0.3 

-0.1 

0.1 

-1.2 

-0.5 

0.2 

3.4 

0.8 

0.5 

-13 

2 

4 

FIBR-OVRO-POZO 

3.0 

0.5 

0.2 

-0.2 

-0.2 

0.0 

4.5 

1.1 

0.5 

36 

45 

2 

FTOR-OVRO-POZO 

0.8 

0.2 

0.1 

-0.5 

-0.2 

0.1 

3.4 

0.7 

0.3 

15 

38 

2 

FTOR-OVRO-PRES 

0.5 

0.4 

0.0 

-1.5 

-0.2 

0.2 

5.1 

0.0 

0.4 

-6 

38 

2 

GAVI-GRAS-MADC 

1.3 

1.1 

0.7 

-0.3 

-0.2 

0.2 

6.5 

3.8 

1.6 

17 

22 

10 

GAVI-GRAS-VNDN 

0.8 

0.8 

0.5 

-0.7 

-0.7 

0.3 

0.4 

-1.0 

2.0 

65 

63 

12 

GAVI-LACU-MADC 

0.1 

0.2 

0.3 

-0.4 

-0.2 

0.1 

4.7 

2.0 

0.8 

-7 

-0 

17 

GAVI-SNRI-VNDN 

0.1 

0.0 

0.1 

-0.7 

-0.7 

0.4 

2.0 

1.1 

1.1 

45 

40 

14 

GRAS-LOSP-MADC 

0.3 

0.1 

0.2 

-1.7 

-1.6 

0.5 

5.4 

3.0 

1.8 

-25 

-31 

9 

GRAS-LOSP-VHDN 

-0.3 

-0.4 

0.2 

-1.6 

-1.7 

0.7 

-0.4 

-1.6 

1.7 

116 

116 

15 

HAPY-HOPP-LOVE 

-1.3 

-2.0 

0.7 

-4.7 

-4.1 

0.9 

5.6 

2.2 

3.6 

-9 

-4 

11 

HAPY-HOPP-SNPA 

2.4 

1.9 

0.8 

-4.7 

-4.1 

0.9 

3.0 

-0.1 

3.6 

-11 

-9 

5 

HAPY-LOVE-SAFE 

-0.2 

-0.1 

0.6 

-3.0 

-3.0 

0.6 

7.7 

.  4.1 

2.6 

18 

33 

9 

HAPY-SCLA-SNPA 

1.3 

1.0 

1.1 

-3.1 

-2.7 

2.6 

-8.9 

-10.7 

11.4 

5 

9 

25 

HOPP-LOVE-MUNS 

1.8 

2.1 

0.7 

-3.7 

-4.0 

1.4 

-10.6 

-13.7 

5.1 

42 

49 

8 

HOPP-MUNS-SNP2 

1.3 

1.6 

0.7 

1.0 

0.5 

1.1 

-0.8 

-3.1 

4.6 

3 

69 

165 

HOPP-SNP2-SNPA 

1.4 

0.7 

0.8 

-8.2 

-7.6 

2.3 

7.2 

4.8 

7.0 

-25 

-27 

8 

JPLl-NIGU-PEAR 

2.7 

0.8 

0.6 

-1.2 

-0.8 

0.2 

5.7 

-3.6 

2.0 

-14 

132 

5 

JPLl-NIGU-PVER 

-0.1 

-0.3 

0.2 

-1.1 

-0.8 

0.2 

2.7 

0.9 

0.7 

-5 

13 

7 

JPLl-PEAR-SAFE 

1.6 

1.1 

0.4 

-1.1 

-0.1 

0.5 

10.3 

0.8 

1.7 

-28 

92 

8 

LACU-MADC-MUNS 

0.2 

-0.1 

0.2 

-0.5 

-0.2 

0.2 

5.7 

1.6 

0.8 

9 

28 

11 

LACU-MUNS-SOLI 

0.5 

0.2 

0.2 

-0.8 

-0.4 

0.2 

3.9 

0.5 

1.1 

15 

29 

8 

LOSP-MADC-POZO 

0.4 

0.3 

0.3 

-0.1 

0.0 

0.2 

3.5 

1.1 

1.1 

9 

50 

18 

LOVE-PEAR-SAFE 

1.6 

0.0 

0.3 

-2.1 

-1.3 

0.6 

4.0 

-2.7 

2.2 

-8 

-9 

6 

MONP-PINl-SOLJ 

0.6 

0.5 

0.2 

-1.1 

-0.1- 

0.2 

.6.5 

3.4 

0.8. 

-33 

105  - 

-  6-  -  - 

MUNS-SNP2-SOLI 

1.7 

1.4 

0.5 

-0.8 

-0.4 

0.2 

4.3 

2.3 

1.4 

12 

23 

6 

NIGU-PINl-SOLJ 

1.2 

0.4 

0.1 

-0.2 

0.3 

0.2 

4.9 

2.0 

0.9 

-1 

-35 

7 

SCLA-SNP2-SNPA 

-2.2 

-2.2 

3.1 

-12.4  - 

■11.9 

5.8 

-8.0 

-9.6 

15.8 

126 

124 

19 

SCLA-SNP2-SOLI 

1.2 

0.8 

0.8 

-4 . 8 

-4.3 

2.2 

1.9 

-0.6 

7.4 

-11 

-10 

14 

Obs . ,  observed;  Res.,  residual;  Unc. ,  uncertainty  after  scaling  by  2.0. 
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Fig.  7.  Principal  axes  of  ihe  horizontal  strain  rale  tensors  in  the  area  around  Vandenbcrg,  calculated  for  the  residual  velocity 
field  shown  in  Figure  5.  In  each  Delaunay  triangle,  the  inward  pointing  arrows  represent  compression:  outward  pointing  arrows 
represent  extension.  If  neither  principal  strain  rate  is  larger  in  magnitude  than  its  uncertainty  or  if  the  orientation  of  the  axes  has 
an  uncertainty  greater  than  30°,  then  the  axes  are  not  plotted.  Values  and  uncertainties  for  all  triangles  are  given  in  Table  6. 


There  is  independent  geodetic  evidence  that  we  have  neglected 
important  structures  near  the  northeast  edge  of  our  network.  Sites 
OVRO  and  MOJA  arc  near  the  Eastern  California  Shear  Zone, 
which  Irom  terrestrial  geodetic  evidence  accumulates  approxi¬ 
mately  8  mm/yr  of  fault-parallel  velocity  across  it  [Savoge  ei  ai, 
1990],  A  rate  of  this  magnitude  is  sufficient  to  explain  the  fault- 
parallel  components  of  the  residual  velocities  we  e.stimate  at  these 
two  stations. 

In  addition  to  the  fault-parallel  component  of  residual  velocity, 
ihe  sites  northeast  of  the  SAP  have  a  significant  component  of 
velocity  to  the  SW.  perpendicular  to  the  trace  of  the  SAP.  This 
lault-normal  velocity  cannot  be  explained  by  adjusting  the 
locking  depths  of  our  model  faults.  Por  example,  there  is  fault- 
ncrrmal  compression  of  5.7  ±  0.5  mm/yr  between  OVRO  and 
BLAN.  The  existence  of  this  comprc.ssion  is  consistent  with  the 
formation  of  the  anticlinal  structures  that  trap  petroleum 
[Callaway.  1971]  and  the  thrust  mechanism  of  the  19, S3  Coalinga 
earthquake  [.Stein  and  King.  1984].  To  our  know  ledge,  this  is  the 
first  delcrmination  of  the  total  rate  of  shortening  across  these 
Structures. 

Southern  Coast  Ranges 

The  residual  velocities  in  the  .southern  Coast  Ranges  (Pigurc  5) 
are  predominantly  parallel  to  the  SAP.  One  possible  explanation 


for  the.se  residuals  is  that  the  locking  depth  for  the  SAP  is  too 
shallow  in  the  simple  reference  model.  As  discussed  above,  it 
seems  more  likely  that  the  modeled  depth  is  too  deep,  in  which 
case  the  residual  fault-parallel  velocities  are  more  likely  the  result 
of  .strain  accumulation  on  .some  other  structure.  The  San  Simeon 
strand  of  the  Hosgri-San  Gregorio  fault  sy.stem  [Racific  Gas  and 
electric.  1988],  a  plausible  candidate  for  such  strain  accumula¬ 
tion.  lies  between  our  stations  at  Point  Piedras  Blancas  (BLAN) 
and  Black  Hill  (BLHL).  near  Morro  Bay.  The.sc  two  stations 
exhibit  a  relative  residual  velocity  of  1.7  =:  (1.7  mm/yr,  a  result 
only  marginally  different  from  no  deformation.  Despite  its  un¬ 
expected  left  lateral  direction  (S57E  ±  42°)  direction,  this  vector 
is  not  sufficiently  precise  to  exclude  activity  on  the  fault.  Even  if 
the  fault  were  active,  these  stations  wiiuld  caplure  very  little 
motion  between  them  if  the  fault  were  locked  to  great  depth, 
because  both  stations  are  quite  close  to  the  fault  (less  than  5  and 
15  km  lor  BLAN  and  BLHL.  respectively).  For  example,  if  the 
lault  were  locked  to  25  km  depth,  and  slipping  at  2  mm/vr  bek'w 
that  depth,  the  arc  tangent  dislocation  calculation  would  predict 
le.ss  than  0.5  mm/yr  of  motion  between  BL.VN  and  BLHL. 
.Although  such  strain  accumulation  would  be  right-lateral,  it 
would  be  difficult  to  detect  with  our  geodetic  measurements. 
Nonetheless,  there  is  a  slight  dillerence  between  our  short-term 
geodetic  left-lateral  rate  and  the  long-term,  right-lateral  rates  of 
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Fig.  8.  Principal  axes  of  the  horizontal  strain  rale  tensors  in  the  area  around  Palos  Verdes,  calculated  for  the  residual  velocity 
field  shown  in  Figure  4.  Plotting  conventions  as  in  Figure  7.  Note  that  the  scale  for  the  strain  rale  is  different  from  that  of  the 
previous  figure  by  a  factor  of  5. 


1.25  ♦  0.5  mm/yr  interred  from  trenching  and  the  1-3  mm/yr 
inferred  from  offset  terraces  and.  drainages  {Pacific  Cas  and 
Electric,  1988].  Given  such  a  small  ratio  of  signal  to  noise,  how¬ 
ever,  the  two  types  of  estimates  are  compatible,  unless  the  fault  is 
slipping  near  the  surface. 

The  significant  fault-parallel  residual  velocities  of  sites  POZO 
and  MADC  indicate  a  shortcoming  in  the  reference  model.  As 
discussed  above,  these  residuals  would  be  increa.sed  by  a  model 
with  a  locking  depth  less  than  25  km.  To  decrease  their 
magnitude  would  require  an  implausibly  deep  locking  depth  on 
the  SAF.  An  alternative  explanation  is  strain  accumulation  on  a 
previously  unrecognized  onshore  structure,  or  on  an  offshore 
structure,  such  as  the  Hosgri  fault. 

In  the  San  Luis  trilateration  network,  fault-normal  comprc.ssion 
has  been  inferred  by  Harris  and  Sei>all  [1987]  using  data 
collected  between  1959  and  1984.  Although  such  motion  is 
kinematically  compatible  with  the  notion  that  the  San  Andreas 
fault  is  “weak”  [Zoback  et  nl.,  1987:  Mount  and  Suppe,  1987],  it 
has  not  been  observed  by  Lisowski  el  al.  [1991]  for  other  USCS 
trilateration  networks,  or  by  Shen  and  Jackson  [1993]  from  the 
combination  of  GPS  and  early  triangulation  data.  Doni’  [1993] 
finds  fault-normal  extension  for  the  same  network  from  a  data  set 
that  includes  USGS  trilateration  data  from  1980  to  1990.  but  the 
extension  rate  is  not  statistically  significant.  Our  CPS  estimates 
provide  independent  evidence  on  this  apparent  contradiction. 


Although  BLHL  is  the  only  station  common  to  both  our  GPS 
network  and  the  trilateration  network  of  Harris  and  Segall  [1987] 
and  Dong  [1993],  an  approximate  comparison  is  possible. 
Between  BLHL  and  Chiches  (near  POZO),  Harris  and  Segall 
[1987]  find  a  fault-normal  velocity  of  3  ±  3  mm/yr.  The  fault- 
normal  (N49°E)  component  of  our  residual  velocity  between 
BLHL  and  POZO  is  0.5  ±  1.2  mm/yr.  which  is  not  significantly 
different  from  zero  or  from  the  analysis  of  the  1980-1990  USGS 
trilateration  data  [Dong,  1993].  A  less  direct  comparison  involves 
lines  spanning  the  SAF,  such  as  BLHL  to  FIBR,  almost  50  km  off 
the  fault  on  the  east  side.  For  this  line,  w'e  find  3.6  ±  1 .0  mm/yr  of 
residua!  fault-normal  comprc.ssion.  which  is  less  than  the  6.1  * 
1.7  mm/yr  rate  of  residual  fault-normal  shortening  between 
BLHL  and  a  point  10  km  east  of  the  SAF  estimated  by  Harris  and 
Segall  [1987]  but  greater  than  the  value  inferred  by  Dong  [1993]. 
A  balanced  cross  section  including  this  area,  but  extending 
offshore  across  the  trace  of  the  Hosgri  fault  system,  yields  6-13 
mm/yr  of  shortening  [Nainson  and  Davis.  1990],  apparently 
indicating  more  deformation  than  measured  in  our  network,  but 
di.siributcd  over  a  larger  region. 

Santa  Maria  Fold  and  Thrust  Belt 

Of  the  five  stations  in  the  Santa  Maria  Fold  and  Thrust  Belt, 
LOSP  and  MADC  exhibit  residual  motions  significantly  different 
from  zero  at  95'%  confidence  (Figure  5).  In  particular,  the  residual 
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Fig.  9.  Avcrap  rolalion  rales  infcrTcd  from  Ihc  residual  velocity  field  for  ihe  area  around  Vandenberg.  The  rales  are  shown  as 
gray  fans,  in  which  each  fold  denotes  20  nrad/yr  (l.lVm.y.).  The  while  wedge  represents  plus  or  minus  one  standard  deviation, 
after  the  formal  uncertainties  have  been  scaled  by  a  factor  of  two.  If  Ihe  uncertainly  on  Ihe  rotation  rale  exceeds  88  nrad/yr 
(5°/m.y.),  Ihe  rolalion  rale  is  not  shown.  Values  and  uncertainties  for  all  triangles  are  given  in  Table  6. 


velocity  of  MADC  with  respect  to  VNDN  implies  2.3  ±  0.8 
mm/yr  of  shortening  and  3.1  ±0.8  mm/yr  of  right-lateral  strike 
slip  in  the  frame  defined  by  the  local  (N60'W)  strike  of  the  folds. 
The  shortening  component  can  be  interpreted  in  terms  of  active 
shortening  on  the  WNW-ESE  trending  folds  in  the  area.  The  rate 
of  shortening  is  compatible  with  the  geological  rate  of  2-5  mm/yr 
estimated  from  a  balanced  cross  section  [Namson  and  Davis, 
1990]  but  smaller  than  the  6  ±  1  mm/yr  estimated  by  a 
comparison  of  GPS  and  historical  survey  data  [Feigl  et  al.,  1990]. 
The  latter  study,  however,  assumed  uniform  strain  in  the  SMFTB, 
no  net  rotation,  and  no  strain  accumulating  from  the  locked  SAP. 
■All  three  of  these  hypotheses  are  suspect  in  light  of  the  spatially 
varying  .strain  rates  (Figure  7),  the  suggestion  of  rotation  (Figure 
9).  and  the  4  mm/yr  of  relative  motion  between  VNDN  and 
MADC  in  our  fault  model  (Figure  3). 

The  amount  of  right-lateral  .strike  slip  shear,  however,  is  larger 
than  expected,  given  the  lack  of  laulls  active  in  the  Quaternary 
[California  Division  Mines  and  Geology,  1992].  One  explanation 
would  be  unmodelcd  .strain  accumulation  on  the  SAF;  in  light  of 
the  evidence  lor  a  shallower  locking  depth  for  the  SAF.  this 
explanation  seems  unlikely.  Allernalivciy.  strain  accumulation  on 
the  ollshorc  Ho.sgri  fault  [Hull,  1978;  1981]  may  be  indicated, 
although  modeling  this  lault  seems  to  require  an  implausibly  deep 
liKking  depth  (>  25  km)  and/or  a  rapid  slip  rate  (>  10  mm/yr)  to 


explain  the  fault-parallel  residual  velocities  at  POZO  and  MADC. 
A  geometrically  simple  .solution  would  be  right-lateral  shear 
accumulating  on  a  NW  striking  fault  between  VNDN  and  the 
Southern  Coast  Ranges  (LOSP,  BLHL,  MADC,  POZO,  and 
BLAN),  but  there  is  little  geological  evidence  for  Quaternary 
activity  on  such  an  onshore  structure  [Sylvester  and  Darrow, 
1979;  Clark  el  al.,  1984;  Pacific  Gas  and  Electric,  1988],  An 
additional  po.ssibilily  is  that  our  model  of  elastic  strain  accumula¬ 
tion  on  the  SAF  is  too  simple  because  it  neglects  viscoelastic 
clfects  through  the  earthquake  cycle.  But  since  the  uncertainties 
are  a  substantial  fraction  ol  the  residual  velocities  in  this  region, 
we  po.stpone  interpretation  of  our  results  in  this  area  until  w-e  have 
more  data.  A  program  to  continue  and  density  the  measurements 
is  currently  underway. 

Santa  Barbara  Channel 

Acro.ss  Ihc  Santa  Barbara  Channel,  we  find  up  to  4.9  ±  1.0 
mm/yr  of  residual  shortening  (CENT-LACU  in  Figure  5).  with 
the  principal  compressive  strain  oriented  NE-SW  (Figure  7). 
Similarly,  a  combination  ol  trilateraiic»n  and  GPS  measurements 
indicates  northeast  trending  convergence  of  6.4  ±  0.9  mm/yr 
between  1970  and  1988  [l.arsen  ci  at.,  1993],  Part  of  this  defor¬ 
mation  may  be  attributable  to  lault-bend  folding  on  the  offshore 
Oak  Ridge  and  Blue  Bottle  trends  described  by  Shaw  and  Suppe 
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Fig.  10.  Average  rotation  rates  inferred  Irom  the  residua!  velocity  field  for  the  area  around  Palos  Verdes.  Plotting  conventions 
as  in  Figure  9. 


[1991].  In  the  central  third  of  the  channel,  in  the  area  bounded  by 
Santa  Cruz  Island  (CENT),  La  Cumbre  Peak  (LACU),  and 
Gaviota  Peak  (GAVI),  we  find  more  left-lateral  shear  across  a'^ 
vertical  plane  striking  ENE  than  we  do  to  the  east,  where  the 
strain  is  essentially  uniaxial.  Larsen  el  al.  [1993]  suggest  a 
change  in  the  deformation  pattern  between  these  two  regions, 
with  uniaxial  strain  to  the  cast  and  left-lateral  strike  slip  to  the 
west.  Our  observations  are  consistent  with  a  change  betw'een  the 
eastern  and  central  thirds  of  the  channel,  although  we  find  some 
dilatational  strain,  in  addition  to  strike  slip,  in  the  central  section. 

Ventura  Basin 

The  rapid  convergence  observed  across  the  Santa  Barbara 
Channel  continues  to  the  east  into  ihe  Ventura  basin,  where  the 
close  spacing  of  the  CPS  stations  allow's  us  to  locate  large 
gradients  in  velocity.  The  convergence  rate  across  the  basin  is  5-8 
mm/yr  [Donnellan  et  al.,  1993,  this  issue].  The  convergence  takes 
place  within  the  narrow  basin,  resulting  in  a  maximum  shear  strain 
rate  of  0.6  ±  0.1  /rrad/yr,  with  the  azimuth  of  maximum  compres¬ 
sion  oriented  approximately  N-S,  averaged  over  the  east  and  central 
parts  of  the  basin  [Donnellan  et  al. ,  this  issue].  This  strain  rate  is 

almost  twice  that  associated  with  the  San  Andreas  fault  in  this 
region.  The  areal  dilatation  rate.  -5  ±  1  x  Ifi-^  ym'.  is  larger  in 
magnitude  than  any  other  observed  in  California.  The  deforma¬ 
tion  pattern  in  this  area  is  consistent  with  slip  on  thrust  faulLs, 
interred  Irom  scismological  and  geological  observations,  loading 


the  upper  5-1,0  km  of  the  faults  that  bound  the  basin  [Donnellan 
et  al..  1993].  The  residual  velocities  of  5-8  mm/yr  acro.ss  the 
basin  arc  less  than  half  those  proposed  previously  [Yeats,  1983]. 
One  way  of  reconciling  the  geologically  inferred  rates  of  slip  on 
the  faults  bounding  the  basin  with  the  geodetically  observed 
convergence  rates  is  if  the  faults  maintain  their  (relatively  sleep) 
near-surface  dips  to  great  depths,  rather  than  llaltcning  out  into 
sub-horizontal  dccollements  [Yeats  and  Huftile,  1992]. 
Alternatively,  the  fault  geometries  used  previously  could  be 
correct,  with  an  error  in  the  inferred  lime  of  deformation 
(P.  Molnar,  written  communicalion,  1992). 

Los  Angeles  Basin 

The  residual  velocity  of  Palos  Verdes  (PVER)  with  respect  to 
the  Jet  Propulsion  Laboratory  (JPLl),  north  of  Pasadena,  implies 
(1.2  ±  1.2  mm/yrof  right-lateral  shear  and  5.0  ±  1.2  mm/yr  of 
shortening  in  a  coordinate  system  defined  by  the  local  strike  of 
the  SAP  (N75‘’W).  A  balanced  geological  cross  .section  provides  a 
minimum  rale  of  3. 8-6. 8  mm/yr  of  shortening  over  the  last  2-4 
mry:  along  a  line  including,  but  longer  than,  the  .segment  between 
PVER  and  JPLl  [T.  /..  L)a\is  el  al..  1989].  Assuming  uniform 
strain,  the  .section  between  PVER  and  JPLl  is  shortening  al  2.5- 
4.5  mm/yr.  consistent  with  our  measurement.  There  are  two 
interpretations  of  our  measurement  being  near  the  upper  bound  of 
the  geological  determination.  This  may  mean  that  the  geologic 
rates  are  near  the  upper  limit  of  the  published  rates  and  that  the 
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struclurcs  accumulating  strain  arc  nearly  completely  spanned  by 
this  baseline.  Alternatively,  strain  could  be  accumulating  on  the 
frontal  fault  of  the  San  Gabriel  Mountains,  as  proposed  by  Bird 
and  Rosenslock  [1984],  although  at  a  rate  lower  than  ihey 
proposed.  Identification  of  the  features  responsible  for  accommo¬ 
dating  the  shortening  should  be  possible  with  geodetic  measure¬ 
ments  currently  underway  [Hudnut  et  al.,  19*51;  Slien  el  al., 
1991). 

Possible  Time  Dependeni  Complications  ~- 

Thc  period  of  our  observations  between  1986  and  1991  may  be 
considered  to  be  interseismic.  During  this  inters'al,  the  only 
important  (Af  >  6)  earthquake  to  occur  in  the  vicinity  of  our  net¬ 
work  was  the  (A/m)  6.9  Loma  Prieta  event  [U.SGS  Staff,  1990; 
Kanainori  and  Satake,  1990].  Although  it  produced  co.scismic  slip 
of  approximately  2  m  [l.isowski  et  al..  1990],  it  probably  did  not 
displace  our  benchmarks,  all  but  one  of  which  were  over  100  km 
away.  The  only  station  in  our  network  likely  to  have  been 
displaced  by  the  co.seismic  motion  is  at  Fort  Ord,  which  moved 
49  ±  4  mm  to  the  .southea.st  [Clark  et  al.,  1990].  Although  such  a 
displacement  would  ordinarily  be  re.solvable  by  GPS.  it  does  not 
appear  in  our  result  because  our  GPS  observations  at  Fort  Ord 
involve  different  bench  marks  before  (FTOR)  and  after  (FT02) 
the  earthquake.  In  the  multisession  solutions,  the  positions  of  the 
two  monuments  are  not  tied  together.  The  velocities  of  all  the 
marks  at  Fort  Ord  (both  GPS  and  VLBl)  are  assumed  to  be  equal 
in  the  preseismic  and  posiseismic  periods. 

The  Whittier  Narrows  earthquake  =  6-0)  [Hauksson  et  al., 
1988]  event  occurred  on  October  1,  1987,  during  our  survey  in 
the  Santa  Maria  area.  Although  we  made  observations  following 
the  earthquake,  we  w'ere  unable  to  find  stations  which  both  had 
been  previously  occupied  and  w’ere  sufficiently  close  to  the 
epicenter.  The  stations  at  jPL  and  Palos  Verdes,  the  core  sites 
clo.scst  to  the  epicenter,  are  too  far  away  to  have  measurable  co¬ 
.seismic  displacement. 

A  more  .serious  time-dependent  complication  in  our  analysis 
may  result  from  the  simplicity  of  the  model  which  we  use  to 
remove  approximately  the  effects  of  deep  slip  on  the  SAF.  We 
know  that  the  Earth  is  not  an  elastic  half-space,  and  we  suspect 
that  faults  do  not  slip  with  uniform  velocity  below  some  locking 
depth.  During  the  interseismic  period,  viscoelastic  effects  arc 
likely  to  be  important.  We  have  chosen  to  ignore  the.se  effects  in 
our  simple  reference  model  because  the  parameters  to  be  used  in  a 
more  realistic  model  arc  not  well  constrained.  Some  of  the 
residual  velocity  Held  may  well  result  from  ignoring  these  effects. 

A  time-dependent  velocity  field  during  the  interseismic 
interval  would  result  either  from  temporally  varying  slip  on  the 
fault  plane,  with  more  rapid  slip  immediately  following  an  earth¬ 
quake  [c.g..  Thatcher.  1983],  or  by  viscoelastic  relaxation  in  an 
intracru.stal  or  subcrustal  asthenospheric  layer  [c.g..  I.i  and  Rice. 
1987].  Savage  [1990]  has  shown  that  the  surface  deformation  for 
an  elastic  lithosphere  overlying  a  viscoelastic  half-space  can  be 
simulated  by  an  appropriate  variation  of  slip  rate  on  a  fault 
embedded  in  an  elastic  half-space.  For  cither  model,  early  in  the 
earthquake  cycle,  the  fault-parallel  velocities  arc  greater  than  the 
geologic  average,  w  hile  late  in  the  earthquake  cycle,  the  velocities 
fall  behind  the  geologic  average.  For  some  choices  of  paramelcrs. 
early  in  the  earthquake  cycle,  sites  near  the  fault  have  higher 
N'elocilies  than  those  farther  away;  the  plot  of  velocity  versus 
distance  from  the  fault  is  no  longer  monotonic,  as  it  is  for  the 
simple  arc  tangent  model.  In  these  models,  the  deep  slip  velocity 


that  would  be  inferred  from  a  simple  inversion  of  velocities  near 
the  fault  would  not  match  the  geologically  inferred  rale,  except  by 
coincidence.  In  addition,  the  slip  rale  inferred  from  observed 
velocities  near  the  fault  would  differ  from  that  inferred  from 
velocities-obscrx'edTarlher  away. 

Some  of  the  residual  velocities  that  w'c  llnd  after  subtracting 
the  predictions  of  a  simple  model  of  faults  slipping  at  their 
geological  rales  from  our  measured- velocities  may  well  be 
evidence  for  time-dependent  velocities  such  as  lho.se  addrc.s.sed  by 
5’oi'nge-[199()].  Before  using  our  ob.served  velocities  to  infer 
reliably  the  existence  of  slip  on  other  structures,  it  will  be 
necessary  to  address  this  problem. 

Resolution  of  the  San  Andreas  Discrepancy 

The  kernel  of  the  San  Andreas  di.scrcpancy  argument  is  that  not 
all  of  the  deformation  at  the  plate  boundary  is  accommodated  on 
the  SAF  [Minster  and  .Iordan.  1984.  1987].  The  part  of  the 
Pacific-North  America  relative  motion  that  is  accommodated 
neither  on  the  SAF  nor  in  the  Basin  and  Range  has  been  named 
the  “modined  discrepancy’’  [.4rgu5  a/i(f  G'or</o/i,  1991],  It  is  the 
deformation  we  expect  to  observe  in  California.  Since  it  is  ba.sed 
on  a  plate  motion  model,  this  vector  is  a  geologic  estimate,  and 
thus  averages  over  many  seismic  cycles.  Our  geodetically 
obseryed  vectors,  on  the  other  hand,  span  much  less  than  the 
typical  recurrence  time  for  a  great  earthquake.  To  account  for  the 
elastic  strain  accumulation  during  this  inter-seismic  lime  period, 
we  consider  the  velocity  predicted  by  the  reference  model  of 
cla.stic  strain  accumulation  due  to  deep  slip  on  known  faults.  The 
velocity  estimates  from  these  three  sources,  geodetic  obserx'ation, 
elastic  modeling,  and  geological  estimation,  can  be  compared  to 
describe  the  partitioning  of  the  Pacific-North  America  motion. 

The  first  question  that  we  address  is  how  much  of  the  motion 
between  the  Pacific  and  North  American  plates  occurs  on  shore. 
In  Figure  I  la,  we  compare  the  geodetic  estimate  and  two 
geological  values  from  plate  motion  models  for  the  motion  of 
Vandenberg  (VNDN)  with  respect  to  North  America.  The 
geodetically  ob.served  velocity  of  VNDN  is  within  the  95% 
confidence  ellipse  for  the  values  predicted  by  both  the  original 
and  the  rescaled  (slower)  NUVEL-1  plate  motion  models.  From 
this  we  infer  that  the  bench  mark  at  VNDN  is  riding  on  the 
Pacific  plate  to  within  3  mm/yr.  Although  the  geodetic 
observation  is  closer  to  the  rescaled  NUVEL-I  velocity  than  to 
the  original  NUVEL-1  value,  we  cannot  reject  the  original  value 
with  any  reasonable  degree  of  confidence.  Similarly,  the  NUVEL 
uncertainties  are  large  enough  to  allow  some  1-2  mm/yr  of 
motion  offshore  of  Vandenberg.  The  relative  motion  between 
VNDN  and  the  .southern  Channel  Islands  (Figure  4)  suggests  that 
this  is  the  case,  although  the  magnitudes  of  their  relative 
velocities  are  less  than  the  uncertainties  in  the  plate  rates.  Since 
we  have  no  .stations  in  the  Pacific  Ocean  offshore  of  Vandenberg, 
the  only  other  measurement  bearing  on  this  problem  is  the 
3.0  ±  0.3  mm/yr  of  change  in  distance  between  VNDN  and  Kokec 
Park,  Hawaii,  a  .station  presumably  on  the  Pacific  plate.  Although 
this  result  is  suggestive  of  deformation  offshore  of  VNDN.  it  may 
also  be  due  to  volcanic  deformation  at  Kokee,  hK'ai  instabilities  at 
either  site,  or  tropospheric  biases  in  Hawaii.  The  value  of  3 
mm/yr.  then,  is  best  interpreted  as  the  uncertainly  with  which  we 
can  state  lhat  VNDN  rides  on  the  Pacific  plate  and  that  the 
"missing'’  motion  occurs  befween  VNDN  and  North  America. 

The  question  then  becomes,  how  is  this  motion  distributed?  To 
address  this  question,  we  compare  the  (rbserved  and  modeled 
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VNDN  w.r.t.  N.  America 
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East  (mm/yr) 

Fig.  II.  Panition  of  ihe  relalive  velocily  belwcen  the  Nonh  America  and  Pacific  plales.  (a)  Three  cslim_ales  for  the  tola! 
velocity  of  ihe  Pacific  plate  with  respect  to  North  America:  original  NUVEL-I  (large  while  ellipse),  rescaled  NUVEL-I  (gray 
ellipse),  and  our  geodetic  estimate  for  Vandenberg  (VNDN).  (6)  Three  vector  sums  for  ihe  partition  of  this  total  velocily,  plotted 
with  respect  to  Vandenberg:  geodetic  observations  (light  gray  ellipses,  with  station  names  labeled  in  roman  text  at  right), 
geologic  estimates  (open  arrows  with  open  confidence  regions  and  features  labeled  in  italic  text  at  left),  and  predictions  of  the 
reference  elastic  model  (small  solid  dots  connected  by  dashed  lines).  The  thick  solid  lines  connecting  the  geodetic  observations 
with  the  model  predictions  arc  thus  residual  velocities.  Finally,  the  prediction  of  the  rc.scaled  NUVEL-1  model  for  North 
America  with  respect  to  VNDN  is  shown  at  lower  right  as  a  large,  dark  gray  ellipse  [DcMcis  cl  at.,  I9P0;  C.  DeMets,  personal 
communication,  I99.T|.  Geologic  estimates  include  the  conventional  value  for  slip  on  Ihe  San  Andreas  fault  (SAF)  from  Minster 
and  Jordan  11684];  extension  in  the  Basin  and  Range  (B&R)  from  the  geologic  estimate  (model  A)  of  Minster  and  Jordan 
[1687|:  and  right-lateral  strike-slip  motion  in  the  Eastern  California  Shear  Zone  (ECSZ)  from  Dokka  and  Travis  [1990a|.  .All 
confidence  regions  (ellipses  and  quadrilateral)  are  95%. 

velocities  in  Figure  11b.  In  Figure  1  lb,  all  of  the  velocilies  arc  The.sc  mismatches  cannot  easily  be  corrected  by  fine  tuning  the 

plotted  relative  to  an  origin  fixed  at  VNDN.  The  gray  ellip.ses  arc  fault  model,  for  example,  by  decreasing  the  locking  depth  on  the 

the  observed  geodetic  velocities  for  a  series  of  stations  in  S.AF.  As  described  above,  an  additional  fault  west  and/or  .south  of 

California.  The  solid  dots  connected  by  a  dashed  line  represent  MADC  seems  necessary  to  bring  the  elastic  model  into  accord 

the  velocities  of  the  .same  sites  predicted  by  the  reference  elastic  with  the  observations.  But  in  the  absence  of  known  geological 

dislocation  model.  The  solid  lines  connecting  the  ob.scrvcd  and  .structures  accumulating  strike-slip  motion  between  MADC  and 

modeled  vectors  are  residual  velocities.  They  arc  .small,  but  VNDN,  either  structures  offshore  of  VNDN  or  time-dependent 

significant.  The  velocities  calculated  from  the  reference  model  all  effects  are  also  plausible  explanations. 

fall  outside  the  95%  confidence  ellipses  for  the  observed  valuc.s.  The  left-most  solid  dot  on  the  plot  represents  the  velocity  of 
The  fault-parallel  and  fault-normal  components  of  the  residuals  the  “Pacific  plate",  relative  to  VNDN.  predicted  by  our  reference 

are  similar  in  magnitude  and  increase  to  the  east.  In  addition,  the  elastic  dislocation  model.  The  reference  model,  using  a  locking 

differences  in  observed  velocity  betw'een  LOSP-VNDN.  MADC-  depth  of  25  km  on  the  central  section  of  the  SAF,  predicts  almo.st 

LOSP,  and  FIBR-MADC  arc  larger  than  predicted  by  the  model.  1  mm/yr  of  elastic  strain  accumulation  offshore  of  VNDN. 
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Vector  sums  w.r.t.  VNDN 


East  (mm/yr) 


Fig.  11.  (continued) 

Finally,  the  small  ellipse  labeled  ‘'North  America”  represents  the  plot  relative  to  VNDN.  It  is  interesting  that  our  elastic  disicxration 

observed  velocity  of  the  ensemble  of  sites  defining  the  North  model,  which  includes  deep  slip  on  the  SAF,  predicts  a  velocity  at 

American  plate  with  respect  to  VNDN.  It  is  compatible  with  the  OVRO  that  is  comparable  in  rate  (31  mm/yr)  but  significantly 

prediction  of  the  rescaled  NUVEL-1  plate  motion  model  (large  different  in  orientation  (S47E'’)  from  the  geologic  rate  for  the 

dark  gray  ellipse),  as  described  above.  We  note,  however,  that  the  SAF.  This  mismatch  is  the  result  of  the  effects  of  elastic  strain 

NUVEL-I  uncertainties  are  large  enough  to  include  the  possi-  caused  by  the  change  of  strike  of  the  SAF  in  the  Big  Bend  region, 

bility  of  1-2  mm/yr  of  deformation  offshore  of  VNDN,  which  we  as  well  as  motion  on  the  Oarlock  fault. 

have  neglected  in  this  comparison.  There  arc  two  geological  estimates  for  the  velocity  to  the  east 

A  third  constraint  on  how  the  deformation  is  partitioned  comes  of  our  network.  The  first  applies  to  a  path  leaving  OVRO  to  the 

from  geological  estimates,  which  we  summarize  by  constructing  northeast  and  repre.sents  the  extensional  deformation  in  the  Basin 

an  integral  along  a  path  crossing  the  plate  boundary  region.  and  Range.  Relative  to  North  America,  which  we  plot  as  the 

Where  the  path  crosses  an  active  structure,  the  geologically  rescaled  NUVEL-1  prediction,  the  Basin  and  Range  moves 

inferred  displacement  rale  and  the  appropriate  strike  arc  added  as  8.0  ±  3.3  mm/yr  at  N64°W  *  10°  according  to  the  geologic  data  in 

vector  components  to  the  sum.  Motions  due  to  block  rotations  ‘‘model  A”  of  Minster  and  Jordan  [1987].  The  confidence  ellip.se 

should  also  be  included,  although  they  arc  difficult  to  determine  for  this  geologic  estimate  docs  not  overlap  the  geodetically 

from  the  geologic  record.  The  resulting  path  integral  gives  the  oKserved  vclix'ity  for  OVRO.  This  remaining  bit  of  “discrepancy” 

relative  velocity  between  the  endpoints  of  the  path  [Minster  and  •  can.be  eliminated  with  a  more  northerly  azimuth  for  the  Basin 
./(Wu/i.  1984).  The  difference  betw’cen  the  path  integral  cro.ssing  and  Range  vector,  as  indicated  by  a  VLBI  estimate  for  the 

a  plate  boundary  and  the  relative  velocity  vector  between  the  two  velocity  of  the  Sierra  Nevada  with  respect  to  North  America 

plates  represents  the  discrepancy  vector.  The  geologic  velocity  [,4rg«.v  and  Cordon.  1991],  |i  could  also  be  eliminated  by 

estimates  arc  plotted  as  open  arrows  in  Figure  11b.  including  the  fault-normal  compression  we  have  measured  near 

The  conventional  geological  slip  rale  for  the  SAF  is  the  SAF. 

34  ±  3  mm/yr  at  S41  ±  2°E  [Minster  and  Jordan,  1984).  which  we  The  second  geological  estimate  involves  a  path  to  the  south  of 
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OVRO.  crossing  ihe  Easlcrn  California  Shear  Zone  (ECSZ) 
[Dokka  and  Travis,  1990a].  The  geological  eslimale  of  the 
average  slip  rale  across  this  structure  is  6-12  mm/yr  [Dokka  and 
Travis,  1990b],  which  we  take  to  be  a  95%  conlldencc  interval. 
Assuming  a  strike  of  N35'’W±  10°  by  an  eyeball  fit,  this  yields  a 
geologic  vector  quite  similar  to  the  geodetic  estimate  of  8  mm/yr 
at  N35°W  [.9fli’i3ge  ct  al.,  1990].  We  plot  the  former  quantity  as 
ECSZ  relative  to  the  rescaled  NUVEL-1  prediction  for  North 
Amc'rica.  The  95%  confidence  ellipse  of  the  ECSZ  geologic 
velocity  overlaps  the  observed  geodetic  velocity  of  OVRO,  ia  all 
the  geological  deformation  along  this  path  seems  to  be  accounted 
for  by  the  geodetic  ob.servalions. 

The  reference  elastic  model,  however,  is  missing  some  of  the 
motion  captured  by-  the  geodetic  and  geologic  estimates.  At 
OVRO.  for  example,  the  velocity  predicted  by  the  reference 
model  is  7-9  mm/yr  short  of  the  gcodelically  observed  value  and 
the  ECSZ  vector.  In  this  area,  the  reference  model  orcla.stic  strain 
accumulation  is  oversimplined.  It  includes  neither  the  ECSZ  nor 
the  Basin  and  Range  because  we  have  insufficient  information  to 
identify  the  active  structures  near  OVRO.  On  the  other  hand,  we 
can  also  attribute  some  of  the  model’s  shortcoming  to  .strike-slip 
motion  on  structures  outboard  of  MADC  and  to  compressive 
structures  in  the  vicinity  of  the  SAP.  neither  of  which  is  included 
in  the  simple  reference  model.  Within  the  context  of  elastic 
disIcKation  models,  motion  on  these  structures  is  sufficient  to 
explain  the  modified  San  Andreas  discrepancy.  We  leave  more 
quantitative  modeling  using  simple  dislocation  models,  as  well  as 
time-dependent  models,  to  future  papers. 


Conclusions 

We  have  developed  a  new  technique  to  rigorously  combine  7 
years  of  VLBI  and  GPS  measurements  in  a  single  solution.  With 
it.  w'e  have  mapped  the  velocity  field  in  central  and  southern 
California  with  a  precision  and  accuracy  as  good  as  2-3  mm/yr, 
sufficient  to  distinguish  among  some  tectonic  models. 

The  dominant  signal  in  the  ob.served  velocity  field  is  the  shear 
due  to  deep  slip  on  the  San  Andreas  and  other  known  faults.  After 
subtraction  of  a  conventional  model  for  the  strain  accumulation, 
how'ever,  significant  shear  remains  in  the  residual  field  through¬ 
out  the  region  of  our  network.  This  could  be  due  to  a  combination 
of  an  inappropriate  choice  of  model  parameters,  an  oversimplified 
model,  or  strain  accumulation  due  to  the  unmodeled  .structures 
such  as  the  San  Gregorio-Hosgri  fault.  Distinguishing  among 
these  possibilities  .should  be  a  high  priority  for  future  work. 

In  addition  to  unexplained  shear,  the  residual  velocity  field 
exhibits  significant  compressive  strain.  The  Ventura  basin,  Santa 
Barbara  Channel.  Los  Angeles  basin,  and  Santa  Maria  Fold  and 
Thrust  Belt  are  also  undergoing  active  shortening  at  rates  of  up  to 
5*1.5±1.5sl.  and  2  *  1  mm/yr.  respectively.  These  rates  are 
consistent  with  the  geological  estimates,  hut  tend  to  fall  at  the  low 
end  ol  the  allowable  range.  In  the  southern  Coast  Ranges,  our 
results  corroborate  the  conclusion  o(  Slien  and  . I acLson  [1993] 
that  the  rate  of  compression  normal  to  the  trace  of  the  San 
Andreas  lault  is  less  than  has  been  proposed  on  the  basis  of  trilat- 
eration  measurements.  We  also  find  significant  compression  to 
the  norlhea.st  of  the  SAP. 

The  rales  ol  rotation  about  vertical  axes  in  the  residual  geodetic 
wh'cily  field  diller  by  up  to  a  factor  of  2  from  those  inferred 
from  paleomagnclic  declinations.  The  reason  for  the  difference 


probably  involves  deformation  which  is  nonuniform,  either  in 
.spaccTir  in  time. 

Our  estimates  indicate  that  the  "San  Andreas  discrepancy’’  can 
be  rc.so!ved  to  within  the  geological  uncertainlies_oL3  mm/yr  by 
accounting  for  deformation  in  California  between  Vandenberg 
and  the  westernmost  Basin  and  Range.  Some  1-2  mm/yr  of  strain 
accumulation  on  structures  offshore  of  Vandenberg  is  also 
allowed  by  the  uncertainties.  South  of  the  Transverse  ranges,  the 
accounting  must  include  .some  5  mm/yr  of  motion  between  the 
oflshore  islands  and  the  mainland.  The  short-term  residual 
velocity  field  estimated  from  the  geodetic  measurements  is  in 
good  agreement  with  available  long-term  geologic  rates  to  within 
their  uncertainties.  We  therefore  have  no  compelling  evidence  for 
invoking  changes  in  the  rale  of  deformation  over  recent  geologic 
time. 


Apphndi.x:  Data  Analysis 
Estimation  Strategy 

Single-day  analysis  of  the  VI.BTobservdtidns.  The  VLBI  group 
delay  observations  are  analyzed  using  the  CALC/SOLVK  soft¬ 
ware  developed  at  the  NASA  Goddard  Space  Flight  Center  [e.g.. 
Ryan  et  al.,  1993]  and  the  Harvard-Smithsonian  Center  for 
.Astrophysics  [Herring  et  al.,  1981,  1986.  1990].  For  each  day,  we 
use  a  Kalman  filler  to  estimate  station  coordinates,  corrections  to 
lERS  Bulletin  B  values  for  pole  position  and  the  Earth's  axial 
rotation  (UTl-AT).  corrections  to  the  1AU-19S0  nutation  .series, 
jnd  .stochastic  variations  in  the  atmospheric  delay  and  reference 
clock  at  each  site  [Herring  et  al.,  1990].  To  avoid  distorting  the 
results  of  our  multisession  analysis,  we  perform  Ihe  single-day 
solutions  with  weak  constraints  on  all  of  the  estimated 
parameters.  The  VLBI  analysis  uses  760.521  group  delay 
measurements  in  1618  experiments,  with  typical  root-mean- 
.square  (rms)  scatter  of  30  ps  (-10  mm)  for  each  experiment. 

Single-day  analysis  of  the  CiPS  observations.  For  each  session, 
we  use  the  GAM  IT  .software  [Hock  et  ah,  1986:  King  and  Bock, 
1993]  to  analyze  carrier  beat  phase  observations.  Satellite  orbits 
are  obtained  by  numerically  integrating  the  initial  conditions 
u.sing  a  force  model  which  includes  lunar  and  solar  perturbations, 
the  GEM  L2  gravity  model  [I.erch  et  al.,  1985]  truncated  to 
degree  and  order  8,  and  a  three-parameter  model  for  non- 
graviialional  forces.  Using  these  satellite  ephemerides,  along  with 
nominal  values  for  the  site  coordinates,  values  of  Earth  rotation 
parameters  determined  from  the  VLBI  ob.serx'ations.  and  standard 
expressions  for  the  precession  and  nutation.  GAM  IT  computes 
theoretical  values  for  the  carrier  phase  observations  at  both  the  LI 
and  L2  frequencies  for  each  station-satellite  combination.  These 
theoretical  values  are  then  subtracted  from  the  ob.serx-ed  values  to 
form  phase  residuals  and  combined  as  double  differences 
(between  satellites  and  between  sites)  in  a  least  squares  anaivsis 
to  estimate  station  coordinates,  orbital  initial  conditions,  and 
phase  ambiguity  parameters  using  the  algorithm  described  by 
Schaffrin  and  Bock  [1988].  Our  "observable"  is  the  (doubly 
differenced)  linear  combmati.onjLC)  of  LJ  and  L2  that  eliminates 
that  part  of  the  ionospheric  delay  inversely  proportional  to  the 
square  of  the  frequency  [e.g..  Bock  et  al..  19S6].  .Atmospheric 
effects  arc  treated  by  estimating  the  deviation  al  zenith  from  the 
model  ol  Saastemoinen  [1972].  assuming  a  single  set  of  nominal 
meteorological  conditions  (2()'C.  1013.25  mbar.  50%  relative 
humidity)  at  sea  level,  but  accounting  lor  elevation  differences 
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among  ihe  slalions.  A  fuller  dcscriplion  of  the  Ihcorclical  model  is 
given  by  /•'c/'g/  [1991]  and  Murray  [1991],  . 

As  is  ihe  case  for  VLBI,  Ihe  product  of  Ihe  single-session 
analysis  is  a  sei  of  loosely  consirained  parameler  cslimalcs  and 
Iheir  covariance  malrix.  Before  generating  such  a  “loose” 
solution,  however,  we  must  first  resolve  the  phase  ambiguities. 
This  is  accomplished  in  a  series  of  four  intermediate  solutions, 
performed  automatically  and  sequentially  by  GAMIT; 

1.  All  parameters  arc  estimated  using  the  LC  combination  of 
the  Lt  and  L2  phases.  A  reference  frame  is  dcllned  by  imposing 
light  constraints  on  the  station  coordinates,  as  described  below. 

2.  With  the  geodetic  parameters  held  fixed  at  the  values 
obtained  in  solution  1.  the  wide-lane  (L2-L1)  ambiguity 
parameters  are  estimated  and,  if  possible,  constrained  to  integer 
values.  This  solution  employs  the  LI  and  L2  phases  under  an 
ionospheric  constraint  that  varies  from  1  to  8  parts  in  10^'.  If 
precise  pseudoranges  at  both  the  LI  .and  L2  frequencies  are 
available,  they  are  used  to  help  determine  the  wide-lane 
ambiguities. 

3.  With  the  wide-lane  ambiguities  held  fixed  at  the  values 
obtained  from  solution  2,  the  narrow-lane  (LI)  ambiguity 
parameters  are  estimated,  along  with  all  the  geodetic  parameters, 
from  the  LC  data.  As  in  .solution  1,  a  reference  frame  is  defined 
by  imposing  tight  constraints  on  the  station  coordinates. 

4.  With  the  resolved  values  of  both  the  wide-lane  and  narrow'- 


lane  ambiguities  held  fixed  to  the  integer  values  obtained  from 
solinions'2  and  3,  the  geodetic  parameters  are  estimated  from  Ihe 
LC  data. 

Thc.se  four  solutions  are  intermediate  slcp.s,  performed  with  the 
goal  of  xc.solvjng  the  integer- pha.se  ambiguities.  The  algorithms 
are  those  oi'  Doug  and  Hock  [1989],  except  that  the  values  of  the 
geodetic  parameters  estimated  in  solution  1  arc  retained  in 
solution  2.  This,  four-step  procedure  allowsms  to  resolve  more 
than  90%  of  the  ambiguities  in  most  experiments. 

Two  additional  solutions  are  neeeded  to  produce  the  parameter- 
estimates  and  covariance  matrices  for  later  input  into  Ihe  .second 
(mulli.scssion)  stage  of  our  analysis.  They  arc  performed  with 
orbital  and  positional  constraints  that  are  sufficiently  tight  to 
avoid  numerical  singularity,  but  sufficiently  .loose  (1  part  in  lO-"') 
to  avoid  contaminating  the  velocities  to  be  estimated  later.' 

5.  The  LC  solution  1  is  repeated  via  manipulation  of  the  least 
.squares  normal  equations  but  with  loose  constraints  on  the 
geodetic  parameters.  The  ambiguity  parameters  arc  free  to  assume 
real  values. 

6.  Solution  4  with  the  resolved  phase  ambiguities  con.strained 
to  jrticger  values  is  rcpeated.wiih-lQOse-eonstraints  on  the  geodetic 
parameters. 

To  resolve  the  ambiguities  requires  a  well-defined  reference 
frame,  which  we  established  by  constraining  the  station 
coordinates  in  solutions  1. 3.  and  4.  The  values  and  uncertainties 


Fig.  Al.  CPS  pha.se  residuals  with  different  noise  signatures.  The  plotted  quantity  is  Ihe  residual  (obsersed  minus  caleulaled) 
doubly  differenced  ionosphere-free  ("LC”)  combination  of  the  LI  and  L2  carrier  phase,  (a)  For  a  shon  (IS1  km)  line  between 
BLHL  and  CENT,  the  rms  scatter  is  about  O.O.td  cycles,  or  S  mm.  (/>)  For  a  longer  (.tU04  km)  line  betw  een  MOJM  and  WSFM, 
w  iih  more  mullipalhing.  Ihe  rms  is  larger,  about  0,094  cycles  (21  mm).  Both  plots  show  data  from  satellites  PRN  6  and  PRN  1 1 
ci'ilecled  on  March  29,  1990. 
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lor  ihcse  constraints  were  obtained  from  a  previous  multisession 
analysi^of  the  loose  solution  5  with  unresolved  arnbiguities.  Wc 
applied  constraints  to  the  coordinates  of  all  stations,  whether  in 
California  or  elsewhere,  thereby  minimizing  the  effect  of  an  error 
in  the  coordinates  for  any  one  station.  For  the  stations  in 
California,  we  applied  constraints  of  10  mm  for  the  horizontal 
coordinates  and  50  mm  for  the  vertical  coordinates;  for  the  other 
stations,  we  used  50  mm  for  the  horizontal  coordinates  and  100 
mm  lor  the  vertical  coordinates.  In  most  cases,  these  values  are 
several. times  larger  than  the  formaruheertainlies  obtained  from 
our  multisession  analysis.  We  repeat,  however,  that  the.se  tight 
constraints  were  only  impo.sed  on  the  intermediate  solutions  1,  3, 
and  4.  They  were  not  applied  to  the  .solutions  5  and  6  passed  to 
the  multisession  analysis  for  estimating  velocities.  In  some  cases, 
our  procedure  yielded  incorrect  values  for  the  ambiguities,  as 
indicated  by  elevated  values  of  X"  and/or  large  differences  in 
position  between  solutions  5  and  6.  In. these  cases  (21  of  113 
sessions),  we  took  the  con.servative  approach,  and  used  the  real- 
valued  solution  5  for  the  subsequent  estimation  of  station 
velocities  in  the  multisession  analysis. 

A  result  of  the  single-day  analysis  is  shown  in  Figure  Al. 
which  provides  two  examples  of  doubly  differenced,  ionosphere- 
free  phase  residuals  for  March  29,  1990.  The  residuals  between 
the  pair  of  field  sites  (Black  Hill  and  Center,  181  km  apart)  have 
an  rms  scatter  of  0.05  cycles  (10  mm),  while  the  scatter  for 


Mojave  an^Westford  is  4  times  larger.  The  difference  is  probably 
due  to  larger  unmodelcd  atmospheric  fluctuations  and  to  signal 
multipathing  from  highly  rcllective  structures  at  the  VLBI 
observatories.  The  unmodcled  excursions  in  phase  evident  at  the 
trackings  sites  can  systematically  a  ffecTfhe'estrhi^csbf  velocity, 
since  the  antennas  used  at  these  sites,  and  the  sites  themselves, 
have  changed  over  time  (Table  3). 

An  important  feature  in  the  residuafs'is  that  they'are  temporally 
correlated.  We  accounted  for  these  correlations  by  assigning  large 
a  priori  uncertainties  lo‘  ihe  30-s  samples  of  the  pha.se 
observations.  The  uncertainties  arc  10  mm  fdrlhc  LI  and  L2 
carrier  phase  measurements,  equivalent  to  64  mm  in  the  doubly 
differenced  ionosphere-free  combination.  In  a  slightly  different 
approach,  /.arson  and  Agnew  [1991]  accounted  for  the  correla¬ 
tions  by  compressing  the  daia  mto  6-min  normal  points.  This  is 
more  elficient  but  may  not  completely  average  out  the  short- 
period  effects  of  multipath.  The  scatter  in  the  estimates  of  relative 
station  positions  indicates  that  the  appropriate  uncertainty  for  6- 
min  .samples  of  the  doubly  differenced,  iono.sphere-frce  combina¬ 
tion  of  the  phase  measurements  is  about  0.1  cycles  (20  mm) 
[Lars'ffddndA^new.  1 99 1  [il-pivi^n'tTo'n.^S  c-ycles  (7()  mm)  for 
30-s  samples.  With  this  data  weighting,  the  scaled  uncertainties 
lor  positions  estimated  from  a  single  session  are  about  4  mm  in 
the  north,  6  mm  in  the  east,  and  20  mm  in  the  vertical 
coordinates,  provided  that  good  global  tracking  coverage  is 


TABLE  At.  Parameter  Constraints  Used  in  the  Multisession  Analyses 
- ??.t^tTie]er _ _ _  _ A  Priori _ _ _ A  Poste'rori  ** _ Stochastic  ^ 


Individual  California  CPS  Experiments 

Site  positions 

10  m 

1  m 

none 

Satellite  positions 

lOOm 

I  m 

0.1  m/d 

Satellite  velocities 

10  mm/sec 

0.1  mm/s 

0.01  mm  s“’  d“' 

Eclipsed  satellite  positions 

lOOm 

1  m 

10  m/d 

Eclipsed  satellite  velocity  c 

10  mm/sec 

O.I  mm/s 

10  mm  s“i  d“* 

Radialion  parameters 

100% 

1% 

1%/d 

Earth  orientation 

not  estimated 

VLB/  Experiments 

Site  positions 

lOm 

I  m 

none 

Site  velocities 

1  m/yr 

O.I  m/vr 

none 

Radio  source  positions 

100  mas 

0.1  mas 

none 

Earth  orientation*^ 

I(X)  mas 

20  mas 

1 0  mas/d 

Ctlobal  GPS  and  Combined  California  ilxpcrimcnts 

Site  positions 

10  m 

1  m 

none 

Site  velocities 

1  m/vr 

0.1  m/vr 

none 

Satellite  positions 

removed*-' 

Satellite  velocity 

removedx' 

Radiation  parameters 

removed*' 

Earth  orientation 

1 00  mas 

20  mas 

1 0  mas/d 

VLBI  and  GPS  Combination 

Site  positions 

10  m 

I  m 

none 

Site  velocities 

1  m/yr 

0.1  m/vr 

none 

Earth  orientation 

100  mas 

20  mas 

1 0  mas/d 

“  Typical  values.  In  many  cases,  these  values  are  large  because  no  coordinate  system  is  imposed  on  the 
.solution.  Quantities  w  hich  arc  invariant  to  changes  in  the  coordinate  .system,  such  as  length  and  analc.  are 
typically  better  determined  by  two  orders  of  magnitude. 

All  stochastic  parameters  are  modeled  as  random  walks  for  which  the  noise  in  the  process  increases  as  the 

square  root  of  time  (see  ! /erring  ct  ai,  [  1900|  for  dicussion)  _ 

s  This  level  of  constraint  is  fairly  loose  and  almost  uncouples  one  day's  orbit  from  the  next.  A  daily 
perturbation  ol  10  m  is  between  10  and  UKl  limes  larger  than  the  level  we  expect  for  most  satellites  [Herring 
Cl  «/..  10P1 1.  It  helps  to  avoid  Ihe  effect  of  unmodcled  nongraviialional  forces  on  Ihe  satellite,  especially 
w  hen  it  spends  part  of  its  orbit  in  Earth’s  shadow  (in  “eclipse’’). 

Earth  orientation  values  are  given  in  milliarc  seconds  (mas). 

‘  For  these  cases,  the  rows  and  columns  ol  the  covariance  matrix  for  these  parameters  are  remosed  belbrc 
incorporation  cd  the  solution  in  the  Kalman  liller.  This  procedure  is  equivalent  to  allow  ing  the  parameter  to 
adjust  Ireely  w  ilh  no  further  constraints  placed  on  its  value. 
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available,  and  phase  ambiguities  are  resolved  [Larson  and  Agnew, 
1 00 1 ;  A/urray,  1 9Q 1  ] . 

Multisession  analysis.  Wc  combine  the  single-day  VLSI  and 
GPS  solutions  using  a  Kalman  filter  program  (GLOBK) 
developed  for  the  combination  of  VLBI  experiments  [Herring  et 
a/.,  1900;  Herring,  1993;  Dong,  1993].  Our  combined  .solution 
contains  the  positions  and  velocities  of  212  stations  obtained  from 
1797  sessions  which  collected  a  total  of  -5,300,000  VLBI  group 
delay  and  GPS  doubly  differenced  carrier  phase  measurements. 

Because  each  single-day  solution  is  generated  with  loose 
constraints  on  all  of  the  parameters,  we  are  free  in  the  multi¬ 
session  analysis  to  impose  constraints  on  any  of  the  parameters. 
Imposing  these  constraints  is  a  rapid  procedure  with  the  GLOBK 
Kalman  filter,  allowing  us  to  test  the  effects  of  different 
constraints  efficiently.  For  example,  we  can  constrain  the  initial 
conditions  and  nongravitational  force  parameters  of  the  GPS 
satellites  to  define  a  single  dynamical  are  over  several  days  or  to 
allow  the  orbit  to  vary  stochastically  from  day  to  day. 

Table  A1  lists  the  a  priori  and  typical  a  posteriori  values  of  the 
uncertainties  for  all  parameters  estimated  in  the  multisession 
analyses.  For  most  parameters,  the  a  priori  constraints  were  suffi¬ 
ciently  loose  that  they  imposed  no  significant  bias  on  the 


estimates.  The  most  notable  exceptions  are  for  GPS  orbital 
parameters.  For  satellites  hot  eclipsed  during  an  observation 
period,  we  treated  the  satellite  state  vector  (initial  Cartesian 
position  and  velocity)  as  a  random  walk  with  the  stochastic 
prcxre.ss  noise  .set  .so  that  the  position  and  velocity  can  change  each 
day  by  10  cm  and  0.01  mm/s.  respectively.  Solar  radiation 
pressure  scale  factors  (direct  and  Y-bias)  were  also  tightly 
constrained  to  allow  1%  changesper^yrTor  eclipsing  .satellites, 
the  constraints  on  the  state  vectors  were  lixisened  to  10  m/d  and 
10  mm  s"’  d"’,  effectively  decoupling  one  day's  orbit  from  the 

next.  - — -  - - - 

The  strongest  constraints,  however,  were  placed  on  the  station 
positions  and  velocities.  We  varied  the  station  constraints  (Table 
Al)  to  [>erform  three  different  solutions,  which  we  u.sed  for 
different  purposes:  (1)  to  define  a  terrestrial  reference  frame.  (2) 
to  establish  a  set  of  a  priori  coordinates  sufficiently  consistent  for 
resolution  of  GPS  pha.se  ambiguities,  and  (3)  to  provide  the  most 
robu.st  estimate  of  the  station  velocities.  For  succinctness,  we  call 
these  solutions  “frame,”  “coordinate,”  and  “velocity”  solutions, 
respectively. 

Theffirst  step  in  the  frame  solution- used-only  VLBI- data  to 
define  a  terrestrial  reference  frame.  Because  of  the  weak 


Fig.  A2.  Evolution  with  time  of  the  (o)  north,  (b)  east,  and  (c)  venical  components  of  the  relative  position  \ector  between 
Mojave  (MOJA)  and  Vandenberg  (VNDN)  as  determined  by  CPS  (solid  circles)  and  VLBI  (open  squares).  The  velocity 
estimates  from  the  coordinate  .solution  have  been  removed  (2>S.60,  -25.77,  and  2.99  mm/yr  respectively)  and  are  represented  by 
the  hori7ontal  solid  line.  There  have  been  no  offsets  introduced  into  the  horizontal  coordinates.  As  discus.sed  in  the  text,  vertical 
discontinuities  corresponding  to  changes  in  antennas  al  Mojave  have  been  removed  from  the  CPS  results  at  1986.9,  1988. 6,  and 
1 989.5. 
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constraints  imposed  on  the  VLB!  solution,  the  reference  frame  for 
the  VLBI  analysis  is  initialiy  only  weakly  constrained.  The 
weakly  defined  frame  has  effectively  a  twelfth-order  rank 
dcllcicncy;  three  translations,  three  rotations,  three  translation 
rates,  and  three  rotation  rates.  To  re.solve  the  rank  detlency  in  rate, 
we  determined  the  three  translation  and  three  rotation  rates 
between  the  weakly  constrained  frame  and  our  realization  of  a 
North  American  fixed  reference  frame  by  minimizing  the 
horizontal  velocities  in  the  North  American  frame  of  10  North 
American  VLBI  sites  (Haystack.  Westford.  Richmond.  Gilmore 
Creek,  Nome,  Yuma,  Greenbank.  Pictown,  Flagstaff,  and 
Platicville).  To  resolve  the  rank  defiency  in  position,  we 
estimated  three  rotation  and  three  translation  parameters  bv 
minimizing  the  differences  between  the  horizontal  coordinates  of 
the  same  10  VLBI  sites  in  our  frame  and  the  ITRF-IWl  reference 
Irame.  This  latter  step  does  not  affect  the  velocity  solution. 

In  the  next  step,  combining  the  VLBI  and  CPS  analyses,  we 
adopted  this  reference  frame  by  fixing  the  horizontal  velocities  of 
three  sites  (Fairbanks.  Westford.  and  Wettzcll)  to  their  VLBI 
values  and  constraining  the  relative  positions  and  velocities  of 
one  or  more  VLBI  and  GPS  monuments  at  these  same  locations. 
This  realization  of  the  North  American  plate  is  independent  of 
any  global  plate  motion  model. 

For  the  coordinate  solution,  we  applied  constraints  on  the 
hc'rizontal  components  ol  the  station  position  and  velocity 
parameters.  This  imposed  the  “ties’’  (vectors)  between  the 
multiple  benchmarks  and  antenna  phase  centers  located  at  the 


same  site.  In  the  coordinate  solution,  ties  were  imposed  in  both 
horizontal  position  and  velocity.  The  ties  in  position  have  been 
e.stablished  by  terrestrial  surveys.  The  tics  in  velocity  allow  no 
relative  motion  between  neighboring  benchmarks.  For  example, 
we  assumed  that  the  five  different  monuments  at  Westford, 
Massachusetts,  arc  connected  to  the  same  rigid  bkx:k. 

For  the  velocity  solution,  wx  imposed  tics  in  horizontal 
velocity  but  not  position.  The  only  exception  was  for  Mojave, 
where  we  tied  the  horizontal  ccx>rdinales  of  the  VLBI  and  GPS 
antennas  together.  Although  avoiding  the  ties  in  position  omits 
some  information,  it  avoids  the  (unfortunately  likely)  possibility 
of  contaminating  the  estimates  of  the  station  velocities  with  an 
erroneous  tie.  The  velocity  field,  shown  in  Figure  2  and  Table  4. 
yields  the  relative  velocities  we  interpret  in  the  text.  As  discu.sscd 
below,  the  velocity  .solution  turns  out  to  be  quite  similar  to  the 
coctrdinate  .solution. 

Discussion  of  Geodetic  Quality 

To  a.s.se.ss  the  accuracy  ot  the  velocity  lield  given  in  Table  4. 
we  need  to  know  the  temporal  and  spatial  nature  of  the  error 
.spectra  of  pexsition  determinations.  In  general,  ihe.se  spectra  arc 
unkno^^■n  and  diltlcult  to  determine.  We  can.  however,  examine 
several  .statistics  from  our  analysis  in  an  attempt  to  assess  the 
likely  magnitudes  of  the  errors.  The  components  of  the  error 
spectra  most  ncce.ssary  for  a.s.sessing  the  velocity  field  uncertainty 
arc  those  with  temporal  scales  comparable  to  the  total  duration 
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Fig.  A.1.  Evolulion  with  lime  of  ihe  (a)  nonh.  (bj  east,  and  (c)  vertical  components  of  ihe  relative  position  vector  between 
Vandenberg  (VNDN)  and  Santa  Cruz  Island  (CENT).  The  velocity  estimates  from  the  coordinate  solution  have  been  removed 
(I.d4.  -2.64,  and  2.34  mm/yr,  respectively)  and  are  represented  by  the  horizontal  solid  line. 


(-5  years)  and  .separation  (-1  year)  of  our  experiment.s,  and  at 
.spatial  .scales  comparable  to  the  baseline  lengths  in  California 
(50-500  km).  Since  there  is  no  redundancy  in  the  measurements 
at  these  time  intervals,  we  are  forced  to  extrapolate  from  .statistics 
for  shorter  time  intervals.  In  doing  so,  we  assume  that  the  noise 
proces.ses  are  stationary.  This  assumption  is  difficult  to  ju.stify 
because  there  has  been  an  evolulion  of  both  the  number  and  type 
of  receivers  and  satellites  over  the  duration  of  our  data  set.  There 
is  greater  spatial  than  temporal  redundancy  in  our  data,  hut  the 
most  reliable  statistics  are  computed  using  the  entire  network. 
The  assumption  of  .spatial  homogeneity  is  also  difficult  to  justify 
since  the  velocities  between  sites  separated  by  larger  distances  are 
more  likely  to  be  affected  by  the  projection,  through  the  orbiLs,  of 
the  error  at  the  tracking  sites  than  velocities  between  clo.scly 
spaced  sites. 

Short-  and  long-term  precision.  In  order  to  assess  the  precision 
of  our  estimates  of  relative  station  position,  wc  evaluated,  their 
scatter  from  day  to  day  and  year  to  year.  To  estimate  the  .scatter, 
we  performed  a  .scries  of  multisession  analyses  in  which  we 
allowed  the  coordinates  of  California  stations  to  vary  stochasti¬ 
cally  from  day  to  day  or  experiment  to  experiment.  In  lhe.se 
analytics,  the  coordinates  of  all  stations  outside  California  were 
tightly  constrained  at  their  values  from  the  coordinate  solution. 
For  the  GPS  experiments,  a  Kalman  filter  back  .solution  was 


neces.sary  in  order  to  apply  the  appropriate  stochastic  constraints 
to  the  orbital  parameters. 

Figure  A2  shows  the  behavior  of  the  VLBI  and  GPS  estimates 
for  the  north,  cast,  and  vertical  coordinates  of  Mojave  with 
respect  to  Vandenberg,  the  most  frequently  observed,  and  one  of 
the  longest  California  baselines  (350  km).  The  weighted  means 
from  the  VLBI  and  CPS  estimates  differ  by  2.8  mm  in  the  north 
and  9.1  mm  in  the  east,  slightly  larger  than  wc  would  expect  from 
errors  in  the  terrestrial  survey  tics.  The  height  differences  are 
much  larger,  as  di.scu.ssed  below. 

To  quantify  the  scalier,  we  prc.scnt  two  staii.slics:  the  weighted 
root  mean  square  (wrms)  scatter  (often  called  repeatability),  the 
value  of  Vx^/  (where  x~ 's  the  chi-.square  statistic,  and  /is  the 
number  of  degrees  of  freedom).  In  the  north  component,  the  wrms 
about  the  weighted  mean  is  6.9  mm  for  150  VLBI  .sc.ssions  and 
12.0  mm  for  81  GPS  sessions.  In  the  ca.st  component,  the  wrms  is 

6.3  mm  lor  VLBI  and  18.0  mm  for  CPS;  in  height  27  mm  for 
VLBI  and  49  mm  for  GPS.  For  both  techniques.  Figure  A2  show's 
evidence  of  both  short-term  and  long-term  scatter  exceeding  the 
formal  uncertainties  of  the  estimates.  For  the  VLBI  sessions  the 
values  of  V^/  are  2.2  in  the  north.  2.5  in  the  east,  and  1.8  in  the 
height;  for  the  GPS  experiments  the  values  arc  4.0  in  the  north. 

3.3  in  the  east,  and  2.2  in  the  height.  These  values  are  higher  than 
normal  lor  GPS  and  VLBI  baselines,  probably  due  to  unmodeled 
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tropospheric  effects  al  Vandenberg,  a  coastal  site  subject  to 
atmospheric  gradients. 

Figure  A3  shows  the  behavior  of  the  estimated  coordinates  of 
site  Center  (CENT)  on  Santa  Cruz  Island,  with  respect  to 
Vandenberg,  101  km  distant.  For  41  sessions  of  data,  the  wrms 
scatter  is  5.0  mm  in  the  north  component,  6.8  mm  in  the  east,  and 
22.0  mm  in  height.  The  smaller  scatter  in  the  east  and  vertical 
components,  compared  with  the  Vandenberg— Mojave  baseline,  is 
typical  for  the  more  closely  spaced  sites  in  our  GPS  core  network. 
For  this  baseline,  phase  ambiguities  were  resolved  more 
consistently  than  for  the  longer  baselines. 

Table  A2  summarizes  the  stati.siics  measuring  the  short-term 
precision  of  the  GPS  and  the  long-term  precision  of  both  the  GPS 
and  VLBI  experiments.  For  each  GPS  experiment,  we  list  the 
'weighted  rms  and  Vz“/  for  the  north,  ea.st,  and  vertical 
components  ol  the  observed  baselines.  We  also  measure  the 
.scaiier’s  dependence  on  baseline  length  by  the  slope  with  which 
the  scatter  increases  as  a  fraction  of  interstation  distance.  For  all 
experiments  except  1  and  18,  the  short-term  scatters  are  Ic.ss  than 
-6  mm  in  the  north  component,  -10  mm  in  the  east,  and  -30  mm 
in  ihe  vertical,  with  average  values  of  4  mm  in  the  north.  7  mm  in 
the  cast,  and  20  mm  in  the  vertical.  Most  values  of  the  V^/ 
slalistic  for  the  individual  experiments  are  between  0.8  and  1.5, 
indicating  that  the  assigned  measurement  uncertainty  is 
appropriate  for  the  short-term  scatter  in  the  estimates.  For  the 
major  experiments  that  included  widely  spaced  sites,  the  ba.selinc 
dependent  term  of  the  scalier  is  less  than  1  part  in  10^  in  length. 


The  multisetfsion  analysis  has  produced  good  GPS  repeatability 
even  for  those  experiments  with  few  regional  sites  (7,  8,  9,  15,  16, 
17)  or  poor  global  tracking  (2,  3,  7,  8.  9).  The  only  experimenLs 
with  wrms  scatters  in  the  horizontal  coordinates  greater  than  10 
mm  were  1  and  18.  Experiment  1  (December  1986)  had  no  sites 
outside  California,  so  the  uncertainties  for  all  the  baseline 
estimates  are  anomalously  large.  The  scatter  in  Experiment  18 
may  have  been  increa.scd  by  incomplete  modeling,  for  receivers 
with  different  sampling  limes,  of  variations  of  the  satellite  clocks 
impo.scd  by  Selective  Availability  [Feigletai.  1991], 

The  la.st  two  entries  in  Table  A2  give  the  long-term  repeata¬ 
bilities  lor  both  GPS  and  VLBI.  For  VLBI  these  values  are  com¬ 
puted  from  the  scatters  in  each  24-hour  se.ssion.  The  values  of  the 
wrms  scatter  are  similar  to  those  obtained  for  the  Mojave- 
Vandenberg  ba.seline  (7  mm  for  the  north  component.  6  mm  for 
the  ea.st.  and  38  mm  for  the  vertical),  and  the  is  about  1.5. 
For  GPS.  we  computed  the  mean  values  and  their  uncertainties 
from  each  3-  to  5-day  experiment  and  u.sed  these  values  to  de¬ 
termine  the  long-term  .scatter  (Figure  A4).  The  wrms  values  are 
similar  to  the  averages  for  short-term  scatter  (4  mm  in  the  north.  6 
mm  in  the  ca.si,  and  23  mm  in  the  vertical),  but  the  values  of  the 
Vk*  /  stati.siic  arc  about  twice  those  for  the  individual  experi¬ 
ments.  This  suggests  that  there  is  a  systematic  comp<inent  of  the 
error  within  each  experiment  that  is  comparable  to  the  short-term 
.scatter  and  that  the  noi.se  spectrum  is  likely  to  be  red;  that  is.  the 
noi.se  power  is  larger  at  low  frequencies  than  al  high  frequencies. 
The  obtained  from  the  long-term  scatters  is  about  1.5  for 
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TABLE  A2.  Short-Term  and  Long-Term  Scalier  of  CPS  and  VLBI  Experiments 


Experiment 

Year 

Number  of 
Lines" 

North 

East 

Up 

wrms^, 

mm 

ppb 

wrms, 

mm 

ppb 

v^/ 

wrms, 

mm 

ppb 

^If 

00 

86.5 

118 

2.5 

-1 

1.0 

1.9 

4 

0.5 

1.3.4 

17 

0.6 

01 

87.0 

35 

8.0 

3.3 

1.2 

1X4 

34 

1.0 

31.6 

25 

1.5 

02 

87.1 

67 

4.6 

2 

0.9 

7.8 

26 

0.8 

26.2 

42 

1.2 

0.4 

87.1 

115 

2.9 

3 

0.6 

5.6 

0.7 

20.0- 

■nF 

0.5 

OT-- 

87.5 

15 

.3.2 

0 

0.7. 

9.9 

-I 

1.2 

26.7 

6 

0.7 

08 

87.8 

.76 

5.3 

3 

0.9 

7.8 

28- 

0.8. 

•33.7 

69 

-  0.9  - 

09 

87.9 

.38 

5.5 

43 

0.8 

4.3 

16 

_  0.6 

18.8 

45 

0.7 

10 

88.3 

46 

2.2 

0 

1.0 

.375 

2 

1.2 

8.8 

1 

0.5 

11 

88.3 

79 

.3.0 

1 

0.8 

4.1 

3 

0.8 

21.1 

3 

0.7 

I.V 

80.3 

36 

2.9 

1 

t.l 

6.4 

11 

1.7 

2.3.1 

-17 

1.2 

14 

89.3 

78 

2.9 

3 

1.6 

.3.7 

1 

1.5 

24.2 

-1 

1.4  - 

15 

89.3 

15 

3.5 

0 

1.6 

X7- 

0. 

1.0 

38.6 

-26 

2.6 

16 

89.3 

36 

.3.1 

5 

1.4 

6.8 

36 

2.2 

20.4 

-4 

.  1.1 

\FI 

90.1 

18 

4.7 

0 

1.5 

4.5 

-3 

1.2 

16.4 

-6 

1.2 

17e 

90.3 

•  29 

4.3 

1 

1.5 

7.0 

10 

1.6 

18.5 

15 

1.1 

ISe-f 

90.3 

210 

.3.6 

6 

1.0 

1.3.2 

37 

2.1 

24.5 

26 

1.0 

VBie 

90.5 

65 

4.5 

to 

2.3 

5.7 

14 

2.1 

26.0 

55 

1.8 

VF2 

90.8 

33 

4.5 

3 

1.8 

4.2 

■■  2 

1.4 

17.4 

49 

0.8 

20^- 

91.2 

1.36 

2.6 

1 

1.0 .  . 

.:=;6.9 

5 

1.2 

12.7 

0.8 

SBI 

91.5 

214 

5.5 

3 

0.6 

■  10.1 

19 

0.6 

23.3 

21  ' 

0.7' 

VB2 

92.5 

153 

4.3 

7 

2.2 

5.5 

16 

2.5 

24.0 

45 

1.8 

All  CPS 

156 

3.6 

2 

2.6 

5.8 

25 

2.5 

23.0 

-11 

2.4 

All  VLBI 

90 

7.1 

-1 

1.4 

5.5 

2 

L5_, 

.  3!U 

r-7 

1.3 

“  For  individual  cxperimenls,  all  baselines  less  than  500  km  in  lengih  and  which  were  measured  al  least  twice 
are  included  in  the  calculations  of  the  wrms  andVj?/.  For  the  last  two  entries,  at  least  three  measuremenLs 
arc  required. 

h  For  individual  experiments,  the  wrms  andVir”/  ‘U'f  computed  about  the  mean  length  for  each  session.  For 
the  last  two  rows,  the  rms  andv^/  are  computed  about  the  best  fitting  slraighl  line.  The  length-dependent 
part  of  the  wrms  is  given  in  parts  per  billion  (ppb)  of  baseline  length. 

For  these  sessions,  biases  are  not  resolved  correctly  on  some  days  (as  judged  by  anomalous  coordinates 
relative  lo  the  solution  with  real-valued  ambiguities)  and  the  solution  with  real-valued  ambiguities  is  used  in 
the  multisession  analyses.  This  occurrs  for  21  sessions,  while  the  biases  are  correctly  resolved  for  the  other 
1  LI  sessions. 

Only  2  of  the  4  days  of  this  experiment  are  included  in  the  analysis. 


North  East 


Length 


Fig.  A4.  Flislogram  of  long  term  scatter  for  the  length,  north,  e,ast.  and  vertical  components  of  the  GPS-eslimaled  inlerstation 
vectors.  The  values  represent  the  wrms  scatter  about  a  straight  line  c.stimaled  for  each  component.  Only  stations  in  California 
observed  over  an  interval  greater  than  1.5  years  are  included. 


VLBI  and  2.5  for  GPS;  and  the  V x"  /  of  the  combined  multi¬ 
session  solution  (which  accounts  rigorously  for  all  the 
correlations  among  the  estimated  quantities)  is  l.d.  These 
statistics  indicate  that  doubling  the  formal  standard  deviations 
from  the  combined  analysis  yields  realistic  uncertainties  if  the 
noi.se  spectrum  is  Hat  between  I  and  5  years.  Though  we  doubt 
that  this  is  strictly  true,  we  have  computed  the  statistic  at  the 


longc.st  po.ssible  sampling  times  while  averaging  as  much  of  the 
high--frcquency  noi.se  as  possible.  — 

Another  approach  to  calibrating  the  uncertainties  is  to  compare 
different  solutions.  For  example,  in  Figure  A5  we  show  the 
differences  in  the  estimated  velocities,  relative  to  Vandenberg, 
between  the  velocity  and  the  coordinate  .solutions.  The  largest 
difference  at  any  of  the  core  sites  is  2.5  mm/yr,  with  the  typical 
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dilTercnce  being  1  mm/yr.  Although  these  dilTerences  could  be 
due  to  an  erroneous  tie  at  any  of  the  common  VLBI  and  GPS 
sites,  they  could  also  indicate  large-scale  systematic  errors  in  the 
GPS  and  VLBI  solutions.  Such  .systematic  errors  might  be  due  to 
incomplete  dynamical  models  of  the  satellite  orbits,  deviations  of 
the  phase  centers  of  the  GPS  receivers,  and  deformation  of  the 
VLBI  antennas.  Doubling  of  the  formal  standard  deviations  helps 
to  account  tor  the.se  effects  in  our  uncertainties. 

As  an  additional  lest,  we  performed  an  analysis  in  which  the 
VLB!  and  GPS  velocities  in  California,  except  at  Mojave,  were 
separate  parameters.  Constraints  on  the  differential  position, 
where  available,  and  velocity  were  applied  at  all  sites  outside 
California.  The  velocities  estimated  at  eight  sites  with  both  VLBI 
and  GPS  observations  are  shown  relative  to  the  joint  GPS-fVLBI 
velocity  solution  in  Figure  A6.  The  dilTerences  in  velocity  esti¬ 
mated  by  the  three  solutions  are  between  1  and  5  mm/yfin  both 
horizontal  components.  Compared  to  the  joint  GPS+VLBI 
solution,  the  GPS-only  estimates  (open  arrows)  show  larger  dif- 
Icrenccs  in  the  cast  component  than  in  north,  while  the  opposite  is 
true  Urr  the  VLBI-only  estimates  (.solid  arrows).  The  differences 
are  signilicant  at  OVRO  and  V'NDN.  At  the  I'ther  six  stations,  the 
conlidence  ellipses  for  the  GPS-only  estimates,  calculated  at  95% 
eonlidence  alter  doubling  the  standard  deviations,  include  the 


VLBI-only  and  the  joint  GPS-VLBI  estimates.  This  comparison 
provides  independent  evidence  that  the  scale  factor  derived  from 
our  repeatability  analysis  is  not  overly  optimistic. 

hxiimales  of  ihe  vertical  component.  For  both  global  and 
regional  networks,  the  relative  vertical  coordinates  of  GPS 
■stations  are  less  accurate  than  the  horizontal  coordinates.  The 
formal  uncertainties  in  the  vertical  coordinates  are  about  a  factor 
of  2  larger  than  for  the  horizontal,  since  the  observations  can  be 
made  only  in  the  hemisphere  above  the  horizon  and  there  is  a 
high  correlation  between  e.stimates  of  atmospheric  parameters  and 
the  vertical  coordinates.  There  is  often  an  additional,  systematic, 
eompi'nent  to  vertical  error  introduced  by  phase  center  variations 
ol  the  antenna.  This  is  most  problematic  when  different  antennas 
arc  used  at  each  end  ol  the  ba.scline.  but  also  important  for  similar 
antennas  observing  at  large  spatial  separations.  Finally,  ties 
between  the  GPS  and  VLBI  antenna  phase  centers  and  ground 
less  certain  for  the  vertical  than  the 
horizontal  directions,  both  in  a  formal  .sense  and  in  the  likelihood 
ol  blunders  in  the  measurements  (see.  e.g..  appendix  A  of  Murray 
[19Q|]  for  a  di.scu.ssion  of  the  suiv-ey  at  Kokee  Park.  Hawaii). 

Our  attempt  to  incorporate  survey  ties  at  the  Mojave  site 
illustrates  the  potential  magnitude  of  the  \erlical  problem.  The 
June  19S6  (experiment  0)  observations  were  made  with  a  Tl  4100 
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Fig.  A6.  Compari.son  of  CPS  (open  arrows)  and  VLBI  (solid  arrow.s)  csiimales  of  slalion  velocities,  shown  with  respect  to  the 
combined  CPS+VLBI  velocity  solution  (ellipses  without  arrows),  in  which  no  ties  were  imposed  between  the  VLBI  and  CPS 
stations  in  California,  except  at  Mojave.  Ellipses  as  in  Figure  A5. 


antenna  .set  on  a  tripod  above  a  ground  monument  (NCMN  1983 
RM  1,  which  we  designate  MOJl).  From  November  1986  to  July 
1988  (experiments  1-1 1 ),  a  T1  4100  antenna  was  mounted  on  the 
corner  of  the  office  building.  In  July  1988.  it  was  replaced  by  a 
FRPA-2  antenna  (MOJF).  which  was  u.sed  for  experiment  14.  In 
February  1989.  a  MiniMac  receiver  and  antenna  were  installed 
nearby,  and  the  Tl  ob.servalions  were  di.scontinued  the  following 
July.  Our  experiments  in  March  and  April  1989  (experiments  13- 
16)  u.sed  both  the  TI/FRPA-2  and  Mini.Mac  observations.  All 
sub.sequent  experiments  used  only  the  MiniMac  ob.se rvations. 
Conventional  or  CPS  surveys  have  been  conducted  by  us.  by  the 
uses,  and  by  Bendix  Field  Engineering,  and  all  have  formal 
uncertainties  of  10  mm  or  less.  The  differences  among  our 
estimates  of  the  relative  heights  of  each  of  the  GPS  antennas  and 
the  (fixed)  VLBI  antenna  (in  the  sense  GPS-VLBI),  after 
applying  the  tie  information  are  19  mm  for  MOJl.  52  mm  for 
MOJA.  1 10  mm  for  MOJF.  and  159  mm  for  MOJM,  each  with  a 
(doubled)  formal  uncertainty  of  about  30  mm.  If  we  exclude 
MOJF  and  MOJM  from  the  comparison  (and  also  JPLl. 
Monument  Peak,  and  Santa  Paula,  which  had  only  a  few 
occupalions),  the  mean  difference  (GPS-VLBI)  in  height  for  the 
California  sites  is  60  mm.  with  all  of  Ihc  values  falling  between 
19  and  74  mm. 


A  similar  consistent  difference  (approximately  60  ±  30  mm)  is 
found  for  the  height  differences  estimated  for  Platlcville,  the  two 
sites  at  Algonquin,  and  the  three  GPS  antennas  used  at  Westford. 
This  consistency  suggests  a  systematic  bias  in  our  estimates  of 
height  for  all  of  the  California  CPS  cxperimenLs.  This  result  is  not 
surprising,  given  the  limited  coverage  of  the  tracking  network  for 
these  experiments  and  the  fact  that  we  did  not  constrain  any  of  the 
height  values.  Estimates  of  height  from  the  global  GPS 
experiments  do  not  show  this  bias.  Two  of  the  North  American 
.sites.  Richmond  and  Yellowknife,  as  well  as  Kokee  and  the  three 
European  sites,  have  much  larger  inconsisiencics  (-340  to  +143 
mm).  Most  of  these  inconsistencies  are  probably  due  to  the 
inlluencc  of  orbital  errors  and  a  changing  and  sometimes  weak 
netw’ork  of  continental-  and  global-scale  GPS  tracking  stations. 
Nevertheless,  we  cannot  rule  out  tic  blunders  in  every  case.  For 
this  reason  we  have  taken  a  conservative  approach  and  not 
appliceL-verlical  ties  between  VLBI  and  GPS  antennas  in  cither 
our  coordinate  or  velocity  .solution. 

Comparison  with  previous  (iPS  and  irilaleralion  results.  In 
January  1987  (experiment  3).  we  reoccupied  a  historical  trilatera- 
tion  network  spanning  the  Santa  Barbara  Channel  to  estimate 
strain  accumulation  by  comparing  line  lengths  determined  by  tri- 
lateration  and  by  CPS  [Larsen,  1991;  Larsen  el  al..  1993,  here- 
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of  the  velocily  solution  eslimaled  in  this  sliidv  ,witbjJwt  .of  LWL  [Urson  and  Webb,  mi-,  Urson,  this 
issue).  The  ploned  vectors  are  this  study  minus  LWL,  shown  relative  to  LACU.  We  show  95%  confidence  ellipses  fron^  the 
LWL  analysis  to  illustrate  the  relative  sizes  of  the  uncertainties  for  the  sites.  Correct  confidence  ellipses  for  the  differences 
between  the  two  analyses  are  difficult  to  compute  because  we  cannot  determine  reliably  the  degree  of  independence  of  the 
estimates,  which  is  a  function  of  the  data  used  and  the  a  priori  weight  assumed  for  both  data  and  parameters  \Done  IW!)  The 
two  data  sets  differ  by  several  CPS  experiments,  as  well  as  all  the  VLBI  data,  as  described  in  this  appendix.  We  also  use  CPS 
data  sampled  every  .^0  s.  rather  than  every  6  min.  as  used  by  LWU  but  the  effects  of  this  difference  are  partially  mitigated  by 
our  larger  assumed  measurement  uncertainly  and  by  correlations  between  the  .70-s  data  samples... 


liller  relerred  to  as  LAH],  This  network  was  reoccupied  again  in 
1991  (SBl),  with  results  reported  by  Larson  and  Webb  [1992] 
and  Larson  [this  issue]  (hereafter  referred  to  as  LWL).  As  part  of 
their  study,  they  present  an  independent  analysis,  using  the 
GIPSY  software  developed  at  JPL  [Lichten  and  Border,  1987; 
l.ic/iten,  1990],  of  most  of  the  data  we  collected  at  several  of  the 
core  sites  in  the  vicinity  of  the  Santa  Barbara  Channel. 

For  the  sites  common  to  both  studies,  we  calculate  the  differ¬ 
ences  between  the  velocities  estimated  by  our  joint  GPS+VLBI 
velocity  solution  and  that  estimated  by  LWL  using  GPS  data 
alone.  In  Figure  A7,  these  differences  are  shown  with  respect  to 
La  Cumbre  (LACU),  near  the  center  of  the  network  of  common 
stations.  The  difference  in  the  nonh  component  is  insignificantly 
different  from  zero  with  95%  confidence  for  all  sites.  For  the  east 
component,  the  difference  is  larger  (3.8  mm/yr  at  the  core  .station 
VNDN)  and  significant  at  four  stations  (VNDN.  SOU,  NIGU, 
and  BRSH).  The  consistent  sense  of  the  difference  in  the  east 
component  (toward  LACU  in  this  plot)  suggests  the  possibility  of 
a  systematic  bias,  which  we  attribute  to  problems  in  the  reference 
Irame  adopted  by  LWL  in  their  simpler.  GPS-only  strategy.  The 
problem  is  most  apparent  at  the  “noncore’’  stations  SOU,  NIGU. 
and  BRSH.  for  which  both  studies  must  rely  heavily  on  a  small 
number  (3,  4,  and  5.  respectively)  ol  occupations,  including  the 
early  experiment  0  in  1986.  when  the  number  of  global  tracking 
stations  was  quite  limited  (Table  3  and  Larson  et  al.  [1991]).  The 
most  significant  difference  between  the  two  strategies  occurs  at 
VNDN,  a  site  where  our  data  set  is  very  dilferent  from  the  one 
used  by  LWL.  At  VNDN.  our  data  .set  includes  150  VLBI  exper¬ 
iments  and  24  GPS  occupations  (in  experiments  15,  VFl.  and 
VF2).  w'hich  w'cre  not  analyzed  by  LWL,  while  their  data  .set  in¬ 
cludes  seven  USGS  occupations  of  the  MOJA-VNDN-OVRO 
triangle  not  included  in  our  analysis. 


While  this  .sensitivity  of  velocily  vectors  to  solution  .strategy  is 
disturbingly  large  for  regions  as  large  as  our  entire  network,  the 
clfects  are  less  serious  lor  relative  velocities  between  neighboring 
stations,  such  as  the  lines  spanning  the  Santa  Barbara  Channel. 
This  can  be  seen  by  comparing  the  rale  of  change  of  distance  dUdt 
between  stations  for  our  GPS  analysis  for  1986-1992,  for  the 
CPS  analysis  of  LWL  for  1986-1991.  and  for  the  comparison  of 
GPS  measurements  in  1987-1988  with  trilateration  measurements 
in  1970—1971  by  LAH.  Sites  LACU  and  GAVI  were  common  to 
all  three  studies,  while  two  trilateration  sites  (DEVL  and  CHAF) 
have  been  tied  to  more  convenient  GPS  sites  within  5  km  (CENT 
and  SOLI,  respectively).  In  our  compari.son.  we  assume  that  there 
is  no  motion  between  ihe.se  pairs  of  stations. 

The  rales  of  change  in  line  length  for  these  three  studies  are 
presented  in  Table  A3.  While  there  is  a  hint  of  more  east-west 
extension  in  the  LWL  results  than  in  ours,  there  is  no  significant 
dilference  (at  90%  confidence  level)  in  dlidi  between  the  sites 
common  to  our  and  the  LWL  analysis,  \^’ith  three  exceptions, 
there  arc  no  significant  differences  among  the  GPS-GPS  and 
GPS-EDM  line  length  changes.  The  discrepancies  arc  for  the 
lines  DEVL-GAVI.  where  both  our  and  the  LWL  results  differ  by 
more  than  two  standard  deviations  from  the  LAH  results,  and 
DEVL-SNRI,  where  the  LWL  and  LAH  results  differ.  Because 
the  GPS  velocities  for  DEVL  and  SNRI.  and  the  LAH  results 
Irom  GPS-EDM  are  all  ba.scd  on  only  two  occupations,  we 
hesitate  to  place  much  significance  on  these  discrepancies.  • 

That  we  estimate  e.sscntially  the  same  rates  in  the  two  different 
inlcmtls;  with  dilferent  techniques,  increa.ses  our  confidence  in 
the  stability  of  the  benchmarks,  as  well  as  in  our  procedures  for 
data  reduction.  It  al.so  suggests  that  the  rate  ol  deformation  acro.ss 
the  channel  has  not  varied  much  over  the  last  20  years. 
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TABLE  A3.  Rales  of  Change  in  Line  Length  for  ihe  Santa  Barbara  Channel  Estimated  from 
GPS-CPS  and  CPS-Trilalcralion  Comparisons 


dUdt, 
mm/vr  - 

Differences, 

mm/vr 

Line 

GPS-GPS‘2 

GPS-CPS* 

CPS-EDM9 

C  -  J 

c-  h 

a  -h 

CENT-CAVl 

DEVL-CAVI 

0.9  ±  1.0 
-L7±  1.1 

2.3*  1.3 
-1.1  ±  1.5 

4.6  *2.1 

3.7  *  2.3 
-  6.3  a  2.4 

2.3  a  2.5 
5.7  a  2.6 

,-1.4  a  1.6 
-  0.6  a  1.9 

■  CENT-LACU 
DEVL-LACU 

-  ^.0  ±  0.6 
-5.6  s  0.9 

-4.0*  1.0 
-5.6  a  1.2 

-6.5  a  1.7 

--i5a  1.8 
-0.9*  1.9 

-2.5  a  2.0 
-0.9  a  2.1 

0.0  a  1.2 

0.0  a  1.5 

CENT-SOLl 

DEVL-SOLI 

-5.1  ±.  1.0 
-4.1  ±  1.3 

^.1  a  1.5 
-2.3  a  1.7 

-5.1  a  1.7  ■ 

0.0  a  2.0 
-1.0  a  2.1 

-1.0  a  2.3 
-2.8  a  2.4 

-l.Oa  1.8 
-1.8  a  2.1 

CENT-SNRI 

DEVL-SNRl 

2.2  ±  1.4 
0.4  ±  1.6 

4.7  a  1.7 
L9a  1.9 

0.4  a  1.1 

-1.8  a  1.8 
0.0*  1.9 

-4.3  a  2.0 
-1.5  a  2.2 

-2.5  a  2.2 
-1.5  a  2.5 

LACU-SOLI 

f.4  =  0.8 

1.6*  1.6 

0.8  a  1.4 

-0.6  a  1.6 

41.8  *2.1 

41.2  a  1.6 

LACU-GAVI 

0.6  a  1.3 

l.Oa  1.2 

0.4  a  1.5 

-0.2  a  2.0 

0.6  a  1.9 

-0.4  a  1.8 

SNRI-GAVI 

-1.4*  1.1 

-LI  a  1.2 

4.3  a 

.  .  5.7  a  2.3 

5.4  a  2:3:-: 

*;-4L3*  1.6  ... 

SNRI-LACU 

-4.0*  1.2 

-.3.1  a  1.3 

^.6  a  2.3 

-0.6  a  2.6 

-1.5  a  2.6 

-0.9  a  1.8 

“  This  study, 

*  iMrson  and  Webb  ( 1 992], 

. — 

. -  ■ 

^  iMrscn  cl  al.  [  1 993|. 


Acknowledgments.  The  global  VLBI  data  used  in  our  analysis  were 
collected  by  the  staffs  of  many  obserealorics  around  the  world  and  made 
available  to  us  primarily  through  the  efforts  of  the  Goddard  Space  Flight 
Center  of  the  National  Aeronautics  and  Space  Administration  (NASA),  Ihe 
National  Geodetic  Survey  (NGS)  and  Geodynamics  Research  Laboratory 
(GDL)  of  National  Oceanic  and  Atmospheric  Administration  (NOAA), 
and  the  U.S.  Naval  Observatory  (USNO).  The  global  GPS  observations 
were  made  available  primarily  through  the  efforts  of  GDL  and  the 
Canadian  Centre  for  Surveying.  The  establishment  in  early  1987  of  our 
core  CPS  network  was  made  possible  by  collaboration  offered  by  Will 
Prescott  of  the  U.S,  Geological  Survey  (USCS)  and  Bill  Strange  of  NGS. 
Most  of  the  Oilifomia  CPS  measurements  were  made  by  the  authors  with 
invaluable  assistance  from  ourcolleague.s,  including  but  not  limited  to  Rob 
Clayton,  Olafur  Gundmunds.son,  Ken  Hudnut,  Dave  Johnson,  Louise 
Kellog,  Scott  King,  Leslie  Sender,  Hong-Kie  Thio,  and  Bruce  Worden 
(Caltech);  John  Goff  and  Peter  Puster  (MIT);  Jeff  Gcnrich,  Glen 
Sasagawa.  Catherine  Hedlin,  Michael  Hcdlin,  Hadley  Johnson,  Nancy 
King  (Scripps);  Abe  Cheng,  Peng  Fang,  Bob  Ce,  Manabu  Hashimoto, 
Naoshi  Hirata,  Mercedes  Kim,  and  Li-Yu  Sung  (UCLA);  James  Stowell 
and  Bruce  Stephens  (U.NAV'CO).  We  thank  the  understanding  landowners, 
including  the  late  Carey  Stanton,  the  Nature  Conservancy,  and  many 
others,  for  allowing  us  to  access  the  bench  marks  on  their  properly.  The 
maps  in  the  figures  of  this  paper  were  generated  using  the  public  domain 
Generic  MappingTools  (GMT)  software  [W'c.v.sc/  and  Smiih.  1991|.  This 
research  was  sponsored  primarily  by  the  National  Science  Foundation 
(grants  EAR-86-18165,  EAR-86- 1 8.39.3.  EAR-86- 1 85 1.3,  EAR-86-1S542. 
E.AR-92-()8447)  with  supplementary  funding  from  the  NASA  Crustal 
Dynamics  Project  (grants  NAG5-4-i7.  NAG5-459.  NAC5-5.38,  NAG5- 
7.37.  NAC5-842,  NAC5-90.5,  NAG.5-1917,  NAGW-21.3.3,  NACW-264I). 
the  U.S.  Geological  Survey  (l4-08-()(t()l-G167.3.  1 4-()8-0(X) I -G 1950, 
I4.34-92-G2185,  14-.34-92-G2196),  the  Southern  California  Earthquake 
Center  (USGS  Cooperative  Agreement  I4-()8-(KK)001-A-()899),  and  the 
.-\ir  Force  Office  of  Scientific  Research  (grants  .^FOSR-89-()40()  and 
.AFOSR-90-0.3.39).  Kurt  FeigI  was  supported  by  an  NSF  Graduate 
Fellowship  and  an  Air  Force  Lriboralory  Graduate  Student  Fellowship. 
Mark  Murray  and  Kristine  Larson  were  supported  in  pan  bv  NASA 
Graduate  Student  Fellowships.  This  is  paper  .38  of  the  Southern  California 
Earthquake  Center. 


Refhrkn('f:s 

Argus,  D.  F.,  and  R;  G.  Gordon.  Pacific-Nonh  American  plalc  motion 
from  very  long  ba.seline  interferometry  compared  with  motion  inferred 
from  magnetic  anomalies,  transform  faults,  and  earthquake  slip  vectors, 

J.  Ccophys.  Res. ,  05,  1 7.3 1 5- 1 7,324,  1 990. 

Argus,  D.  F..  and  R.  G.  Gordon,  Current  Sierra  Nevada-North  America 
motion  from  very  long  baseline  interferometry:  implications  for  the 
kinematics  of  the  we.stern  United  Slates,  Geology,  /O,  1085-1088,  1991. 
Atwater,  T.,  Implications  of  plate  tectonics  for  the  Cenozoic  evolution  of 
western  North  America,  Geol.  Soc.  Am.  Bulk,  SI,  35I.3-.35.36,  1970. 
Bent,  A.  L.,  and  D.  Helmberger,  Seismic  characteristics  of  earthquakes 
along  the  offshore  extension  of  the  Western  Transverse  Ranges, 
California,  Bull.  Scismol.  .Soc.  Am..  81.  399-422,  1991. 

Bird,  P.,  and  R.  W.  Roscnstock,  Kinematics  of  present  crust  and  mantle 
flow  in  southern  California,  Geol.  Soc.  .Am.  Bull.,  95,  946-957,  1984. 
Blcwiti,  G.,  Carrier-phase  ambiguity  resolution  for  the  Global  Positioning 
System  applied  to  geodetic  baselines  up  to  2(K)0  km, 7.  Geophys.  Res., 
94,  l().l87-10.203,'l9S9. 

Blcwiti,  G.,  M,  B.  Heflin.  F.  H.  Webb,  U.  J.  Lindquisler,  and  R.  P.  Malla, 
Global  coordinates  with  centimeter  accuracy  in  the  International 
Terrestrial  Reference  Frame  using  GPS,  Geophys.  Res.  Lett.,  19,  85.3— 
856,  1992. 

BIcw  itt,  G..  M.  B.  Heflin.  K.  J.  Hurst,  D.  C.  Jefferson,  F.  H.  Webb,  and  J. 
F.  Zumberge,  Ab.solule  far-field  displacements  from  the  28  June  1992 
Landers  earthquake  .sequence.  Nature.  .161,  ,340-.342,  1993. 

Bock.  Y.,  Continuous  monitoring  of  crustal  deformation,  GPS  World,  40- 
47,  June,  1991. 

Bock,  Y..  and  S.  Shimada.  Continuously  monitoring  GPS  networks  for 
deformation  measurements,  in  Global  Positioning  System:  .An 
Overview,  Proceedings  of  the  International  .Association  of  Geodesy 
Symposium  102.  edited  by  Y.  Bock  and  N.  Leppard,  pp.  40-56. 

,  Springer-yerlag,  New  York,  1990. 

Bock,  Y.,  S.  A.  Gourevilch.  C.  C.  Counselman  III.  R.  W.  King,  and  R.  I. 
Abbot,  Interferometric  analysis  of  CPS  phase  observations.  Manus. 
GVod.. //,  282-2,88.  19,86. 

BiK'k.  Y..  D.  C.  Agnew,  P.  Fang.  J.  F.  Gcnrich.  B.  H.  Hager,  T.  A,  Herring, 

K.  W.  Hudnut,  R.  W.  King,  S.  Larsen,  J.  B.  Minster.  K.  Stark.  S. 


21,710 


Fcigi.  et  al.:  California  Crustal  Deformation  Measurements 


Wdnowski,  and  F.  K.  Wsall,  Dcicclion  of  cruslal  deformation  from  the 
Landers  eanhquake  using  continuous  geodetic  measuremenis.  Nature, 
36/,  237-340.  1003. 

Bowie,  W.,  Earth  movements  in  California.  Spec.  Puhl.  No.  /06,  U.S. 
Coast  and  Ccod.  Survey.  1024.  (Reprinted  in  //eperts  on  Geodetic 
Measurements  of  Crustal  Movement,  l‘)06-7l  U.S.  National  Geodetic 
Survey.  RcK-kville.  Md..  1073.) 

Bowie,  W.,  Comparison  of  old  and  new  iriangulation  data  In  California. 
V.  S.  Coast  and  Geodetic  .Survey  .Spec.  Publicalioti  No.  151,  1028. 
Reprinted  in  Reports  on  Geodetic  Measurements  of  Crustal  Movement. 
1006-71,  U.  S.  National  Geodetic  Survey,  Rockville.  MD,  107.3. 
Burford.  R.  O..  and  P.  W.  Harsh,  Slip  on  the  San  Andreas  fault  in  central 
California  from  alinement  arravs.  Bull.  Sci.smol.  Soc.  Am.,  70,  12.3.3- 
1261,  1080. 

California  Division  of  Mines  and  Geology.  Preliminary  fault  activity  map 
of  California.  Cilif.  Dep.  of  Conserv..  Sacramento,  1002. 

Callaway.  D.  C..  Petroleum  potential  of  the,  San  Joaquin  Basin,  California, 
in  Future  Petroleum  Provinces  of  the  United  Statc.s:  Their  Geology  and 
Potential,  edited  by  I.  H.  Crum,  vol.  I,  pp.  230-25.3,  American 
.Association  of  Petroleum  Geologists.  Tulsa.  Okla..  1071. 

Cande,  S.  C.,  and  D.  V.  Kent.  A  new  geomagnetic  polarity  time  scale  for 
the  late  Cretaceous  and  Cenozoic../.  Gcophvs.  Re.s.,  07,  13.017-1.3,051, 
1002. 

Carter,  W.  E.,  D.  S.  Robertson,  and  J.  R.  MacKay,  Geodetic  radio 
interferometric  surveying:  Applications  and  results,  J.  Gcophvs.  Res., 
00.4577-4588,1085. 

Cheng,  A.,  D.  D.  Jackson,  and  M.  Matsu’ura,  Aseismic  cruslal  deformation 
in  the  Transverse  Ranges  of  southern  California,  Tcetonophysies,  144, 
150-180,  1087. 

Clark,  M.  .M.,  et  al..  Preliminary  slip-rale  table  and  map  of  late  Quaternary 
faults  of  California,  U.  S.  Gcol.  Surv.  Open  File  Report  i34-l06,  1084. 
Clark,  T.  .A.,  ct  al„  Precision  geodesy  using  the  Mark  III  very-long- 
baseline  interferometry  system,  IEEE  Tran.s.  Gcosei.  Remote  Sens.,  GE- 

2.3,  4.38-440.  1 085. 

Clark,  T.  A..  D.  Cordon,  W.  E.  Himwich,  C.  Ma,  A.  Mallama,  and  J.  W. 
Ryan,  Determinations  of  relative  site  motions  in  the  western  United 
States  using  Mark  III  very  long  baseline  interferometry,  2.  Geophvs. 
/?«.,  02.  12.741-12.750,  1087. 

Clark,  T.  A.,  C.  Ma.  J.  M.  Sauber.  J.  W.  Ryan,  D.  Cordon,  D.  B.  Shaffer, 
D.  S.  Capretle,  and  N.  R.  Vandenberg,  Geodetic  measurement  of 
deformation  in  the  Loma  Priela,  California  earthquake  with  very  long 
baseline  interferometry,  Gcophys.  Res.  Lett.,  17.  2015-2018,  lOOO. 
Coates,  R.  J.,  H.  Frey.  J.  Boswonh,  and  G.  D,  Mead,  Space  age  geodesy: 
The  NASA  crustal  dynamics  project,  IEEE  Trans.  Gcosei.  Remote 
Sens.,  GE-23.  .358-368.  1085. 

Grouch.  J.  K.,  S.  B.  Bachman,  and  J.  T.  Shay,  Post-Miocene  comprcssional 
tectonics  along  the  central  California  margin,  in  Tectonics  and 
.Sedimentation  Along  the  California  Margin,  Puhl.  .?5,  edited  by  J.  K. 
Crouch  and  S.  B.  Bachman,  pp.  .37-54,  Pacific  Section,  Society  of 
Economic  Paleontologists  and  Mineralogists,  Bakersfield,  Calif.,  1084, 
Davis.  J.  C..  Statistics  and  Data  .Analysis  in  Geology,  2nd  cd..  646  pp.. 
John  Wiley,  New  York,  10, S6. 

Da\  is.  J.  L.,  W.  H.  Prescott,  J.  L.  Svarc,  and  K.  J.  Wendt,  Assessment  of 
Global  Positioning  System  measurements  for  studies  of  crustal 
deformation.  J.  Geophvs.  Res.,  04,  1 .3..36.3-1 3,650,  1080. 

Davis,  T.  L..  J.  Namson.  and  R.  F.  Yerkes,  A  cross  section  of  the  Los 
Angeles  area:  Seismically  active  fold  and  thrust  bell,  the  1087  Whittier 
Narrows  eanhquake  and  eanhquake  hazard.2.  Geophvs.  Res..  04.  0644- 
0664,  1080. 

Dehlinger,  P..  and  B.  A.  Boll.  Eanhquakes  and  associated  tectonics  in  a 
part  of  coastal  central  California,  Bull.  Scismol.  Soe.  Am.,  77.  '>056- 

207.3.  1088. 

DeMels,  C.,  R.  G.  Gordon.  D.  F.  Argus,  and  S.  Stein,  Current  plate 
motions,  Geophvs.  J.  Int..  101.  425-478,  1000. 

Dibblee.  T.  W.,  Jr..  The  Rinconada  and  related  faults  in  the  southern  Coast 
Ranges,  California,  and  their  tectonic  significance,  U.S.  Gcdl.  .Surv. 
Prof  Pap..  SI.  1077. 

Dixon,  T.  H.,  .An  introduction  to  the  Global  Positioning  System  and  s(>mc 
geological  applications.  Rev.  Gcophys..  20,  240-276.  1001. 

Dixon.  T.,  G.  Bleu  iit.  K.  Larson,  D.  Agnew,  B.  Hager.  P.  Kroger,  L. 
Skrumeda.  and  W.  Strange,  GPS  measurements  of  regional  deformation 


in  California:  Some  constraints  on  performance.  Eos  Trans.  AGU,  71, 
1051-1056.  1000. 

Dokka.  R.  K.,  and  C.  J.  Travis,  Late  Cenozoic  strike-slip  faulting  in  the 
Mojave  desert.  California,  Tectonics,  0.  31 1-.340,  1 000a. 

Dokka.  R.  K.,  and  C.  J.  Travis,  Role  of  Eastern  California  Shear  Zone  in 
accommodating  Pacific-North  American  plate  motion,  Geophvs.  Res. 
Utt..  17,  I.32.VI.326,  lOOOb. 

Dong.  D,  The  horizontal  velocity  field  in  southern  California  from  a 
combination  of  terrestrial  and  space-geodetic  data,  Ph.D.  thesis,  157  pp.. 
Mass.  Iasi,  of  Technol.,  Cambridge,  1003. 

Dong,  D.  N.,  and  Y.  Bock,  Global  Positioning  System  network  analysis 
with  phase  ambiguity  resolution  applied  to  cruslal  deformation  studies 
in  California, 2.  Gcophys.  Res.,  04,  3040-3066,  1080. 

Donnellan.  ,A..  A  geodetic  study  of  cruslal  deformation  in  the  Ventura 
Basin  region,  southern  California.  Ph.  D.  thesis.  220  pp.,  Calif.  Inst,  of 
Technol.,  Pasadena,  1002. 

Donnellan,  A.,  B.  H.  Hager,  R.  W.  King,  and  Thomas  A.  Herring, 
Geodetic  measurement  of  deformation  in  the  Ventura  Basin  region, 
southern  California.  2.  Geophvs.  Res.,  this  issue. 

Donnellan,  A.,  B.  H.  Hager,  and  R.  W.  King,  Revision  of  geologic  rates 
in  the  Ventura  Basin  from  space  geodetic  measurements,  Nature,  in 
press,  1993. 

Ebcrhart-Phillips,  D..  M.  Lisowski,  and  M.  D.  Zohack,  Crustal  strain  near 
the.  Big  Bend  of  the  San  Andreas  fanit:  AnaTysisrof  the  Los  Padres- 
Tehachapi  trilaleration  networks,  California,  2.  Gcophvs.  Res.,  05, 
11.30-115.3.1900. 

Eksirom.  G.,  and  P.  England.  Seismic  strain  rales  in  regions  of  distributed 
continental  deformation, 2.  Gcophys.  Res.,  04,  10.2.31-10.257.  1980. 
Ellsworth,  W.  L.,  Earthquake  history,  1769-1989,  in  The  San  Andreas 
Fault  System,  edited  by  R.  E.  Wallace,  U.  S.  Ceol.  Surv.  Prof.  Pap. 
ISIS,  1990. 

Feigl.  K.  L.,  Geodetic  measurement  of  tectonic  deformation  in  central 
California,  Ph.D.  thesis,  222  pp..  Mass.  Inst,  of  Technol..  Cambridge, 
1001. 

Feigl,  K.  L.,  R.  W.  King,  and  T.  H.  Jordan,  Geodetic  measurement  of 
tectonic  deformation  in  the  Santa  Maria  fold  and  thrust  bell.  California, 
2.  Geophvs.  Res.,  OS.  2670-2609.  1000. 

Feigl.  K.  L..  R.  W.  King.  T.  A.  Herring,  and  M.  Rothacher,  A  scheme  for 
reducing  the  effect  of  Selective  availability  on  preci.se  CPS  carrier  phase 
measuremenis.  Gcophys.  Re.s.  Lett..  IS,  1289-1292,  1091. 

Cold,  C.  M.,  Automated  contour  mapping  using  triangular  element  data 
.structures  and  an  inlerpolanl  over  each  irregular  domain,  Compui. 
Graphics.  H.  170-175,  1075. 

Hager.  B.  H.,  R.  W,  King,  and  M.  H.  Murray.  Measurement  of  crustal 
deformation  using  the  Global  Positioning  System.  Annu.  Rev.  Earth 
Planet.  Sci.,  10,  35 1-382.  1001. 

Hall.  C.  A.,  Jr.,  Origin  and  development  of  the  Lompoc-Santa  Maria  pull 
apart  basin  and  its  relation  to  the  San  Simeon-Hosgri  strike  slip,  western 
Olifornia,  .Spec.  Rep.  Cohf  Div.  Mines  and  Gcol.,  I.U,  1078. 

Hall.  C.  A.,  Jr.,  San  Luis  Obispo  transform  fault  and  middle  Miocene 
rotation  of  the  western  Transverse  Ranges.  California,  2.  Geophvs.  Re.s.. 
86.1015-1031,1081. 

Harland,  W.  B.,  A.  V,  Cox.  P.  G.  Llewellyn.  C.  A.  C.  Picklon,  A.  G. 
Smith,  and  R.  Walters.  .A  Geologic  Time  Seale,  Cambridge  University 
Press,  New  York,  1 082. 

Harri.s,  R.  A.,  and  P.  Segall.  Detection  of  a  locked  zone  at  depth  on  the 
Parkfield  California  .segment  of  the  San  Andreas  fault.  2.  Geophvs.  Re.s., 
<>2.7045-7062,1087. 

Hauksson.  E.,  and  L.  M.  Jones.  The  July  10,86  Oceanside  earthquake 
.sequence  in  the  continental  borderland,  southern  California,  Bull. 
.Sci.smol.  Soc.  Am..  7S.  1885-1005,  10,88. 

Hauksson,  E.,  el  al..  The  1087  Whittier  Narrows  earthquake  in  the  Los 
.Angeles  metropolitan  area.  .Scicnec.  230,  1400-1412.  1088. 

Herring,  T.  .A..  GLOBK:  Global  Kalman  filter  VLSI  and  GPS  analysis 
-  program  version- 3.4-,  Mass.  Inst-,  of  Technol..  Cambridge.  1 003. 

Herring,  T.  A.,  el  al..  Geodesy  by  radio  interferometry:  Intercontinental 
distance  delcrminatiern  with  subdecimeler  precision. ./.  Geophvs.  Res.. 
86.1647-1651.1081. 

Herring,  T.  A.,  el  at..  Geodesy  by  radio  interferometry:  Ex  idence  for 
contemporary  plate  motion,./.  Gcophys.  Re.s..  01.  8341-8347,  1086. 


Fl  ir.L  KT  AL.;  CAI.HORNIA  CRUSTAL  DtKORMATION  MLASUREMKNTS 


21,711 


Herring,  T.  A.,  J.  L.  Davis,  and  I.  I.  Shapiro,  Geodesy  by  radio  asironomy: 
ihe  application  ol  Kalman  filtering  to  very  long  baseline  inlerltromelry, 
J.  Gcophys.  Res..  <J5,  12.561-12.581,  l‘)QO. 

Herring.  T.  A..  D.  Dong,  and  R.  W.  King,  Sub-milliarcsecond 
determination  of  pole  position  using  Global  Positioning  System  data, 
Giophys.  Res.  I.eil..  IS.  l8d.T-18Q6,  IQQI. 

Hilgen,  F.  J.,  Extension  of  the  astronomically  calibrated  (polarity)  lime 
scale  to  the  Miocene/Pliocene  boundary.  Earth  Planet.  Sci.  Lett.,  107, 
lOdl. 

Hornafius,  J.  S..  B.  P.  Luyendyk,  R.  R.  Terres,  and  M.  J.  Kamerling, 
Timing  and  extent  of  Ncogene  tectonic  rotation  in  the  western 
Transverse  Ranges.  California.  Gcol.  Soc.  Am.  Rail..  07.  1476-1487, 
IW6. 

Hudnul,  K.,  J.  Canias,  L.  Cotton,  D.  Jackson,  R.  Packard,  R.  Smith.  G. 
Tayner.  and  W.  Young.  The  southern  California  inter-county  CPS 
project  (abstract).  Eos  Trans.  AGU,  72  (44),  Fall  Meeting  suppl.,  115. 
|W1. 

Hutton,  L.  K..  L.  M.  Jones,  E.  Hauksson,  and  D.  D.  Given, 
Seismotectonics  of  southern  California,  in  Ncotectonics  of  North 
.America  Decade  Map  Volume  I.  edited  by  D.  B.  Slemmons,  E.  R. 
Engdahl.  M.  D.  Zoback,  and  D.  D.  Blackwell,  pp.  I.TV152,  Geological 
Society  of  America,  Boulder.  Colo.,  IWt. 

Jackson.  J..  and  P.  Molnar,  Active  faulting  and  block  rotations  in  the 
wc.siern  Transverse  Ranges,  California,  J.  Gcophys.  Res.,  05,  2207.T- 
22087,  1990.  (Correction' 7.  Geophys.  Res.,  96,  220.7,  1991.) 

Kanamori.  H.,  and  K,  Satake,  Broadband  study  of  the  1989  Loma  Prieta 
earthquake,  Geophys.  Re.s.  Lett..  17,  1179-1 182,  1990. 

King,  N.  E..  P.  Segall,  and  W.  Prescott,  Geodetic  measurements  near 
Parkfield,  California,  1959-1984,  J.  Geophys.  Res.,  02,  2727-2766, 
1987. 

King,  R.  W.,  and  Y.  Bock,  Documentation  for  the  MIT  GPS  analysis 
software:  GAMIT.  Mass.  Inst,  of  Tcchnol.,  Cambridge,  199.7. 

Larsen,  S.,  Crustal  deformation  in  .southern  California.  Ph.  D.  thesis.  .751 
pp.,  Calif.  Inst,  of  Tcchnol.,  Pasadena,  |99|. 

Larsen.  S.  C.,  D.  C.  Agnew,  and  B.  H.  Hager.  Strain  accumulation  in  the 
Santa  Barbara  Channel:  1970-1988,7.  Gcophys.  Res.,  08,  21 19-21.74, 
199.7. 

Larson,  K.  M..  Precision,  accuracy,  and  tectonics  from  the  Global 
Positioning  System,  Ph.D.  thesis,  269  pp.,  Univ.  of  Calif.,  San  Diego, 
1990a. 

Larson,  K.  M.,  Evaluation  of  CPS  estimates  of  relative  positions  from 
centra!  California.  1986-1988,  Gcophys.  Res.  Lett.,  17,  24.7.7-2436. 
19906. 

Larson,  K.  M..  Application  of  the  Global  Positioning  System  to  crustal 
defomation  measurement,  .7,  Results  from  the  southern  California 
borderlands. 7.  Gcophys.  Res.,  this  issue. 

Larson,  K.  M..  and  D.  C.  Agnew,  Application  of  the  Global  Positioning 
System  to  crustal  deformation  measurement.  I,  Precision  and  accuracy. 
7.  Gcophys.  Res..  96.  16,547-16.565,  1991. 

Larson.  K.  M.,  and  F.  H.  Webb,  Deformation  in  the  Santa  Barbara  Channel 
from  GPS  measurements  1987-1991,  Geophys.  Res.  Lett.,  10,  1491- 
1494,  1992. 

Larson,  K.  M..  F.  H.  Webb,  and  D.  C.  Agnew,  Appliailion  of  the  Global 
Positioning  System  to  crustal  deformation  measurement,  2.  The 
inlluence  of  errors  in  orbit  determination  networks, ./.  Geophys.  Re.s.. 
96,16.567-16.584.1991. 

Legg.  M.  R..  V.  Wong,  and  F.  Saure/..  Geologic  structure  and  tectonics  of 
the  inner  continental  borderland  of  northern  Baja  California,  in  The  Gulf 
and  Peninsular  Province  of  the  Culifornias.  edited  by  J.  P.  Dauphin  and 
B.  R.  T.  Simoneit.  A.APG  Mem.  47.  Tulsa.  Okla..  1989. 

Lerch,  F.  J.,  S.  M.  Klosko.  C.  B.  Patel,  and  C.  A.  Wagner,  A  gravity  model 
for  crustal  dynamics  (GEM-L2),  7.  Geophys.  Res.,‘H),  9.701-0.711. 
1985. 

Li.  V.  C..  and  J.  R.  Rice.  Crustal  deformation  in  great  Californian 
earthquake  cycles. 7,  Geophys.  Res..  92.  1 1.53.7—1 1.551,  1987. 

Lichten.  S.  M..  Estimation  and  filtering  for  high  precision  GPS 
applications.  Manuser.  Geod..  15.  159-176.  1990. 

Lichten,  S.  M..  and  J.  S.  Border.  Strategies  for  high-precision  Global 
Positioning  System  orbit  determination.  7.  Geophys.  Rcs..o2.  12.741- 
12,762.  1987.' 

Lindqwister.  U.  J..  G.  Blewitt.  J.  Zumberge.  and  F.  Webb.  Few  millimeter 


precision  for  baselines  in  the  California  permanent  CPS  geodetic  Array, 
Geophys.  Res.  Lett.,  IS,  1 135-1 138,  1991. 

Lisowski,  M..  W.  H.  Prc.scott,  J.  C.  Savage,  and  M.  L.  Johnston,  Geodetic 
estimate  of  coseismic  slip  during  the  1989  Li'ma  Prieta,  California, 
earthquake,  Geophys.  Res.  Lett.,  17,  1437-1440,  1990. 

Lisowski.  M.,  J.  C.  Savage,  and  W.  H.  Prescott,  The  velocity  field  along 
the  San  Andreas  Fault  in  central  and  .southern  California,  7.  Gcophys. 
Kc.r..  9.7,8.369-8.789,  1991. 

Luyendyk,  B.  P.,  A  model  for  Ncogene  crustal  rotations,  Iranstension,  and 
Iransprcssion  in  southern  California.  Gcol.  Soc.  Amcr.  Bull..  10.1.  1528- 
15.76,  1991.  . 

Malvern, -L.  E.*.  Introduction  to  Mechanics  of  a  Continuous  Medium,  713 
pp..  Prenticc-Hall,  Englewood  Cliffs.  N.J..  1969. 

McDougall.  I.,  F.  H.  Brown.  T.  E.  Cerling,  and  J.  W.  Hillhouse,  A 
reappraisal  of  the  geomagnetic  polarity  time  scale  to  4  Ma  using  data 
from  the  Turkana  Basin.  E;ist  Africa.  Geophys.  Res.  Lett.,  10,  2349- 
2.352,  1992. 

Min.sfer.  J.  B.,  and  T.  H.  Jordan.  Vector  constraints  on  Quaternary 
deformation  of  Ihe  western  United  States  cast  and  west  of  the  San 
Andreas  fault.  Tectonics  and  Sedimentation  Along  the  California 
Margin.  Puhl.  .18,  edited  by  J.  K.  Crouch  and  S.  B.  Bachman,  pp.  1-16. 
Pacific  Section.  Society  of  Economic  Paleontologists  and  Mineralogist.s. 
Bakersfield,  Calif.,  1984. 

Minster,  J.  B.,  and  T.  H.  Jordan,  Vector  constraints  on  western  U.S. 
deformations  from  space  geodesy,  neotectonics  and  'plate  motions,  7.- 
Geophys.  Res.  02,  4798-4804.  1987. 

Mount,  V..  and  J.  Suppe.  Stale  of  slrests  near  the  San  Andreas  fault: 
Implications  for  wrench  tectonics.  Gco/ogv’,  15,  114.3-1146,  1987. 

Murray,  M.  H.,  Global  Positioiitng' System  measurements  of  crustal 
deformation  in  central' California,  Ph.D.  thesis,  310  pp.,  Mass.  Inst,  of 
Tcchnol.,  Cambridge,  1991. 

Namson.  J..  and  T.  Davis.  Structural  transect  of  the  w'csiern  Transverse 
Range.s,  California:  Implications  for  lithospheric  kinematics  and  seismic 
risk  evaluation.  Geology,  16,  675-679.  1988. 

Namson,  J..  and  T.  L.  Davis,  Late  Cenozoic  fold  and  thrust  belt  of  the 
southern  Coasl  Ranges  and  Santa  Maria  Basin,  California,  AAPG  Bull., 
74.467-192.1990. 

Okada.  Y..  Surface  deformation  to  shear  and  tensile  faults  in  a  half-space. 
Bull.  Sci.smol.  .Soc.  Am.,  75,  1 1 35- 1154,  1 985. 

Pacific  Gas  and  Electric  Company,  Final  Report  of  the  Diablo  Canyon 
Long-Term  Safely  Program  for  the  Diablo  Canyon  Power  Plant,  Docket 
numbers  50-275  and  50-323,  submitted  to  the  Nuclear  Regulatory 
Commission,  San  Francisco,  1988. 

Pacheco,  J..  and  J.  Nahelek,  Source  mechanisms  of  three  moderate 
California  earthquakes  of  July  \9S6.  BulL  Scismol.  Soc.  .Am.,  78,  1907- 
1929,  1988. 

Page.  B.  M.,  Geology  of  Ihe  Coasl  Range  of  California,  Bull.  Calif  Die. 
.Mines  Gcol.,  00.  255-276,  1966. 

Page,  B.  .M.,  The  southern  Coasl  Ranges,  in  The  Geotcctonic  Development 
of  California,  edited  by  W.  G.  Ernst,  pp.  329-417,  Prentice  Hall. 
Englewood  Cliffs,  N.J.,  1981. 

Pre.scoil,  W.  H.,  J.  C.  Savage,  and  W.  T.  Kinoshita,  Strain  accumulation  in 
the  western  United  States  between  1970  and  1978, 7.  Geophys.  Res..  84. 
542,V54.75.  1979. 

Ryan.  J.  W..  C.  Ma.'and  D.  S.  Caprett,  N.ASA  Space  Geodesy  Program — 
GSFC  Data  Analysis — 1992:  Final  report  of  the  Crustal  Dynamics 
Project  VLBI  Geodetic  Results  1979-91,  Rep. 104572,  NASA, 
\k'ashington.  DC..  1993. 

Saaslamoinen.  J.,  .Atmospheric  correction  for  the  troposphere  and  the 
stratosphere  in  radio  ranging  of  satellites,  in  The  Use  of  .Artificial 
.Satellites  for  Geodesy,  Geophys.  Mono^r.  Ser..  vol.  15.  edited  by  S.  W. 
Henriksen,  A.  Mancini.  and  B.  H.  Chovilz,  pp.  247-251.  AGU, 
Washington,  D.  C..  1972. 

Sauber,  J.,  Geodetic  measurement  of  deformation  in  California.  Ph.D. 
thesis.  Mass.  Inst.  ofTechnol..  Cambridge.  21 1  pp..  1988. 

Sayagcc  J,_C.,  Equivalent  strike-slip  earthquake  cycles  in  half-space  and 
lilhosphere-a.sihenosphere  Earth  models.  7.  Geophys.  Res..  05.  487.7- 
4879.  199(1, 

Savage.  J.  C.,  and  R.  O.  Burford.  Geodetic  determination  of  relalixe  plate 
motion  in  central  California.  7.  Gcophys.  Res..  78.  832-845,  1973. 

Savage.  J.  C..  and  W.  H.  Prescott.  Geodetic  control  and  the  1927  Lxampoc. 


21,712 


FEICL  ET  AL.:  CALIFORNIA  CRUSTAL  DEFORMATION  MEASUREMENTS 


California  Earthquake.  Hull.  Sci-;mot.  .Soc.  Am.,  (>8,  I69y-!70.1,  1978. 
Savage,  J.  C.,  M.  Lisowski,  and  W.  H.  PreseoII,  An  apparent  .shear  zone 
trending  nonh-northwesi  across  the  Mojave  desert  into  Owens  Valley, 
ea.stern  California,  Gcophys.  Res.  LeU.,  17,  21 1.V21 16,  1990, 

Schaffrin,  B.,  and  Y.  Bock,  A  unified  scheme  for  processing  GPS  phase 
observations.  Bull.  CM.  62,  142-160,  1988. 

Scoiti,  O.,  A,  Nur,  and  R.  Estevez,  Distributed  deformation  and  block 
rotation  in  three  dimensions,  J.  Gcophys.  Res.,  06  P ‘>'>5-1  ■> '’4T 
1991. 

Segall,  P.,  and  M.  V.  Matthews,  Displacement  calculations  from  geodetic 
data  and  the  testing  of  geophysical  deformation  models,  J.  Gcophys 
Res.,02,  14,954-14,966,  1988. 

Shaw,  J.  H.,  and  J.  Suppe,  Slip  rales  on  active  blind  thrust  faults  in  the 
Santa  Barbara  Channel,  Western  Transverse  Ranges  Fold  and  Thrust 
Belt  (abstract),  Eos,  Trans.  Amcr.  Gcophys.  Un.,  72  (44).  Fall  Meeting 
Suppl.,44.7,  1991. 

Shen,  Z.,  Regional  tectonic  deformation  in  southern  California,  inferred 
from  terrestrial  geodesy  and  the  Global  Positioning  System,  Ph.  D. 
thesis,  16.2  pp.,  Univ.  Calif.,  Los  Angeles,  1991. 

Shen,  Z.  K.,  and  D.  D.  Jackson.  Global  Positioning  System  reoccupation  of 
early  trilateration  sites:  tectonic  deformation  of  the  Southern  Coast 
Ranges, 7.  Gcophys.  Res.,  08,  99.21-9946,  199.2. 

Shen,  Z.  K.,  Ge,  B,  X.,  and  D.  D.  Jackson,  Los  Angeles  Basin  Tectonics 
and  the  1991  Sierra  Madre  Earthquake  (abstract),  £os  Trans.  AGU.  72 
'9PPl->  119.  1991. 

Sieh,  K.  E.,  and  R.  H.  Jahns,  Holocene  activity  of  the  San  Andreas  fault  at 
Wallace  Creek.  Gcol.  .8oc.  Am.  Hull.,  OS,  88.2-896,  1984. 

Sieh,  K.,  el  al..  Near-field  investigations  of  the  Landers  earthquake 
sequence.  April  to  July,  \002,  Science,  260,  I7I-I76,  199.2. 

Snay,  R.  A.,  M.  W.  Cline,  and  E.  L.  Timmerman,  Regional  deformation  of 
the  Earth  model  for  the  San  Diego  Region,  California,  7.  Gcophvs  Res 
.VS.  5009-5024,  198.2. 

Snay,  R.  A.,  M.  W.  Cline,  and  E.  L.  Timmerman,  Horizontal  crustal 
deformation  models  for  California  from  historical  geodetic  data,  R.  S  oc 
N.  Z.  Bull.,  24,  1 .2 1  - 1 40.  1986.  ■  ' 

Snay,  R.  A.,  M.  W.  Cline,  and  E.  L.  Timmerman.  Project  REDEAM: 
Models  for  historical  horizontal  deformation,  NO.AA  Tech.  Rep  NOS 
125  NGS-42,m7. 

Stein,  R.  S.,  and  0.  C.  P.  King,  Seismic  potential  revealed  by  surface 
folding:  198.2  Coalinga,  California  earthquake.  Science,  224,  869-872 
1984. 

Stein.  R.  S.,  and  R.  S.  Yeats,  Hidden  earthquakes.  Set.  Am.,  260  (6)  48- 
57,  1989. 

Sylvester,  A.  C.,  and  A.  C.  Darrow,  Structure  and  neolectonics  of  the 
western  Santa  Ynez  fault  system  in  southern  California,  Tecionophvsics 
52.  ,289-}05,  1979. 

Thatcher,  W,,  Nonlinear  strain  buildup  and  the  earthquake  cvclc  on  the 
San  Andreas  faull.7.  Gcophys.  Res.,  88.  589.2-5902,  198.2. 

U.S.  Geologic.ll  Survey  Stalf,  The  Loma  Prieta,  California,  earthquake:  An 
anticipated  event.  Science,  247,  286-29.2,  1990. 


Ward,  S.  N.,  Pacific-North  America  plate  motions:  New  results  from  very 
long  baseline  interferometry,  7.  Gcophys.  Res.,  05.  2 1,965-2 1  981 
1990.  ’  ’ 

Watson,  D.  F.,  ACORD:  Automatic  contouring  of  raw  data,  Compui. 
Gcosci..  8,  97- 1 0 1 .  1 982. 

Webb.  F.  H.,  Geodetic  measurement  of  deformation  in  the  offshore  of 
southern  California,  Ph.D.  thesis,  211  pp.,  Calif.  Inst,  of  Technol.. 
Pasadena,  1991. 

Weldon,  R.,  and  E,_Humphrcy.s,  A  kinematic  model  of  southern  California, 
7cc/o/»ci.  5,  .2.2-48,  1986.  -  .  ' 

Wessel,  P..  and  W.  H.  F.  Smith,  Free  software  helps  map  and  display  data. 

F.os,  Trans.  AGU,  72,  A4\,  \00\.  . 

Whttten,  C.  A^.  Crustal  movement  in  California  and  Nevada,  Eos  'Trans. 
Am.  Gcophys.  Union,  .17,  1956.  (Reprinted  in  Reports  on  Geodetic 
Measurements  of  Crustal  Movement,  1006-71  U.  S.  National  Geodetic 
Survey,  RtK-kville,  Md.,  197.2.) 

Yeats,  R.'S..  Quaternary  flake  tectonics  of  the  California  Transverse 

Ranges, <7co/ogy,  0,  16-20.  1981.  - - 

Yeats,  R.  S..  Large-scale  Quaternary  detachments  in  Ventura  Basin. 

California,  7.  Gcophys.  Res.,  88.  569-58.2,  198.2. 

Yeats,  R.  S.,  and  G.  J.  Huflile,  Alternate  model  for  convergence  acro.ss  the 
Ventura  Basin,  California,  Eos,  Trans.  AGU.  7.?,  (4.2),  Fall  Meeting 
suppl..  590.  1992. 

Yerkes,  R.  F..  .ind  W.  H.„K._(^e,-Late  Quaternary  deformation  in  the 
Western  Transverse  Ranges,  California,-f7.5.^6'c<?/.  .Vurv  Circ  7008  ' 
1979.  .  "  ’■ 

Zoback.  M.  D..  el  al..  New  evidence  on  the  state  of  strc.ss  of  the  San 
Andreas  fault  system.  Science.  288,  1 105-1  Ill.  1987. 


D.  C.  .Agnew  and  Y.  Bock.  Institute  of  Geophysics  and  Planetary 
Physics,  Scripps  Institution  of  Oceanography,  La  Jolla",  CA  9209.2. 

D.  Dong,  A.  Donnellan,  and  F.  H.  Webb,  Jet  Propulsion  Laboratory, 
48(K)  Oak  Grove  Dr.,  Pasadena.  CA  9!  109. 

K.  L.  FeigI,  UPR  2.24,  Observaloire  Midi  Pyrenees,  14  Ave,  E.  Belin, 
.21400  Toulouse.  France. 

B.  H.  Hager,  T.  A.  Herring,  T.  H.  Jordan,  and  R.  W.  King,  Department 
of  Earth,  Atmospheric  and  Planetary  Sciences.  54-620,  Massachusetts 
Institute  of  Technology,  Cambridge.  .MA  021.29. 

D.  D.  Jackson  and  Z.  Shen,  Department  of  Earth  and  Space  Sciences, 
University  of  California,  Los  Angeles,  CA  90024, 

S.  Larsen.  Lawrence  Livermore  National  Laboratory,  L208,  Livermore 
CA  94550. 

K.  M.  Larson,  Department  of  Aerospace  Engineering,  University  of 
Colorado.  Campus  Box  431,  Boulder,  CO  80.209.0431. 

M.  H.  Murray,  US  Geological  Survey,  .245  Middlefield  Rd..  Menlo 
ParkVMS  977.  CA  94025. 

(Received  December  18,  1992; 
revised  August  16,  199.2; 
accepted  August  24,  1 99,2.) 


